Dream: powerful differential expression analysis for repeated measures designs

Personalized Dietary Advice to Increase Protein Intake in Older Adults Does Not Affect the Gut Microbiota, Appetite or Central Processing of Food Stimuli in Community-Dwelling Older Adults: A Six-Month Randomized Controlled Trial

Expert groups argue to raise the recommended daily allowance for protein in older adults from 0.8 to 1.2 g/kg/day to prevent undernutrition. However, protein is thought to increase satiety, possibly through effects on gut microbiota and central appetite regulation. If true, raising daily protein intake may work counterproductively. In a randomized controlled trial, we evaluated the effects of dietary advice aimed at increasing protein intake to 1.2 g/kg adjusted body weight/day (g/kg aBW/day) on appetite and gut microbiota in 90 community-dwelling older adults with habitual protein intake <1.0 g/kg aBW/day (Nintervention = 47, Ncontrol = 43). Food intake was determined by 24-h dietary recalls and gut microbiota by 16S rRNA sequencing. Functional magnetic resonance imaging (fMRI) scans were performed in a subgroup of 48 participants to evaluate central nervous system responses to food-related stimuli. Both groups had mean baseline protein intake of 0.8 ± 0.2 g/kg aBW/day. At 6 months’ follow-up this increased to 1.2 ± 0.2 g/kg aBW/day for the intervention group and 0.9 ± 0.2 g/kg aBW/day for the control group. Microbiota composition was not affected, nor were appetite or brain activity in response to food-related stimuli. Increasing protein intake in older adults to 1.2 g/kg aBW/day does not negatively impact the gut microbiota or suppress appetite.

Cell type-specific histone acetylation profiling of Alzheimer's disease subjects and integration with genetics

We profile genome-wide histone 3 lysine 27 acetylation (H3K27ac) of 3 major brain cell types from hippocampus and dorsolateral prefrontal cortex (dlPFC) of subjects with and without Alzheimer's Disease (AD). We confirm that single nucleotide polymorphisms (SNPs) associated with late onset AD (LOAD) show a strong tendency to reside in microglia-specific gene regulatory elements. Despite this significant colocalization, we find that microglia harbor more acetylation changes associated with age than with amyloid-β (Aβ) load. In contrast, we detect that an oligodendrocyte-enriched glial (OEG) population contains the majority of differentially acetylated peaks associated with Aβ load. These differential peaks reside near both early onset risk genes (APP, PSEN1, PSEN2) and late onset AD risk loci (including BIN1, PICALM, CLU, ADAM10, ADAMTS4, SORL1, FERMT2), Aβ processing genes (BACE1), as well as genes involved in myelinating and oligodendrocyte development processes. Interestingly, a number of LOAD risk loci associated with differentially acetylated risk genes contain H3K27ac peaks that are specifically enriched in OEG. These findings implicate oligodendrocyte gene regulation as a potential mechanism by which early onset and late onset risk genes mediate their effects, and highlight the deregulation of myelinating processes in AD. More broadly, our dataset serves as a resource for the study of functional effects of genetic variants and cell type specific gene regulation in AD.

Molecular states during acute COVID-19 reveal distinct etiologies of long-term sequelae

Post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are debilitating, clinically heterogeneous and of unknown molecular etiology. A transcriptome-wide investigation was performed in 165 acutely infected hospitalized individuals who were followed clinically into the post-acute period. Distinct gene expression signatures of post-acute sequelae were already present in whole blood during acute infection, with innate and adaptive immune cells implicated in different symptoms. Two clusters of sequelae exhibited divergent plasma-cell-associated gene expression patterns. In one cluster, sequelae associated with higher expression of immunoglobulin-related genes in an anti-spike antibody titer-dependent manner. In the other, sequelae associated independently of these titers with lower expression of immunoglobulin-related genes, indicating lower non-specific antibody production in individuals with these sequelae. This relationship between lower total immunoglobulins and sequelae was validated in an external cohort. Altogether, multiple etiologies of post-acute sequelae were already detectable during SARS-CoV-2 infection, directly linking these sequelae with the acute host response to the virus and providing early insights into their development.

Multi-omics identify falling LRRC15 as a COVID-19 severity marker and persistent pro-thrombotic signals in convalescence

Patients with end-stage kidney disease (ESKD) are at high risk of severe COVID-19. Here, we perform longitudinal blood sampling of ESKD haemodialysis patients with COVID-19, collecting samples pre-infection, serially during infection, and after clinical recovery. Using plasma proteomics, and RNA-sequencing and flow cytometry of immune cells, we identify transcriptomic and proteomic signatures of COVID-19 severity, and find distinct temporal molecular profiles in patients with severe disease. Supervised learning reveals that the plasma proteome is a superior indicator of clinical severity than the PBMC transcriptome. We show that a decreasing trajectory of plasma LRRC15, a proposed co-receptor for SARS-CoV-2, is associated with a more severe clinical course. We observe that two months after the acute infection, patients still display dysregulated gene expression related to vascular, platelet and coagulation pathways, including PF4 (platelet factor 4), which may explain the prolonged thrombotic risk following COVID-19.

Longitudinal immune profiling after radiation-attenuated sporozoite vaccination reveals coordinated immune processes correlated with malaria protection

Background

Identifying immune processes required for liver-stage sterilizing immunity to malaria remains an open problem. The IMRAS trial comprised 5x immunizations with radiation-attenuated sporozoites resulting in 55% protection from subsequent challenge.

Methods

To identify correlates of vaccination and protection, we performed detailed systems immunology longitudinal profiling of the entire trial time course including whole blood transcriptomics, detailed PBMC cell phenotyping and serum antigen array profiling of 11 IMRAS radiation-attenuated sporozoite (RAS) vaccinees at up to 21 timepoints each.

Results

RAS vaccination induced serum antibody responses to CSP, TRAP, and AMA1 in all vaccinees. We observed large numbers of differentially expressed genes associated with vaccination response and protection, with distinctly differing transcriptome responses elicited after each immunization. These included inflammatory and proliferative responses, as well as increased abundance of monocyte and DC subsets after each immunization. Increases in Vδ2 γδ; T cells and MAIT cells were observed in response to immunization over the course of study, and CD1c+ CD40+ DC abundance was significantly associated with protection. Interferon responses strongly differed between protected and non-protected individuals with high interferon responses after the 1^st immunization, but not the 2^nd-5^th. Blood transcriptional interferon responses were correlated with abundances of different circulating classical and non-classical monocyte populations.

Conclusions

This study has revealed multiple coordinated immunological processes induced by vaccination and associated with protection. Our work represents the most detailed immunological profiling of a RAS vaccine trial performed to date and will guide the design and interpretation of future malaria vaccine trials.

GlycoEnzOnto: a GlycoEnzyme pathway and molecular function ontology

MOTIVATION

The 'glycoEnzymes' include a set of proteins having related enzymatic, metabolic, transport, structural and cofactor functions. Currently, there is no established ontology to describe glycoEnzyme properties and to relate them to glycan biosynthesis pathways.

RESULTS

We present GlycoEnzOnto, an ontology describing 403 human glycoEnzymes curated along 139 glycosylation pathways, 134 molecular functions and 22 cellular compartments. The pathways described regulate nucleotide-sugar metabolism, glycosyl-substrate/donor transport, glycan biosynthesis and degradation. The role of each enzyme in the glycosylation initiation, elongation/branching and capping/termination phases is described. IUPAC linear strings present systematic human/machine-readable descriptions of individual reaction steps and enable automated knowledge-based curation of biochemical networks. All GlycoEnzOnto knowledge is integrated with the Gene Ontology biological processes. GlycoEnzOnto enables improved transcript overrepresentation analyses and glycosylation pathway identification compared to other available schema, e.g. KEGG and Reactome. Overall, GlycoEnzOnto represents a holistic glycoinformatics resource for systems-level analyses.

AVAILABILITY AND IMPLEMENTATION

https://github.com/neel-lab/GlycoEnzOnto.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational analysis of serum-derived extracellular vesicle miRNAs in juvenile sheep model of single stage Fontan procedure

Patients with single ventricle heart defects requires a series of staged open-heart procedures, termed Fontan palliation. However, while lifesaving, these operations are associated with significant morbidity and early mortality. The attendant complications are thought to arise in response to the abnormal hemodynamics induced by Fontan palliation, although the pathophysiology underlying these physicochemical changes in cardiovascular and other organs remain unknown. Here, we investigated the microRNA (miRNA) content in serum and serum-derived extracellular vesicles (EVs) by sequencing small RNAs from a physiologically relevant sheep model of the Fontan operation. The differential expression analysis identified the enriched miRNA clusters in (1) serum vs. serum-derived EVs and (2) pre-Fontan EVs vs. post-Fontan EVs. Metascape analysis showed that the overexpressed subset of EV miRNAs by Fontan procedure target liver-specific cells, underscoring a potentially important pathway involved in the liver dysfunction that occurs as a consequence of Fontan palliation. We also found that post-Fontan EV miRNAs were associated with senescence and cell death, whereas pre-Fontan EV miRNAs were associated with stem cell maintenance and epithelial-to-mesenchymal transition. This study shows great potential to identify novel circulating EV biomarkers from Fontan sheep serum that may be used for the diagnosis, prognosis, and therapeutics for patients that have undergone Fontan palliation.

lmerSeq: an R package for analyzing transformed RNA-Seq data with linear mixed effects models

Background

Studies that utilize RNA Sequencing (RNA-Seq) in conjunction with designs that introduce dependence between observations (e.g. longitudinal sampling) require specialized analysis tools to accommodate this additional complexity. This R package contains a set of utilities to fit linear mixed effects models to transformed RNA-Seq counts that properly account for this dependence when performing statistical analyses.

Results

In a simulation study comparing lmerSeq and two existing methodologies that also work with transformed RNA-Seq counts, we found that lmerSeq was comprehensively better in terms of nominal error rate control and statistical power.

Conclusions

Existing R packages for analyzing transformed RNA-Seq data with linear mixed models are limited in the variance structures they allow and/or the transformation methods they support. The lmerSeq package offers more flexibility in both of these areas and gave substantially better results in our simulations.

The effect of the menstrual cycle on the circulating microRNA pool in human plasma: a pilot study

STUDY QUESTION

Do ovarian hormone changes influence the levels of cell-free or circulating microRNA (cf-miRNA) across the menstrual cycle?

SUMMARY ANSWER

This exploratory study suggests that fluctuations in hormonal levels throughout the menstrual cycle may alter cf-miRNAs levels.

WHAT IS KNOWN ALREADY

cf-miRNA levels vary with numerous pathological and physiological conditions in both males and females and are regulated by exogenous and endogenous factors, including hormones.

STUDY DESIGN, SIZE, DURATION

A prospective, monocentric study was conducted between March and November 2021. Since this was a pilot study, the sample size was based on feasibility as well as previous similar human studies conducted in different tissues. A total of 20 participants were recruited for the study.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We conducted an exploratory study where blood samples were collected from 16 eumenorrheic females in the early follicular phase, the ovulation phase and the mid-luteal phase of the menstrual cycle. The levels of oestrogen, progesterone, LH and FSH were measured in serum by electrochemiluminescence. The levels of 174 plasma-enriched miRNAs were profiled using a PCR-based panel, including stringent internal and external controls to account for the potential differences in RNA extraction and reverse-transcription stemming from low-RNA input samples.

MAIN RESULTS AND THE ROLE OF CHANCE

This exploratory study suggests that cf-miRNAs may play an active role in the regulation of the female cycle by mediating the expression of genes during fluctuating hormonal changes. Linear mixed-models, adjusted for the relevant variables, showed associations between phases of the menstrual cycle, ovarian hormones and plasma cf-miRNA levels. Validated gene targets of the cf-miRNAs varying with the menstrual cycle were enriched within female reproductive tissues and are primarily involved in cell proliferation and apoptosis.

LARGE SCALE DATA

All relevant data are available from the Mendeley database: LEGER, Bertrand (2022), 'MiRNA and menstrual cycle', Mendeley Data, V1, doi: 10.17632/2br3zp79m3.1.

LIMITATIONS, REASONS FOR CAUTION

Our study was conducted on a small participant cohort. However, it was tightly controlled for endogenous and exogenous confounders, which is critical to ensure robust and reproducible cf-miRNA research. Both adjusted and non-adjusted P-values are presented throughout the article.

WIDER IMPLICATIONS OF THE FINDINGS

Measures of ovarian hormones should be rigorously included in future studies assessing cf-miRNA levels in females and used as time-varying confounders. Our results reinforce the importance of accounting for female-specific biological processes in physiology research by implementing practical or statistical mitigation strategies during data collection and analysis.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the Clinique romande de réadaptation, Sion, Switzerland. S.L. was supported by an Australian Research Council (ARC) Future Fellowship (FT10100278). D.H. was supported by an Executive Dean's Postdoctoral Research Fellowship from Deakin University. The authors declare no competing interests.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis.

SARS-CoV-2 infection produces chronic pulmonary epithelial and immune cell dysfunction with fibrosis in mice

A subset of individuals who recover from coronavirus disease 2019 (COVID-19) develop post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal tissue samples. The mouse-adapted SARS-CoV-2 strain MA10 produces an acute respiratory distress syndrome in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute to clinical recovery phases. At 15 to 120 days after virus clearance, pulmonary histologic findings included subpleural lesions composed of collagen, proliferative fibroblasts, and chronic inflammation, including tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal up-regulation of profibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early antifibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

A translational genomics approach identifies IL10RB as the top candidate gene target for COVID-19 susceptibility

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes that reduce COVID-19 host susceptibility is a critical next step. Using a translational genomics approach that integrates COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX), and perturbagen signatures, we identified IL10RB as the top candidate gene target for COVID-19 host susceptibility. In a series of validation steps, we show that predicted GReX upregulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes and that in vitro IL10RB overexpression is associated with increased viral load and activation of disease-relevant molecular pathways.

Influenza vaccination and single cell multiomics reveal sex dimorphic immune imprints of prior mild COVID-19

Viral infections can have profound and durable functional impacts on the immune system. There is an urgent need to characterize the long-term immune effects of SARS-CoV-2 infection given the persistence of symptoms in some individuals and the continued threat of novel variants. Here we use systems immunology, including longitudinal multimodal single cell analysis (surface proteins, transcriptome, and V(D)J sequences) from 33 previously healthy individuals after recovery from mild, non-hospitalized COVID-19 and 40 age- and sex-matched healthy controls with no history of COVID-19 to comparatively assess the post-infection immune status (mean: 151 days after diagnosis) and subsequent innate and adaptive responses to seasonal influenza vaccination. Identification of both sex-specific and -independent temporally stable changes, including signatures of T-cell activation and repression of innate defense/immune receptor genes (e.g., Toll-like receptors) in monocytes, suggest that mild COVID-19 can establish new post-recovery immunological set-points. COVID-19-recovered males had higher innate, influenza-specific plasmablast, and antibody responses after vaccination compared to healthy males and COVID-19-recovered females, partly attributable to elevated pre-vaccination frequencies of a GPR56 expressing CD8+ T-cell subset in male recoverees that are "poised" to produce higher levels of IFNγ upon inflammatory stimulation. Intriguingly, by day 1 post-vaccination in COVID-19-recovered subjects, the expression of the repressed genes in monocytes increased and moved towards the pre-vaccination baseline of healthy controls, suggesting that the acute inflammation induced by vaccination could partly reset the immune states established by mild COVID-19. Our study reveals sex-dimorphic immune imprints and in vivo functional impacts of mild COVID-19 in humans, suggesting that prior COVID-19, and possibly respiratory viral infections in general, could change future responses to vaccination and in turn, vaccines could help reset the immune system after COVID-19, both in an antigen-agnostic manner.

Ancient dental calculus preserves signatures of biofilm succession and interindividual variation independent of dental pathology

Dental calculus preserves oral microbes, enabling comparative studies of the oral microbiome and health through time. However, small sample sizes and limited dental health metadata have hindered health-focused investigations to date. Here, we investigate the relationship between tobacco pipe smoking and dental calculus microbiomes. Dental calculus from 75 individuals from the 19th century Middenbeemster skeletal collection (Netherlands) were analyzed by metagenomics. Demographic and dental health parameters were systematically recorded, including the presence/number of pipe notches. Comparative data sets from European populations before and after the introduction of tobacco were also analyzed. Calculus species profiles were compared with oral pathology to examine associations between microbiome community, smoking behavior, and oral health status. The Middenbeemster individuals exhibited relatively poor oral health, with a high prevalence of periodontal disease, caries, heavy calculus deposits, and antemortem tooth loss. No associations between pipe notches and dental pathologies, or microbial species composition, were found. Calculus samples before and after the introduction of tobacco showed highly similar species profiles. Observed interindividual microbiome differences were consistent with previously described variation in human populations from the Upper Paleolithic to the present. Dental calculus may not preserve microbial indicators of health and disease status as distinctly as dental plaque.

Cholesterol and matrisome pathways dysregulated in astrocytes and microglia

The impact of apolipoprotein E ε4 (APOE4), the strongest genetic risk factor for Alzheimer's disease (AD), on human brain cellular function remains unclear. Here, we investigated the effects of APOE4 on brain cell types derived from population and isogenic human induced pluripotent stem cells, post-mortem brain, and APOE targeted replacement mice. Population and isogenic models demonstrate that APOE4 local haplotype, rather than a single risk allele, contributes to risk. Global transcriptomic analyses reveal human-specific, APOE4-driven lipid metabolic dysregulation in astrocytes and microglia. APOE4 enhances de novo cholesterol synthesis despite elevated intracellular cholesterol due to lysosomal cholesterol sequestration in astrocytes. Further, matrisome dysregulation is associated with upregulated chemotaxis, glial activation, and lipid biosynthesis in astrocytes co-cultured with neurons, which recapitulates altered astrocyte matrisome signaling in human brain. Thus, APOE4 initiates glia-specific cell and non-cell autonomous dysregulation that may contribute to increased AD risk.

Exposure to the Amino Acids Histidine, Lysine, and Threonine Reduces mTOR Activity and Affects Neurodevelopment in a Human Cerebral Organoid Model

Evidence of the impact of nutrition on human brain development is compelling. Previous in vitro and in vivo results show that three specific amino acids, histidine, lysine, and threonine, synergistically inhibit mTOR activity and behavior. Therefore, the prenatal availability of these amino acids could be important for human neurodevelopment. However, methods to study the underlying mechanisms in a human model of neurodevelopment are limited. Here, we pioneer the use of human cerebral organoids to investigate the impact of amino acid supplementation on neurodevelopment. In this study, cerebral organoids were exposed to 10 mM and 50 mM of the amino acids threonine, histidine, and lysine. The impact was determined by measuring mTOR activity using Western blots, general cerebral organoid size, and gene expression by RNA sequencing. Exposure to threonine, histidine, and lysine led to decreased mTOR activity and markedly reduced organoid size, supporting findings in rodent studies. RNA sequencing identified comprehensive changes in gene expression, with enrichment in genes related to specific biological processes (among which are mTOR signaling and immune function) and to specific cell types, including proliferative precursor cells, microglia, and astrocytes. Altogether, cerebral organoids are responsive to nutritional exposure by increasing specific amino acid concentrations and reflect findings from previous rodent studies. Threonine, histidine, and lysine exposure impacts the early development of human cerebral organoids, illustrated by the inhibition of mTOR activity, reduced size, and altered gene expression.

Comparative computational RNA analysis of cardiac-derived progenitor cells and their extracellular vesicles

Stem/progenitor cells, including cardiac-derived c-kit+ progenitor cells (CPCs), are under clinical evaluation for treatment of cardiac disease. Therapeutic efficacy of cardiac cell therapy can be attributed to paracrine signaling and the release of extracellular vesicles (EVs) carrying diverse cargo molecules. Despite some successes and demonstrated safety, large variation in cell populations and preclinical/clinical outcomes remains a problem. Here, we investigated this variability by sequencing coding and non-coding RNAs of CPCs and CPC-EVs from 30 congenital heart disease patients and used machine learning methods to determine potential mechanistic insights. CPCs retained RNAs related to extracellular matrix organization and exported RNAs related to various signaling pathways to CPC-EVs. CPC-EVs are enriched in miRNA clusters related to cell proliferation and angiogenesis. With network analyses, we identified differences in non-coding RNAs which give insight into age-dependent functionality of CPCs. By taking a quantitative computational approach, we aimed to uncover sources of CPC cell therapy variability.

The transcriptome from asexual to sexual in vitro development of Cystoisospora suis (Apicomplexa: Coccidia)

The apicomplexan parasite Cystoisospora suis is an enteropathogen of suckling piglets with woldwide distribution. As with all coccidian parasites, its lifecycle is characterized by asexual multiplication followed by sexual development with two morphologically distinct cell types that presumably fuse to form a zygote from which the oocyst arises. However, knowledge of the sexual development of C. suis is still limited. To complement previous in vitro studies, we analysed transcriptional profiles at three different time points of development (corresponding to asexual, immature and mature sexual stages) in vitro via RNASeq. Overall, transcription of genes encoding proteins with important roles in gametes biology, oocyst wall biosynthesis, DNA replication and axonema formation as well as proteins with important roles in merozoite biology was identified. A homologue of an oocyst wall tyrosine rich protein of Toxoplasma gondii was expressed in macrogametes and oocysts of C. suis. We evaluated inhibition of sexual development in a host-free culture for C. suis by antiserum specific to this protein to evaluate whether it could be exploited as a candidate for control strategies against C. suis. Based on these data, targets can be defined for future strategies to interrupt parasite transmission during sexual development.

Contribution of Age, Brain Region, Mood Disorder Pathology, and Interindividual Factors on the Methylome of Human Microglia

BACKGROUND

Transcriptome studies have revealed age-, disease-, and region-associated microglial phenotypes reflecting changes in microglial function during development, aging, central nervous system homeostasis, and pathology. The molecular mechanisms that contribute to these transcriptomic changes are largely unknown. The aim of this study was to characterize the DNA methylation landscape of human microglia and the factors that contribute to variations in the microglia methylome. We hypothesized that both age and brain region would have a large impact on DNA methylation in microglia.

METHODS

Microglia from postmortem brain tissue of four different brain regions of 22 donors, encompassing 1 patient with schizophrenia, 13 patients with mood disorder pathology, and 8 control subjects, were isolated and assayed using a genome-wide methylation array.

RESULTS

We found that human microglial cells have a methylation profile distinct from bulk brain tissue and neurons, and age explained a considerable part of the variation. Additionally, we showed that interindividual factors had a much larger effect on the methylation landscape of microglia than brain region, which was also seen at the transcriptome level. In our exploratory analysis, we found various differentially methylated regions that were related to disease status (mood disorder vs. control). This included differentially methylated regions that are linked to gene expression in microglia, as well as to myeloid cell function or neuropsychiatric disorders.

CONCLUSIONS

Although based on relatively small samples, these findings suggest that the methylation profile of microglia is responsive to interindividual variations and thereby plays an important role in the heterogeneity of microglia observed at the transcriptome level.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Evolutionary insights into primate skeletal gene regulation using a comparative cell culture model

The evolution of complex skeletal traits in primates was likely influenced by both genetic and environmental factors. Because skeletal tissues are notoriously challenging to study using functional genomic approaches, they remain poorly characterized even in humans, let alone across multiple species. The challenges involved in obtaining functional genomic data from the skeleton, combined with the difficulty of obtaining such tissues from nonhuman apes, motivated us to consider an alternative in vitro system with which to comparatively study gene regulation in skeletal cell types. Specifically, we differentiated six human (Homo sapiens) and six chimpanzee (Pan troglodytes) induced pluripotent stem cell lines (iPSCs) into mesenchymal stem cells (MSCs) and subsequently into osteogenic cells (bone cells). We validated differentiation using standard methods and collected single-cell RNA sequencing data from over 100,000 cells across multiple samples and replicates at each stage of differentiation. While most genes that we examined display conserved patterns of expression across species, hundreds of genes are differentially expressed (DE) between humans and chimpanzees within and across stages of osteogenic differentiation. Some of these interspecific DE genes show functional enrichments relevant in skeletal tissue trait development. Moreover, topic modeling indicates that interspecific gene programs become more pronounced as cells mature. Overall, we propose that this in vitro model can be used to identify interspecific regulatory differences that may have contributed to skeletal trait differences between species.

Primates display a range of skeletal morphologies and susceptibilities to skeletal diseases, but the molecular basis of these phenotypic differences is unclear. Studies of gene expression variation in primate skeletal tissues are extremely restricted due to the ethical and practical challenges associated with collecting samples. Nevertheless, the ability to study gene regulation in primate skeletal tissues is crucial for understanding how the primate skeleton has evolved. We therefore developed a comparative primate skeletal cell culture model that allows us to access a spectrum of human and chimpanzee cell types as they differentiate from stem cells into bone cells. While most gene expression patterns are conserved across species, we also identified hundreds of differentially expressed genes between humans and chimpanzees within and across stages of differentiation. We also classified cells by osteogenic stage and identified additional interspecific differentially expressed genes which may contribute to skeletal trait differences. We anticipate that this model will be extremely useful for exploring questions related to gene regulation variation in primate bone biology and development.

Placenta and fetal brain share a neurodevelopmental disorder DNA methylation profile in a mouse model of prenatal PCB exposure

Polychlorinated biphenyls (PCBs) are developmental neurotoxicants implicated as environmental risk factors for neurodevelopmental disorders (NDDs). Here, we report the effects of prenatal exposure to a human-relevant mixture of PCBs on the DNA methylation profiles of mouse placenta and fetal brain. Thousands of differentially methylated regions (DMRs) distinguish placenta and fetal brain from PCB-exposed mice from sex-matched vehicle controls. In both placenta and fetal brain, PCB-associated DMRs are enriched for functions related to neurodevelopment and cellular signaling and enriched within regions of bivalent chromatin. The placenta and brain PCB DMRs overlap significantly and map to a shared subset of genes enriched for Wnt signaling, Slit/Robo signaling, and genes differentially expressed in NDD models. The consensus PCB DMRs also significantly overlap with DMRs from human NDD brain and placenta. These results demonstrate that PCB-exposed placenta contains a subset of DMRs that overlap fetal brain DMRs relevant to an NDD.

Large-Scale Meta-Longitudinal Microbiome Data with a Known Batch Factor

Data contamination in meta-approaches where multiple biological samples are combined considerably affects the results of subsequent downstream analyses, such as differential abundance tests comparing multiple groups at a fixed time point. Little has been thoroughly investigated regarding the impact of the lurking variable of various batch sources, such as different days or different laboratories, in more complicated time series experimental designs, for instance, repeatedly measured longitudinal data and metadata. We highlight that the influence of batch factors is significant on subsequent downstream analyses, including longitudinal differential abundance tests, by performing a case study of microbiome time course data with two treatment groups and a simulation study of mimic microbiome longitudinal counts.

A model of persistent post SARS-CoV-2 induced lung disease for target identification and testing of therapeutic strategies

COVID-19 survivors develop post-acute sequelae of SARS-CoV-2 (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal samples. Mouse-adapted SARS-CoV-2 MA10 produces an acute respiratory distress syndrome (ARDS) in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute disease through clinical recovery. At 15-120 days post-virus clearance, histologic evaluation identified subpleural lesions containing collagen, proliferative fibroblasts, and chronic inflammation with tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal upregulation of pro-fibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early anti-fibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

Gene essentiality in cancer cell lines is modified by the sex chromosomes

Human sex differences arise from gonadal hormones and sex chromosomes. Studying the direct effects of sex chromosomes in humans is still challenging. Here we studied how the sex chromosomes can modulate gene expression and the outcome of mutations across the genome by exploiting the tendency of cancer cell lines to lose or gain sex chromosomes. We inferred the dosage of the sex chromosomes in 355 female and 408 male cancer cell lines and used it to dissect the contributions of the Y and X Chromosomes to sex-biased gene expression. Furthermore, based on genome-wide CRISPR screens, we identified genes whose essentiality is different between male and female cells depending on the sex chromosomes. The most significant genes were X-linked genes compensated by Y-linked paralogs. Our sex-based analysis identifies genes that, when mutated, can affect male and female cells differently and reinforces the roles of the X and Y Chromosomes in sex-specific cell function.

Genetic analysis of the human microglial transcriptome across brain regions, aging and disease pathologies

Microglia have emerged as important players in brain aging and pathology. To understand how genetic risk for neurological and psychiatric disorders is related to microglial function, large transcriptome studies are essential. Here we describe the transcriptome analysis of 255 primary human microglial samples isolated at autopsy from multiple brain regions of 100 individuals. We performed systematic analyses to investigate various aspects of microglial heterogeneities, including brain region and aging. We mapped expression and splicing quantitative trait loci and showed that many neurological disease susceptibility loci are mediated through gene expression or splicing in microglia. Fine-mapping of these loci nominated candidate causal variants that are within microglia-specific enhancers, finding associations with microglial expression of USP6NL for Alzheimer's disease and P2RY12 for Parkinson's disease. We have built the most comprehensive catalog to date of genetic effects on the microglial transcriptome and propose candidate functional variants in neurological and psychiatric disorders.

Sex Differences in the Human Brain Transcriptome of Cases With Schizophrenia

BACKGROUND

While schizophrenia differs between males and females in the age of onset, symptomatology, and disease course, the molecular mechanisms underlying these differences remain uncharacterized.

METHODS

To address questions about the sex-specific effects of schizophrenia, we performed a large-scale transcriptome analysis of RNA sequencing data from 437 controls and 341 cases from two distinct cohorts from the CommonMind Consortium.

RESULTS

Analysis across the cohorts identified a reproducible gene expression signature of schizophrenia that was highly concordant with previous work. Differential expression across sex was reproducible across cohorts and identified X- and Y-linked genes, as well as those involved in dosage compensation. Intriguingly, the sex expression signature was also enriched for genes involved in neurexin family protein binding and synaptic organization. Differential expression analysis testing a sex-by-diagnosis interaction effect did not identify any genome-wide signature after multiple testing corrections. Gene coexpression network analysis was performed to reduce dimensionality from thousands of genes to dozens of modules and elucidate interactions among genes. We found enrichment of coexpression modules for sex-by-diagnosis differential expression signatures, which were highly reproducible across the two cohorts and involved a number of diverse pathways, including neural nucleus development, neuron projection morphogenesis, and regulation of neural precursor cell proliferation.

CONCLUSIONS

Overall, our results indicate that the effect size of sex differences in schizophrenia gene expression signatures is small and underscore the challenge of identifying robust sex-by-diagnosis signatures, which will require future analyses in larger cohorts.

Comparison of approaches to transcriptomic analysis in multi-sampled tumors

Intratumoral heterogeneity is a well-documented feature of human cancers and is associated with outcome and treatment resistance. However, a heterogeneous tumor transcriptome contributes an unknown level of variability to analyses of differentially expressed genes (DEGs) that may contribute to phenotypes of interest, including treatment response. Although current clinical practice and the vast majority of research studies use a single sample from each patient, decreasing costs of sequencing technologies and computing power have made repeated-measures analyses increasingly economical. Repeatedly sampling the same tumor increases the statistical power of DEG analysis, which is indispensable toward downstream analysis and also increases one's understanding of within-tumor variance, which may affect conclusions. Here, we compared five different methods for analyzing gene expression profiles derived from repeated sampling of human prostate tumors in two separate cohorts of patients. We also benchmarked the sensitivity of generalized linear models to linear mixed models for identifying DEGs contributing to relevant prostate cancer pathways based on a ground-truth model.

A Markov random field model for network-based differential expression analysis of single-cell RNA-seq data

BACKGROUND

Recent development of single cell sequencing technologies has made it possible to identify genes with different expression (DE) levels at the cell type level between different groups of samples. In this article, we propose to borrow information through known biological networks to increase statistical power to identify differentially expressed genes (DEGs).

RESULTS

We develop MRFscRNAseq, which is based on a Markov random field (MRF) model to appropriately accommodate gene network information as well as dependencies among cell types to identify cell-type specific DEGs. We implement an Expectation-Maximization (EM) algorithm with mean field-like approximation to estimate model parameters and a Gibbs sampler to infer DE status. Simulation study shows that our method has better power to detect cell-type specific DEGs than conventional methods while appropriately controlling type I error rate. The usefulness of our method is demonstrated through its application to study the pathogenesis and biological processes of idiopathic pulmonary fibrosis (IPF) using a single-cell RNA-sequencing (scRNA-seq) data set, which contains 18,150 protein-coding genes across 38 cell types on lung tissues from 32 IPF patients and 28 normal controls.

CONCLUSIONS

The proposed MRF model is implemented in the R package MRFscRNAseq available on GitHub. By utilizing gene-gene and cell-cell networks, our method increases statistical power to detect differentially expressed genes from scRNA-seq data.

Acute COVID-19 gene-expression profiles show multiple etiologies of long-term sequelae

Two years into the SARS-CoV-2 pandemic, the post-acute sequelae of infection are compounding the global health crisis. Often debilitating, these sequelae are clinically heterogeneous and of unknown molecular etiology. Here, a transcriptome-wide investigation of this new condition was performed in a large cohort of acutely infected patients followed clinically into the post-acute period. Gene expression signatures of post-acute sequelae were already present in whole blood during the acute phase of infection, with both innate and adaptive immune cells involved. Plasma cells stood out as driving at least two distinct clusters of sequelae, one largely dependent on circulating antibodies against the SARS-CoV-2 spike protein and the other antibody-independent. Altogether, multiple etiologies of post-acute sequelae were found concomitant with SARS-CoV-2 infection, directly linking the emergence of these sequelae with the host response to the virus.

Novel insights into the molecular mechanisms underlying risk of colorectal cancer from smoking and red/processed meat carcinogens by modeling exposure in normal colon organoids

Tobacco smoke and red/processed meats are well-known risk factors for colorectal cancer (CRC). Most research has focused on studies of normal colon biopsies in epidemiologic studies or treatment of CRC cell lines in vitro. These studies are often constrained by challenges with accuracy of self-report data or, in the case of CRC cell lines, small sample sizes and lack of relationship to normal tissue at risk. In an attempt to address some of these limitations, we performed a 24-hour treatment of a representative carcinogens cocktail in 37 independent organoid lines derived from normal colon biopsies. Machine learning algorithms were applied to bulk RNA-sequencing and revealed cellular composition changes in colon organoids. We identified 738 differentially expressed genes in response to carcinogens exposure. Network analysis identified significantly different modules of co-expression, that included genes related to MSI-H tumor biology, and genes previously implicated in CRC through genome-wide association studies. Our study helps to better define the molecular effects of representative carcinogens from smoking and red/processed meat in normal colon epithelial cells and in the etiology of the MSI-H subtype of CRC, and suggests an overlap between molecular mechanisms involved in inherited and environmental CRC risk.

Dysregulation of mitochondrial and proteolysosomal genes in Parkinson's disease myeloid cells

An increasing number of identified Parkinson's disease (PD) risk loci contain genes highly expressed in innate immune cells, yet their role in pathology is not understood. We hypothesize that PD susceptibility genes modulate disease risk by influencing gene expression within immune cells. To address this, we have generated transcriptomic profiles of monocytes from 230 individuals with sporadic PD and healthy subjects. We observed a dysregulation of mitochondrial and proteasomal pathways. We also generated transcriptomic profiles of primary microglia from brains of 55 subjects and observed discordant transcriptomic signatures of mitochondrial genes in PD monocytes and microglia. We further identified 17 PD susceptibility genes whose expression, relative to each risk allele, is altered in monocytes. These findings reveal widespread transcriptomic alterations in PD monocytes, with some being distinct from microglia, and facilitate efforts to understand the roles of myeloid cells in PD as well as the development of biomarkers.

ELAVL4, splicing, and glutamatergic dysfunction precede neuron loss in MAPT mutation cerebral organoids

Frontotemporal dementia (FTD) because of MAPT mutation causes pathological accumulation of tau and glutamatergic cortical neuronal death by unknown mechanisms. We used human induced pluripotent stem cell (iPSC)-derived cerebral organoids expressing tau-V337M and isogenic corrected controls to discover early alterations because of the mutation that precede neurodegeneration. At 2 months, mutant organoids show upregulated expression of MAPT, glutamatergic signaling pathways, and regulators, including the RNA-binding protein ELAVL4, and increased stress granules. Over the following 4 months, mutant organoids accumulate splicing changes, disruption of autophagy function, and build-up of tau and P-tau-S396. By 6 months, tau-V337M organoids show specific loss of glutamatergic neurons as seen in individuals with FTD. Mutant neurons are susceptible to glutamate toxicity, which can be rescued pharmacologically by the PIKFYVE kinase inhibitor apilimod. Our results demonstrate a sequence of events that precede neurodegeneration, revealing molecular pathways associated with glutamate signaling as potential targets for therapeutic intervention in FTD.

Characterization of the SARS-CoV-2 Host Response in Primary Human Airway Epithelial Cells from Aged Individuals

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19), a global pandemic characterized by an exaggerated immune response and respiratory illness. Age (>60 years) is a significant risk factor for developing severe COVID-19. To better understand the host response of the aged airway epithelium to SARS-CoV-2 infection, we performed an in vitro study using primary human bronchial epithelial cells from donors >67 years of age differentiated on an air–liquid interface culture. We demonstrate that SARS-CoV-2 infection leads to early induction of a proinflammatory response and a delayed interferon response. In addition, we observed changes in the genes and pathways associated with cell death and senescence throughout infection. In summary, our study provides new and important insights into the temporal kinetics of the airway epithelial innate immune response to SARS-CoV-2 in older individuals.

Downregulation of exhausted cytotoxic T cells in gene expression networks of multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56dimCD57+ natural killer (NK) cells and exhausted CD8+ T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8+ T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C.

Single-nucleus transcriptome analysis of human brain immune response in patients with severe COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has been associated with neurological and neuropsychiatric illness in many individuals. We sought to further our understanding of the relationship between brain tropism, neuro-inflammation, and host immune response in acute COVID-19 cases.

METHODS

Three brain regions (dorsolateral prefrontal cortex, medulla oblongata, and choroid plexus) from 5 patients with severe COVID-19 and 4 controls were examined. The presence of the virus was assessed by western blot against viral spike protein, as well as viral transcriptome analysis covering > 99% of SARS-CoV-2 genome and all potential serotypes. Droplet-based single-nucleus RNA sequencing (snRNA-seq) was performed in the same samples to examine the impact of COVID-19 on transcription in individual cells of the brain.

RESULTS

Quantification of viral spike S1 protein and viral transcripts did not detect SARS-CoV-2 in the postmortem brain tissue. However, analysis of 68,557 single-nucleus transcriptomes from three distinct regions of the brain identified an increased proportion of stromal cells, monocytes, and macrophages in the choroid plexus of COVID-19 patients. Furthermore, differential gene expression, pseudo-temporal trajectory, and gene regulatory network analyses revealed transcriptional changes in the cortical microglia associated with a range of biological processes, including cellular activation, mobility, and phagocytosis.

CONCLUSIONS

Despite the absence of detectable SARS-CoV-2 in the brain at the time of death, the findings suggest significant and persistent neuroinflammation in patients with acute COVID-19.

Loss of coordinated expression between ribosomal and mitochondrial genes revealed by comprehensive characterization of a large family with a rare Mendelian disorder

Non-canonical intronic variants are a poorly characterized yet highly prevalent class of alterations associated with Mendelian disorders. Here, we report the first RNA expression and splicing analysis from a family whose members carry a non-canonical splice variant in an intron of RPL11 (c.396 +3A>G). This mutation is causative for Diamond Blackfan Anemia (DBA) in this family despite incomplete penetrance and variable expressivity. Our analyses revealed a complex pattern of disruptions with many novel junctions of RPL11. These include an RPL11 transcript that is translated with a late stop codon in the 3' untranslated region (3'UTR) of the main isoform. We observed that RPL11 transcript abundance is comparable among carriers regardless of symptom severity. Interestingly, both the small and large ribosomal subunit transcripts were significantly overexpressed in individuals with a history of anemia in addition to congenital abnormalities. Finally, we discovered that coordinated expression between mitochondrial components and RPL11 was lost in all carriers, which may lead to variable expressivity. Overall, this study highlights the importance of RNA splicing and expression analyses in families for molecular characterization of Mendelian diseases.

IL10RB as a key regulator of COVID-19 host susceptibility and severity

BACKGROUND

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes and readily available compounds that reduce COVID-19 host susceptibility is a critical next step.

METHODS

We integrate COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX) and perturbargen signatures to identify candidate genes and compounds that reverse the predicted gene expression dysregulation associated with COVID-19 susceptibility. The top candidate gene is validated by testing both its GReX and observed blood transcriptome association with COVID-19 severity, as well as by in vitro perturbation to quantify effects on viral load and molecular pathway dysregulation. We validate the in silico drug repositioning analysis by examining whether the top candidate compounds decrease COVID-19 incidence based on epidemiological evidence.

RESULTS

We identify IL10RB as the top key regulator of COVID-19 host susceptibility. Predicted GReX up-regulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes. In vitro IL10RB overexpression is associated with increased viral load and activation of immune-related molecular pathways. Azathioprine and retinol are prioritized as candidate compounds to reduce the likelihood of testing positive for COVID-19.

CONCLUSIONS

We establish an integrative data-driven approach for gene target prioritization. We identify and validate IL10RB as a suitable molecular target for modulation of COVID-19 host susceptibility. Finally, we provide evidence for a few readily available medications that would warrant further investigation as drug repositioning candidates.

A large-scale genome-wide gene expression analysis in peripheral blood identifies very few differentially expressed genes related to antidepressant treatment and response in patients with major depressive disorder

A better understanding of the biological factors underlying antidepressant treatment in patients with major depressive disorder (MDD) is needed. We perform gene expression analyses and explore sources of variability in peripheral blood related to antidepressant treatment and treatment response in patients suffering from recurrent MDD at baseline and after 8 weeks of treatment. The study includes 281 patients, which were randomized to 8 weeks of treatment with vortioxetine (N = 184) or placebo (N = 97). To our knowledge, this is the largest dataset including both gene expression in blood and placebo-controlled treatment response measured by a clinical scale in a randomized clinical trial. We identified three novel genes whose RNA expression levels at baseline and week 8 are significantly (FDR < 0.05) associated with treatment response after 8 weeks of treatment. Among these genes were SOCS3 (FDR = 0.0039) and PROK2 (FDR = 0.0028), which have previously both been linked to depression. Downregulation of these genes was associated with poorer treatment response. We did not identify any genes that were differentially expressed between placebo and vortioxetine groups at week 8 or between baseline and week 8 of treatment. Nor did we replicate any genes identified in previous peripheral blood gene expression studies examining treatment response. Analysis of genome-wide expression variability showed that type of treatment and treatment response explains very little of the variance, a median of <0.0001% and 0.05% in gene expression across all genes, respectively. Given the relatively large size of the study, the limited findings suggest that peripheral blood gene expression might not be the best approach to explore the biological factors underlying antidepressant treatment.

A streamlined whole blood CyTOF workflow defines a circulating immune cell signature of COVID-19

Mass cytometry (CyTOF) represents one of the most powerful tools in immune phenotyping, allowing high throughput quantification of over 40 parameters at single-cell resolution. However, wide deployment of CyTOF-based immune phenotyping studies are limited by complex experimental workflows and the need for specialized CyTOF equipment and technical expertise. Furthermore, differences in cell isolation and enrichment protocols, antibody reagent preparation, sample staining, and data acquisition protocols can all introduce technical variation that can confound integrative analyses of large data-sets of samples processed across multiple labs. Here, we present a streamlined whole blood CyTOF workflow which addresses many of these sources of experimental variation and facilitates wider adoption of CyTOF immune monitoring across sites with limited technical expertise or sample-processing resources or equipment. Our workflow utilizes commercially available reagents including the Fluidigm MaxPar Direct Immune Profiling Assay (MDIPA), a dry tube 30-marker immunophenotyping panel, and SmartTube Proteomic Stabilizer, which allows for simple and reliable fixation and cryopreservation of whole blood samples. We validate a workflow that allows for streamlined staining of whole blood samples with minimal processing requirements or expertise at the site of sample collection, followed by shipment to a central CyTOF core facility for batched downstream processing and data acquisition. We apply this workflow to characterize 184 whole blood samples collected longitudinally from a cohort of 72 hospitalized COVID-19 patients and healthy controls, highlighting dynamic disease-associated changes in circulating immune cell frequency and phenotype.

Temporal Dynamic Methods for Bulk RNA-Seq Time Series Data

Dynamic studies in time course experimental designs and clinical approaches have been widely used by the biomedical community. These applications are particularly relevant in stimuli-response models under environmental conditions, characterization of gradient biological processes in developmental biology, identification of therapeutic effects in clinical trials, disease progressive models, cell-cycle, and circadian periodicity. Despite their feasibility and popularity, sophisticated dynamic methods that are well validated in large-scale comparative studies, in terms of statistical and computational rigor, are less benchmarked, comparing to their static counterparts. To date, a number of novel methods in bulk RNA-Seq data have been developed for the various time-dependent stimuli, circadian rhythms, cell-lineage in differentiation, and disease progression. Here, we comprehensively review a key set of representative dynamic strategies and discuss current issues associated with the detection of dynamically changing genes. We also provide recommendations for future directions for studying non-periodical, periodical time course data, and meta-dynamic datasets.

muscat detects subpopulation-specific state transitions from multi-sample multi-condition single-cell transcriptomics data

Single-cell RNA sequencing (scRNA-seq) has become an empowering technology to profile the transcriptomes of individual cells on a large scale. Early analyses of differential expression have aimed at identifying differences between subpopulations to identify subpopulation markers. More generally, such methods compare expression levels across sets of cells, thus leading to cross-condition analyses. Given the emergence of replicated multi-condition scRNA-seq datasets, an area of increasing focus is making sample-level inferences, termed here as differential state analysis; however, it is not clear which statistical framework best handles this situation. Here, we surveyed methods to perform cross-condition differential state analyses, including cell-level mixed models and methods based on aggregated pseudobulk data. To evaluate method performance, we developed a flexible simulation that mimics multi-sample scRNA-seq data. We analyzed scRNA-seq data from mouse cortex cells to uncover subpopulation-specific responses to lipopolysaccharide treatment, and provide robust tools for multi-condition analysis within the muscat R package.

Cytotoxic lymphocytes are dysregulated in multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8 + T-cells and CD56 dim CD57 + NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21 , a central coordinator of exhausted CD8 + T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.