Data-driven hypothesis weighting increases detection power in genome-scale multiple testing

BRD9 determines the cell fate of hematopoietic stem cells by regulating chromatin state

ATP-dependent chromatin remodeling SWI/SNF complexes exist in three subcomplexes: canonical BAF (cBAF), polybromo BAF (PBAF), and a newly described non-canonical BAF (ncBAF). While cBAF and PBAF regulate fates of multiple cell types, roles for ncBAF in hematopoietic stem cells (HSCs) have not been investigated. Motivated by recent discovery of disrupted expression of BRD9, an essential component of ncBAF, in multiple cancers, including clonal hematopoietic disorders, we evaluate here the role of BRD9 in normal and malignant HSCs. BRD9 loss enhances chromatin accessibility, promoting myeloid lineage skewing while impairing B cell development. BRD9 significantly colocalizes with CTCF, whose chromatin recruitment is augmented by BRD9 loss, leading to altered chromatin state and expression of myeloid-related genes within intact topologically associating domains. These data uncover ncBAF as critical for cell fate specification in HSCs via three-dimensional regulation of gene expression and illuminate roles for ncBAF in normal and malignant hematopoiesis.

Quorum sensing modulates bacterial virulence and colonization dynamics of the gastrointestinal pathogen Citrobacter rodentium

Quorum Sensing (QS) is a form of cell-to-cell communication that enables bacteria to modify behavior according to their population density. While QS has been proposed as a potential intervention against pathogen infection, QS-mediated communication within the mammalian digestive tract remains understudied. Using an LC-MS/MS approach, we discovered that Citrobacter rodentium, a natural murine pathogen used to model human infection by pathogenic Escherichia coli, utilizes the CroIR system to produce three QS-molecules. We then profiled their accumulation both in vitro and across different gastrointestinal sites over the course of infection. Importantly, we found that in the absence of QS capabilities the virulence of C. rodentium is enhanced. This highlights the role of QS as an effective mechanism to regulate virulence according to the pathogen's spatio-temporal context to optimize colonization and transmission success. These results also demonstrate that inhibiting QS may not always be an effective strategy for the control of virulence.

Detection of oxaliplatin- and cisplatin-DNA lesions requires different global genome repair mechanisms that affect their clinical efficacy

The therapeutic efficacy of cisplatin and oxaliplatin depends on the balance between the DNA damage induction and the DNA damage response of tumor cells. Based on clinical evidence, oxaliplatin is administered to cisplatin-unresponsive cancers, but the underlying molecular causes for this tumor specificity are not clear. Hence, stratification of patients based on DNA repair profiling is not sufficiently utilized for treatment selection. Using a combination of genetic, transcriptomics and imaging approaches, we identified factors that promote global genome nucleotide excision repair (GG-NER) of DNA-platinum adducts induced by oxaliplatin, but not by cisplatin. We show that oxaliplatin-DNA lesions are a poor substrate for GG-NER initiating factor XPC and that DDB2 and HMGA2 are required for efficient binding of XPC to oxaliplatin lesions and subsequent GG-NER initiation. Loss of DDB2 and HMGA2 therefore leads to hypersensitivity to oxaliplatin but not to cisplatin. As a result, low DDB2 levels in different colon cancer cells are associated with GG-NER deficiency and oxaliplatin hypersensitivity. Finally, we show that colon cancer patients with low DDB2 levels have a better prognosis after oxaliplatin treatment than patients with high DDB2 expression. We therefore propose that DDB2 is a promising predictive marker of oxaliplatin treatment efficiency in colon cancer.

Inference on tree-structured subgroups with subgroup size and subgroup effect relationship in clinical trials

When multiple candidate subgroups are considered in clinical trials, we often need to make statistical inference on the subgroups simultaneously. Classical multiple testing procedures might not lead to an interpretable and efficient inference on the subgroups as they often fail to take subgroup size and subgroup effect relationship into account. In this paper, built on the selective traversed accumulation rules (STAR), we propose a data-adaptive and interactive multiple testing procedure for subgroups which can take subgroup size and subgroup effect relationship into account under prespecified tree structure. The proposed method is easy-to-implement and can lead to a more interpretable and efficient inference on prespecified tree-structured subgroups. Possible accommodations to post hoc identified tree-structure subgroups are also discussed in the paper. We demonstrate the merit of our proposed method by re-analyzing the panitumumab trial with the proposed method.

Unbiased discovery of natural sequence variants that influence fungal virulence

Isolates of Cryptococcus neoformans, a fungal pathogen that kills over 112,000 people each year, differ from a 19-Mb reference genome at a few thousand up to almost a million DNA sequence positions. We used bulked segregant analysis and association analysis, genetic methods that require no prior knowledge of sequence function, to address the key question of which naturally occurring sequence variants influence fungal virulence. We identified a region containing such variants, prioritized them, and engineered strains to test our findings in a mouse model of infection. At one locus, we identified a 4-nt variant in the PDE2 gene that occurs in common laboratory strains and severely truncates the encoded phosphodiesterase. The resulting loss of phosphodiesterase activity significantly impacts virulence. Our studies demonstrate a powerful and unbiased strategy for identifying key genomic regions in the absence of prior information and provide significant sequence and strain resources to the community.

Systemic and intrinsic functions of ATRX in glial cell fate and CNS myelination in male mice

Myelin, an extension of the oligodendrocyte plasma membrane, wraps around axons to facilitate nerve conduction. Myelination is compromised in ATR-X intellectual disability syndrome patients, but the causes are unknown. We show that loss of ATRX leads to myelination deficits in male mice that are partially rectified upon systemic thyroxine administration. Targeted ATRX inactivation in either neurons or oligodendrocyte progenitor cells (OPCs) reveals OPC-intrinsic effects on myelination. OPCs lacking ATRX fail to differentiate along the oligodendrocyte lineage and acquire a more plastic state that favors astrocytic differentiation in vitro and in vivo. ATRX chromatin occupancy in OPCs greatly overlaps with that of the chromatin remodelers CHD7 and CHD8 as well as H3K27Ac, a mark of active enhancers. Overall, our data indicate that ATRX regulates the onset of myelination systemically via thyroxine, and by promoting OPC differentiation and suppressing astrogliogenesis. These functions of ATRX identified in mice could explain white matter pathogenesis observed in ATR-X syndrome patients.

Integrated Proteomics Identifies Troponin I Isoform Switch as a Regulator of a Sarcomere-Metabolism Axis During Cardiac Regeneration

Adult mammalian cardiomyocytes have limited proliferative potential, and after myocardial infarction (MI), injured cardiac tissue is replaced with fibrotic scar rather than with functioning myocardium. In contrast, the neonatal mouse heart possesses a regenerative capacity governed by cardiomyocyte proliferation; however, a metabolic switch from glycolysis to fatty acid oxidation during postnatal development results in loss of this regenerative capacity. Interestingly, a sarcomere isoform switch also takes place during postnatal development where slow skeletal troponin I (ssTnI) is replaced with cardiac troponin I (cTnI). In this study, we first employ integrated quantitative bottom-up and top-down proteomics to comprehensively define the proteomic and sarcomeric landscape during postnatal heart maturation. Utilizing a cardiomyocyte-specific ssTnI transgenic mouse model, we found that ssTnI overexpression increased cardiomyocyte proliferation and the cardiac regenerative capacity of the postnatal heart following MI compared to control mice by histological analysis. Our global proteomic analysis of ssTnI transgenic mice following MI reveals that ssTnI overexpression induces a significant shift in the cardiac proteomic landscape. This shift is characterized by an upregulation of key proteins involved in glycolytic metabolism. Collectively, our data suggest that the postnatal TnI isoform switch may play a role in the metabolic shift from glycolysis to fatty acid oxidation during postnatal maturation. This underscores the significance of a sarcomere-metabolism axis during cardiomyocyte proliferation and heart regeneration.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

IL-2-induced Stat3 Signaling is Critical for Effector Treg Cell Programming

Maintenance of immune homeostasis to the intestinal mictrobiota is dependent on a population of effector regulatory T (eTreg) cells that develop from microbiota-reactive induced (i)Treg cells. A cardinal feature of eTreg cells is their production of IL-10, which plays a non-redundant role in immune tolerance of commensal microbes. Here, we identify an unexpected role for IL-2-induced Stat3 signaling to program iTreg cells for eTreg cell differentiation and Il10 transcriptional competency. IL-2 proved to be both necessary and sufficient for eTreg cell development - contingent on Stat3 output of the IL-2 receptor coordinate with IL-2 signaling during early Treg cell commitment. Induction of iTreg cell programming in absence of IL-2-induced Stat3 signaling resulted in impaired eTreg cell differentiation and a failure to produce IL-10. An IL-2 mutein with reduced affinity for the IL-2Rγ (γ c ) chain was found to have blunted IL-2R Stat3 output, resulting in a deficiency of Il10 transcriptional programming that could not be fully rescued by Stat3 signaling subsequent to an initial window of iTreg cell differentiation. These findings expose a heretofore unappreciated role of IL-2 signaling that acts early to program subsequent production of IL-10 by developing eTreg cells, with broad implications for IL-2-based therapeutic interventions in immune-mediated diseases.

KMT2A associates with PHF5A-PHF14-HMG20A-RAI1 subcomplex in pancreatic cancer stem cells and epigenetically regulates their characteristics

Pancreatic cancer (PC), one of the most aggressive and life-threatening human malignancies, is known for its resistance to cytotoxic therapies. This is increasingly ascribed to the subpopulation of undifferentiated cells, known as pancreatic cancer stem cells (PCSCs), which display greater evolutionary fitness than other tumor cells to evade the cytotoxic effects of chemotherapy. PCSCs are crucial for tumor relapse as they possess 'stem cell-like' features that are characterized by self-renewal and differentiation. However, the molecular mechanisms that maintain the unique characteristics of PCSCs are poorly understood. Here, we identify the histone methyltransferase KMT2A as a physical binding partner of an RNA polymerase-associated PHF5A-PHF14-HMG20A-RAI1 protein subcomplex and an epigenetic regulator of PCSC properties and functions. Targeting the protein subcomplex in PCSCs with a KMT2A-WDR5 inhibitor attenuates their self-renewal capacity, cell viability, and in vivo tumorigenicity.

mTORC1 Regulates the Metabolic Switch of Postnatal Cardiomyocytes During Regeneration

The metabolic switch from glycolysis to fatty acid oxidation in postnatal cardiomyocytes contributes to the loss of the cardiac regenerative potential of the mammalian heart. However, the mechanisms that regulate this metabolic switch remain unclear. The protein kinase complex mechanistic target of rapamycin complex 1 (mTORC1) is a central signaling hub that regulates cellular metabolism and protein synthesis, yet its role during mammalian heart regeneration and postnatal metabolic maturation is undefined. Here, we use immunoblotting, rapamycin treatment, myocardial infarction, and global proteomics to define the role of mTORC1 in postnatal heart development and regeneration. Our results demonstrate that the activity of mTORC1 is dynamically regulated between the regenerating and the non-regenerating hearts. Acute inhibition of mTORC1 by rapamycin or everolimus reduces cardiomyocyte proliferation and inhibits neonatal heart regeneration following injury. Our quantitative proteomic analysis demonstrates that transient inhibition of mTORC1 during neonatal heart injury did not reduce protein synthesis, but rather shifts the cardiac proteome of the neonatal injured heart from glycolysis towards fatty acid oxidation. This indicates that mTORC1 inhibition following injury accelerates the postnatal metabolic switch, which promotes metabolic maturation and impedes cardiomyocyte proliferation and heart regeneration. Taken together, our results define an important role for mTORC1 in regulating postnatal cardiac metabolism and may represent a novel target to modulate cardiac metabolism and promote heart regeneration.

Consistency across multi-omics layers in a drug-perturbed gut microbial community

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

Feeding during hibernation shifts gene expression toward active season levels in brown bears (Ursus arctos)

Hibernation in bears involves a suite of metabolical and physiological changes, including the onset of insulin resistance, that are driven in part by sweeping changes in gene expression in multiple tissues. Feeding bears glucose during hibernation partially restores active season physiological phenotypes, including partial resensitization to insulin, but the molecular mechanisms underlying this transition remain poorly understood. Here, we analyze tissue-level gene expression in adipose, liver, and muscle to identify genes that respond to midhibernation glucose feeding and thus potentially drive postfeeding metabolical and physiological shifts. We show that midhibernation feeding stimulates differential expression in all analyzed tissues of hibernating bears and that a subset of these genes responds specifically by shifting expression toward levels typical of the active season. Inferences of upstream regulatory molecules potentially driving these postfeeding responses implicate peroxisome proliferator-activated receptor gamma (PPARG) and other known regulators of insulin sensitivity, providing new insight into high-level regulatory mechanisms involved in shifting metabolic phenotypes between hibernation and active states.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

Sustained CD28 costimulation is required for self-renewal and differentiation of TCF-1+ PD-1+ CD8 T cells

During persistent antigen stimulation, such as in chronic infections and cancer, CD8 T cells differentiate into a hypofunctional programmed death protein 1-positive (PD-1+) exhausted state. Exhausted CD8 T cell responses are maintained by precursors (Tpex) that express the transcription factor T cell factor 1 (TCF-1) and high levels of the costimulatory molecule CD28. Here, we demonstrate that sustained CD28 costimulation is required for maintenance of antiviral T cells during chronic infection. Low-level CD28 engagement preserved mitochondrial fitness and self-renewal of Tpex, whereas stronger CD28 signaling enhanced glycolysis and promoted Tpex differentiation into TCF-1neg exhausted CD8 T cells (Tex). Furthermore, enhanced differentiation by CD28 engagement did not reduce the Tpex pool. Together, these findings demonstrate that continuous CD28 engagement is needed to sustain PD-1+ CD8 T cells and suggest that increasing CD28 signaling promotes Tpex differentiation into more functional effector-like Tex, possibly without compromising long-term responses.

Adjusting for gene-specific covariates to improve RNA-seq analysis

SUMMARY

This article suggests a novel positive false discovery rate (pFDR) controlling method for testing gene-specific hypotheses using a gene-specific covariate variable, such as gene length. We suppose the null probability depends on the covariate variable. In this context, we propose a rejection rule that accounts for heterogeneity among tests by using two distinct types of null probabilities. We establish a pFDR estimator for a given rejection rule by following Storey's q-value framework. A condition on a type 1 error posterior probability is provided that equivalently characterizes our rejection rule. We also present a suitable procedure for selecting a tuning parameter through cross-validation that maximizes the expected number of hypotheses declared significant. A simulation study demonstrates that our method is comparable to or better than existing methods across realistic scenarios. In data analysis, we find support for our method's premise that the null probability varies with a gene-specific covariate variable.

AVAILABILITY AND IMPLEMENTATION

The source code repository is publicly available at https://github.com/hsjeon1217/conditional_method.

Control of protein synthesis through mRNA pseudouridylation by dyskerin

Posttranscriptional modifications of mRNA have emerged as regulators of gene expression. Although pseudouridylation is the most abundant, its biological role remains poorly understood. Here, we demonstrate that the pseudouridine synthase dyskerin associates with RNA polymerase II, binds to thousands of mRNAs, and is responsible for their pseudouridylation, an action that occurs in chromatin and does not appear to require a guide RNA with full complementarity. In cells lacking dyskerin, mRNA pseudouridylation is reduced, while at the same time, de novo protein synthesis is enhanced, indicating that this modification interferes with translation. Accordingly, mRNAs with fewer pseudouridines due to knockdown of dyskerin are translated more efficiently. Moreover, mRNA pseudouridylation is severely reduced in patients with dyskeratosis congenita caused by inherited mutations in the gene encoding dyskerin (i.e., DKC1). Our findings demonstrate that pseudouridylation by dyskerin modulates mRNA translatability, with important implications for both normal development and disease.

VBASS enables integration of single cell gene expression data in Bayesian association analysis of rare variants

Rare or de novo variants have substantial contribution to human diseases, but the statistical power to identify risk genes by rare variants is generally low due to rarity of genotype data. Previous studies have shown that risk genes usually have high expression in relevant cell types, although for many conditions the identity of these cell types are largely unknown. Recent efforts in single cell atlas in human and model organisms produced large amount of gene expression data. Here we present VBASS, a Bayesian method that integrates single-cell expression and de novo variant (DNV) data to improve power of disease risk gene discovery. VBASS models disease risk prior as a function of expression profiles, approximated by deep neural networks. It learns the weights of neural networks and parameters of Gamma-Poisson likelihood models of DNV counts jointly from expression and genetics data. On simulated data, VBASS shows proper error rate control and better power than state-of-the-art methods. We applied VBASS to published datasets and identified more candidate risk genes with supports from literature or data from independent cohorts. VBASS can be generalized to integrate other types of functional genomics data in statistical genetics analysis.

A beginner's guide into curated analyses of open access datasets for biomarker discovery in neurodegeneration

The discovery of surrogate biomarkers reflecting neuronal dysfunction in neurodegenerative diseases (NDDs) remains an active area of research. To boost these efforts, we demonstrate the utility of publicly available datasets for probing the pathogenic relevance of candidate markers in NDDs. As a starting point, we introduce the readers to several open access resources, which contain gene expression profiles and proteomics datasets from patient studies in common NDDs, including proteomics analyses of cerebrospinal fluid (CSF). Then, we illustrate the method for curated gene expression analyses across select brain regions from four cohorts of Parkinson disease patients (and from one study in common NDDs), probing glutathione biogenesis, calcium signaling and autophagy. These data are complemented by findings of select markers in CSF-based studies in NDDs. Additionally, we enclose several annotated microarray studies, and summarize reports on CSF proteomics across the NDDs, which the readers can utilize for translational purposes. We anticipate that this "beginner's guide" will benefit the research community in NDDs, and would serve as a useful educational tool.

Bromodomain-containing Protein 4 regulates innate inflammation via modulation of alternative splicing

Introduction

Bromodomain-containing Protein 4 (BRD4) is a transcriptional regulator which coordinates gene expression programs controlling cancer biology, inflammation, and fibrosis. In the context of airway viral infection, BRD4-specific inhibitors (BRD4i) block the release of pro-inflammatory cytokines and prevent downstream epithelial plasticity. Although the chromatin modifying functions of BRD4 in inducible gene expression have been extensively investigated, its roles in post-transcriptional regulation are not well understood. Given BRD4's interaction with the transcriptional elongation complex and spliceosome, we hypothesize that BRD4 is a functional regulator of mRNA processing.

Methods

To address this question, we combine data-independent analysis - parallel accumulation-serial fragmentation (diaPASEF) with RNA-sequencing to achieve deep and integrated coverage of the proteomic and transcriptomic landscapes of human small airway epithelial cells exposed to viral challenge and treated with BRD4i.

Results

We discover that BRD4 regulates alternative splicing of key genes, including Interferon-related Developmental Regulator 1 (IFRD1) and X-Box Binding Protein 1 (XBP1), related to the innate immune response and the unfolded protein response (UPR). We identify requirement of BRD4 for expression of serine-arginine splicing factors, splicosome components and the Inositol-Requiring Enzyme 1 IREα affecting immediate early innate response and the UPR.

Discussion

These findings extend the transcriptional elongation-facilitating actions of BRD4 in control of post-transcriptional RNA processing via modulating splicing factor expression in virus-induced innate signaling.

The transcriptome response of astronaut leukocytes to long missions aboard the International Space Station reveals immune modulation

Introduction

Spaceflight leads to the deconditioning of multiple body systems including the immune system. We sought to characterize the molecular response involved by capturing changes in leukocyte transcriptomes from astronauts transitioning to and from long-duration spaceflight.

Methods

Fourteen male and female astronauts with ~6-month- long missions aboard the International Space Station (ISS) had 10 blood samples collected throughout the three phases of the study: one pre-flight (PF), four in-flight (IF) while onboard the ISS, and five upon return to Earth (R). We measured gene expression through RNA sequencing of leukocytes and applied generalized linear modeling to assess differential expression across all 10 time points followed by the analysis of selected time points and functional enrichment of changing genes to identify shifts in biological processes.

Results

Our temporal analysis identified 276 differentially expressed transcripts grouped into two clusters (C) showing opposite profiles of expression with transitions to and from spaceflight: (C1) decrease-then-increase and (C2) increase-then-decrease. Both clusters converged toward average expression between ~2 and ~6 months in space. Further analysis of spaceflight transitions identified the decrease-then-increase pattern with most changes: 112 downregulated genes between PF and early spaceflight and 135 upregulated genes between late IF and R. Interestingly, 100 genes were both downregulated when reaching space and upregulated when landing on Earth. Functional enrichment at the transition to space related to immune suppression increased cell housekeeping functions and reduced cell proliferation. In contrast, egress to Earth is related to immune reactivation.

Conclusion

The leukocytes' transcriptome changes describe rapid adaptations in response to entering space followed by opposite changes upon returning to Earth. These results shed light on immune modulation in space and highlight the major adaptive changes in cellular activity engaged to adapt to extreme environments.

Weighted multiple testing procedures in genome-wide association studies

Multiple testing procedures controlling the false discovery rate (FDR) are increasingly used in the context of genome wide association studies (GWAS), and weighted multiple testing procedures that incorporate covariate information are efficient to improve the power to detect associations. In this work, we evaluate some recent weighted multiple testing procedures in the specific context of GWAS through a simulation study. We also present a new efficient procedure called wBHa that prioritizes the detection of genetic variants with low minor allele frequencies while maximizing the overall detection power. The results indicate good performance of our procedure compared to other weighted multiple testing procedures. In particular, in all simulated settings, wBHa tends to outperform other procedures in detecting rare variants while maintaining good overall power. The use of the different procedures is illustrated with a real dataset.

Coding and noncoding transcriptomes of NODULIN HOMEOBOX (NDX)-deficient Arabidopsis inflorescence

Arabidopsis NODULIN HOMEOBOX (NDX) is a plant-specific transcriptional regulator whose role in small RNA biogenesis and heterochromatin homeostasis has recently been described. Here we extend our previous transcriptomic analysis to the flowering stage of development. We performed mRNA-seq and small RNA-seq measurements on inflorescence samples of wild-type and ndx1-4 mutant (WiscDsLox344A04) Arabidopsis plants. We identified specific groups of differentially expressed genes and noncoding heterochromatic siRNA (hetsiRNA) loci/regions whose transcriptional activity was significantly changed in the absence of NDX. In addition, data obtained from inflorescence were compared with seedling transcriptomics data, which revealed development-specific changes in gene expression profiles. Overall, we provide a comprehensive data source on the coding and noncoding transcriptomes of NDX-deficient Arabidopsis flowers to serve as a basis for further research on NDX function.

LIS1 RNA-binding orchestrates the mechanosensitive properties of embryonic stem cells in AGO2-dependent and independent ways

Lissencephaly-1 (LIS1) is associated with neurodevelopmental diseases and is known to regulate the molecular motor cytoplasmic dynein activity. Here we show that LIS1 is essential for the viability of mouse embryonic stem cells (mESCs), and it governs the physical properties of these cells. LIS1 dosage substantially affects gene expression, and we uncovered an unexpected interaction of LIS1 with RNA and RNA-binding proteins, most prominently the Argonaute complex. We demonstrate that LIS1 overexpression partially rescued the extracellular matrix (ECM) expression and mechanosensitive genes conferring stiffness to Argonaute null mESCs. Collectively, our data transforms the current perspective on the roles of LIS1 in post-transcriptional regulation underlying development and mechanosensitive processes.

The Rab5 effector FERRY links early endosomes with mRNA localization

Localized translation is vital to polarized cells and requires precise and robust distribution of different mRNAs and ribosomes across the cell. However, the underlying molecular mechanisms are poorly understood and important players are lacking. Here, we discovered a Rab5 effector, the five-subunit endosomal Rab5 and RNA/ribosome intermediary (FERRY) complex, that recruits mRNAs and ribosomes to early endosomes through direct mRNA-interaction. FERRY displays preferential binding to certain groups of transcripts, including mRNAs encoding mitochondrial proteins. Deletion of FERRY subunits reduces the endosomal localization of transcripts in cells and has a significant impact on mRNA levels. Clinical studies show that genetic disruption of FERRY causes severe brain damage. We found that, in neurons, FERRY co-localizes with mRNA on early endosomes, and mRNA loaded FERRY-positive endosomes are in close proximity of mitochondria. FERRY thus transforms endosomes into mRNA carriers and plays a key role in regulating mRNA distribution and transport.

PIK3CA mutation in endometriotic epithelial cells promotes viperin-dependent inflammatory response to insulin

Endometrial epithelia are known to harbor cancer driver mutations in the absence of any pathologies, including mutations in PIK3CA. Insulin plays an important role in regulating uterine metabolism during pregnancy, and hyperinsulinemia is associated with conditions impacting fertility. Hyperinsulinemia also promotes cancer, but the direct action of insulin on mutated endometrial epithelial cells is unknown. Here, we treated 12Z endometriotic epithelial cells carrying the PIK3CAH1047R oncogene with insulin and examined transcriptomes by RNA-seq. While cells naively responded to insulin, the magnitude of differential gene expression (DGE) was nine times greater in PIK3CAH1047R cells, representing a synergistic effect between insulin signaling and PIK3CAH1047R expression. Interferon signaling and the unfolded protein response (UPR) were enriched pathways among affected genes. Insulin treatment in wild-type cells activated normal endoplasmic reticulum stress (ERS) response programs, while PIK3CAH1047R cells activated programs necessary to avoid ERS-induced apoptosis. PIK3CAH1047R expression alone resulted in overexpression (OE) of Viperin (RSAD2), which is involved in viral response and upregulated in the endometrium during early pregnancy. The transcriptional changes induced by insulin in PIK3CAH1047R cells were rescued by knockdown of Viperin, while Viperin OE alone was insufficient to induce a DGE response to insulin, suggesting that Viperin is necessary but not sufficient for the synergistic effect of PIK3CAH1047R and insulin treatment. We identified interferon signaling, viral response, and protein targeting pathways that are induced by insulin but dependent on Viperin in PIK3CAH1047R mutant cells. These results suggest that response to insulin signaling is altered in mutated endometriotic epithelial cells.

Escherichia coli Leucine-Responsive Regulatory Protein Bridges DNA In Vivo and Tunably Dissociates in the Presence of Exogenous Leucine

Feast-famine response proteins are a widely conserved class of global regulators in prokaryotes, the most highly studied of which is the Escherichia coli leucine-responsive regulatory protein (Lrp). Lrp senses the environmental nutrition status and subsequently regulates up to one-third of the genes in E. coli, either directly or indirectly. Lrp exists predominantly as octamers and hexadecamers (16mers), where leucine is believed to shift the equilibrium toward the octameric state. In this study, we analyzed the effects of three oligomerization state mutants of Lrp in terms of their ability to bind to DNA and regulate gene expression in response to exogenous leucine. We find that oligomerization beyond dimers is required for Lrp's regulatory activity and that, contrary to previous speculation, exogenous leucine modulates Lrp activity at its target promoters exclusively by inhibiting Lrp binding to DNA. We also show evidence that Lrp binding bridges DNA over length scales of multiple kilobases, revealing a new range of mechanisms for Lrp-mediated transcriptional regulation. IMPORTANCE Leucine-responsive regulatory protein (Lrp) is one of the most impactful regulators in E. coli and other bacteria. Lrp senses nutrient conditions and responds by controlling strategies for virulence, cellular motility, and nutrient acquisition. Despite its importance and being evolutionarily highly conserved across bacteria and archaea, several mysteries remain regarding Lrp, including how it actually responds to leucine to change its regulation of targets. Previous studies have led to the hypothesis that Lrp switches between two states, an octamer (8 Lrp molecules together) and a hexadecamer (16 Lrp molecules together), upon exposure to leucine; these are referred to as different oligomerization states. Here, we show that contrary to previous expectations, it is Lrp's propensity to bind DNA, rather than its oligomerization state, that is directly affected by leucine in the cell's environment. Our new understanding of Lrp activity will aid in identifying and disrupting pathways used by bacteria to cause disease.

Generation of left ventricle-like cardiomyocytes with improved structural, functional, and metabolic maturity from human pluripotent stem cells

Decreased left ventricle (LV) function caused by genetic mutations or injury often leads to debilitating and fatal cardiovascular disease. LV cardiomyocytes are, therefore, a potentially valuable therapeutical target. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) are neither homogeneous nor functionally mature, which reduces their utility. Here, we exploit cardiac development knowledge to instruct differentiation of hPSCs specifically toward LV cardiomyocytes. Correct mesoderm patterning and retinoic acid pathway blocking are essential to generate near-homogenous LV-specific hPSC-CMs (hPSC-LV-CMs). These cells transit via first heart field progenitors and display typical ventricular action potentials. Importantly, hPSC-LV-CMs exhibit increased metabolism, reduced proliferation, and improved cytoarchitecture and functional maturity compared with age-matched cardiomyocytes generated using the standard WNT-ON/WNT-OFF protocol. Similarly, engineered heart tissues made from hPSC-LV-CMs are better organized, produce higher force, and beat more slowly but can be paced to physiological levels. Together, we show that functionally matured hPSC-LV-CMs can be obtained rapidly without exposure to current maturation regimes.

Immune cell dynamics deconvoluted by single-cell RNA sequencing in normothermic machine perfusion of the liver

Normothermic machine perfusion (NMP) has emerged as an innovative organ preservation technique. Developing an understanding for the donor organ immune cell composition and its dynamic changes during NMP is essential. We aimed for a comprehensive characterization of immune cell (sub)populations, cell trafficking and cytokine release during liver NMP. Single-cell transcriptome profiling of human donor livers prior to, during NMP and after transplantation shows an abundance of CXC chemokine receptor 1⁺/2⁺ (CXCR1⁺/CXCR2⁺) neutrophils, which significantly decreased during NMP. This is paralleled by a large efflux of passenger leukocytes with neutrophil predominance in the perfusate. During NMP, neutrophils shift from a pro-inflammatory state towards an aged/chronically activated/exhausted phenotype, while anti-inflammatory/tolerogenic monocytes/macrophages are increased. We herein describe the dynamics of the immune cell repertoire, phenotypic immune cell shifts and a dominance of neutrophils during liver NMP, which potentially contribute to the inflammatory response. Our findings may serve as resource to initiate future immune-interventional studies.

pycoMeth: a toolbox for differential methylation testing from Nanopore methylation calls

We present pycoMeth, a toolbox to store, manage and analyze DNA methylation calls from long-read sequencing data obtained using the Oxford Nanopore Technologies sequencing platform. Building on a novel, rapid-access, read-level and reference-anchored methylation storage format MetH5, we propose efficient algorithms for haplotype aware, multi-sample consensus segmentation and differential methylation testing. We show that MetH5 is more efficient than existing solutions for storing Oxford Nanopore Technologies methylation calls, and carry out benchmarking for pycoMeth segmentation and differential methylation testing, demonstrating increased performance and sensitivity compared to existing solutions designed for short-read methylation data.

Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures

Cancer genomes harbor a broad spectrum of structural variants (SVs) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assembled complex rearrangements, including a 1.55-Mbp chromothripsis event, and we uncover a complex SV pattern termed templated insertion (TI) thread, characterized by short (mostly <1 kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50 kbp in size. TI threads occur in 3% of cancers, with a prevalence up to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read-based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

Outcomes Stratification of Head and Neck Cancer Using Pre- and Post-treatment DNA Methylation From Peripheral Blood

PURPOSE

Established prognostic factors for head and neck squamous cell carcinoma (HNSCC) mostly consist of clinical and tumor features assessed before treatment. We report a novel application of DNA methylation in peripheral blood before and after radiation therapy to further improve outcomes stratification.

METHODS AND MATERIALS

Peripheral blood samples from patients with nonmetastatic HNSCC were obtained for methylation analysis 1 week before and 1 month after radiation therapy. Patients were randomized 1:1 to a Discovery Cohort or a Validation Cohort. In the Discovery Cohort, associations between genome-wide methylation change (posttreatment minus pretreatment) and recurrence-free survival (RFS) as well as overall survival (OS) were evaluated using Cox regression. A methylation risk score (MRS) was then constructed from methylation levels at the top associated sites, filtered for residing within the regulatory regions of genes expressed in cells of hematopoietic lineage. The prognostic value of MRS was separately assessed in the Discovery and Validation Cohorts.

RESULTS

Between December 2013 and September 2018, 115 patients participated in this study. Human papilloma virus negative status, oral cavity cancer, gastrostomy tube insertion, and higher neutrophil count before radiation therapy were associated with shorter RFS and OS (P < .05). Genes downstream of the methylation sites comprising MRS are HIF1A, SF1, LGALS9, and FUT5, involved in hypoxia response, blood cell maturation, and immune modulation. High MRS (in the top third) was significantly associated with worse RFS (hazard ratio [HR], 7.1; 95% confidence interval [CI], 1.4-35.5; P = .016) and OS (HR, 15.9; 95% CI, 1.6-153.6; P = .017) in the Discovery Cohort, independent of the aforementioned risk factors. These findings were replicated in the Validation Cohort, for which high MRS also independently predicted worse RFS (HR, 10.2; 95%, CI 2.4-43.4; P = .002) and OS (HR, 3.7; 95% CI, 1.3-10.4; P = .015).

CONCLUSIONS

We successfully trained and validated a signature of DNA methylation in peripheral blood before and after radiation therapy that stratified outcomes among patients with HNSCC, implicating the potential for genomics-tailored surveillance and consolidation treatment.

Pitfalls in the statistical analysis of microbiome amplicon sequencing data

Microbiome data are characterized by several aspects that make them challenging to analyse statistically: they are compositional, high dimensional and rich in zeros. A large array of statistical methods exist to analyse these data. Some are borrowed from other fields, such as ecology or RNA-sequencing, while others are custom-made for microbiome data. The large range of available methods, and which is continuously expanding, means that researchers have to invest considerable effort in choosing what method(s) to apply. In this paper we list 14 statistical methods or approaches that we think should be generally avoided. In several cases this is because we believe the assumptions behind the method are unlikely to be met for microbiome data. In other cases we see methods that are used in ways they are not intended to be used. We believe researchers would be helped by more critical evaluations of existing methods, as not all methods in use are suitable or have been sufficiently reviewed. We hope this paper contributes to a critical discussion on what methods are appropriate to use in the analysis of microbiome data.

TBX3 is dynamically expressed in pancreatic organogenesis and fine-tunes regeneration

BACKGROUND

The reactivation of genetic programs from early development is a common mechanism for injury-induced organ regeneration. T-box 3 (TBX3) is a member of the T-box family of transcription factors previously shown to regulate pluripotency and subsequent lineage commitment in a number of tissues, including limb and lung. TBX3 is also involved in lung and heart organogenesis. Here, we provide a comprehensive and thorough characterization of TBX3 and its role during pancreatic organogenesis and regeneration.

RESULTS

We interrogated the level and cell specificity of TBX3 in the developing and adult pancreas at mRNA and protein levels at multiple developmental stages in mouse and human pancreas. We employed conditional mutagenesis to determine its role in murine pancreatic development and in regeneration after the induction of acute pancreatitis. We found that Tbx3 is dynamically expressed in the pancreatic mesenchyme and epithelium. While Tbx3 is expressed in the developing pancreas, its absence is likely compensated by other factors after ablation from either the mesenchymal or epithelial compartments. In an adult model of acute pancreatitis, we found that a lack of Tbx3 resulted in increased proliferation and fibrosis as well as an enhanced inflammatory gene programs, indicating that Tbx3 has a role in tissue homeostasis and regeneration.

CONCLUSIONS

TBX3 demonstrates dynamic expression patterns in the pancreas. Although TBX3 is dispensable for proper pancreatic development, its absence leads to altered organ regeneration after induction of acute pancreatitis.

Aerosolized nicotine from e-cigarettes alters gene expression, increases lung protein permeability, and impairs viral clearance in murine influenza infection

E-cigarette use has rapidly increased as an alternative means of nicotine delivery by heated aerosolization. Recent studies demonstrate nicotine-containing e-cigarette aerosols can have immunosuppressive and pro-inflammatory effects, but it remains unclear how e-cigarettes and the constituents of e-liquids may impact acute lung injury and the development of acute respiratory distress syndrome caused by viral pneumonia. Therefore, in these studies, mice were exposed one hour per day over nine consecutive days to aerosol generated by the clinically-relevant tank-style Aspire Nautilus aerosolizing e-liquid containing a mixture of vegetable glycerin and propylene glycol (VG/PG) with or without nicotine. Exposure to the nicotine-containing aerosol resulted in clinically-relevant levels of plasma cotinine, a nicotine-derived metabolite, and an increase in the pro-inflammatory cytokines IL-17A, CXCL1, and MCP-1 in the distal airspaces. Following the e-cigarette exposure, mice were intranasally inoculated with influenza A virus (H1N1 PR8 strain). Exposure to aerosols generated from VG/PG with and without nicotine caused greater influenza-induced production in the distal airspaces of the pro-inflammatory cytokines IFN-γ, TNFα, IL-1β, IL-6, IL-17A, and MCP-1 at 7 days post inoculation (dpi). Compared to the aerosolized carrier VG/PG, in mice exposed to aerosolized nicotine there was a significantly lower amount of Mucin 5 subtype AC (MUC5AC) in the distal airspaces and significantly higher lung permeability to protein and viral load in lungs at 7 dpi with influenza. Additionally, nicotine caused relative downregulation of genes associated with ciliary function and fluid clearance and an increased expression of pro-inflammatory pathways at 7 dpi. These results show that (1) the e-liquid carrier VG/PG increases the pro-inflammatory immune responses to viral pneumonia and that (2) nicotine in an e-cigarette aerosol alters the transcriptomic response to pathogens, blunts host defense mechanisms, increases lung barrier permeability, and reduces viral clearance during influenza infection. In conclusion, acute exposure to aerosolized nicotine can impair clearance of viral infection and exacerbate lung injury, findings that have implications for the regulation of e-cigarette products.

Integrated proteomics reveals alterations in sarcomere composition and developmental processes during postnatal swine heart development

The neonatal swine heart possesses an endogenous ability to regenerate injured myocardium through the proliferation of pre-existing cardiomyocyte (CM) populations. However, this regenerative capacity is lost shortly after birth. Normal postnatal developmental processes and the regenerative capacity of mammalian hearts are tightly linked, but not much is known about how the swine cardiac proteome changes throughout postnatal development. Herein, we integrated robust and quantitative targeted "top-down" and global "bottom-up" proteomic workflows to comprehensively define the dynamic landscape of the swine cardiac proteome throughout postnatal maturation. Using targeted top-down proteomics, we were able to identify significant alterations in sarcomere composition, providing new insight into the proteoform landscape of sarcomeres that can disassemble, a process necessary for productive CM proliferation. Furthermore, we quantified global changes in protein abundance using bottom-up proteomics, identified over 700 differentially expressed proteins throughout postnatal development, and mapped these proteins to changes in developmental and metabolic processes. We envision these results will help guide future investigations to comprehensively understand endogenous cardiac regeneration toward the development of novel therapeutic strategies for heart failure.

Small noncoding RNA interactome capture reveals pervasive, carbon source-dependent tRNA engagement of yeast glycolytic enzymes

Small noncoding RNAs fulfill key functions in cellular and organismal biology, typically working in concert with RNA-binding proteins (RBPs). While proteome-wide methodologies have enormously expanded the repertoire of known RBPs, these methods do not distinguish RBPs binding to small noncoding RNAs from the rest. To specifically identify this relevant subclass of RBPs, we developed small noncoding RNA interactome capture (snRIC_2C) based on the differential RNA-binding capacity of silica matrices (2C). We define the S. cerevisiae proteome of nearly 300 proteins that specifically binds to RNAs smaller than 200 nt in length (snRBPs), identifying informative distinctions from the total RNA-binding proteome determined in parallel. Strikingly, the snRBPs include most glycolytic enzymes from yeast. With further methodological developments using silica matrices, 12 tRNAs were identified as specific binders of the glycolytic enzyme GAPDH. We show that tRNA engagement of GAPDH is carbon source-dependent and regulated by the RNA polymerase III repressor Maf1, suggesting a regulatory interaction between glycolysis and RNA polymerase III activity. We conclude that snRIC_2C and other 2C-derived methods greatly facilitate the study of RBPs, revealing previously unrecognized interactions.

Gut Microbiota Remodeling and Intestinal Adaptation to Lipid Malabsorption After Enteroendocrine Cell Loss in Adult Mice

BACKGROUND & AIMS

Enteroendocrine cells (EECs) and their hormones are essential regulators of whole-body energy homeostasis. EECs sense luminal nutrients and microbial metabolites and subsequently secrete various hormones acting locally or at a distance. Impaired development of EECs during embryogenesis is life-threatening in newborn mice and humans due to compromised nutrient absorption. However, the physiological importance of the EEC system in adult mice has yet to be directedly studied. Herein, we aimed to determine the long-term consequences of a total loss of EECs in healthy adults on energy metabolism, intestinal transcriptome, and microbiota.

METHODS

We depleted intestinal EECs by tamoxifen treatment of adult Neurog3^fl/fl; Villin-CreER^T2 male mice. We studied intestinal cell differentiation, food efficiency, lipid absorption, microbiota composition, fecal metabolites, and transcriptomic responses in the proximal and distal small intestines of mice lacking EECs. We also determined the high-fat diet-induced transcriptomic changes in sorted Neurog3^eYFP/+ EECs.

RESULTS

Induction of EEC deficiency in adults is not life-threatening unless fed with a high-fat diet. Under a standard chow diet, mice lose 10% of weight due to impaired food efficiency. Blood concentrations of cholesterol, triglycerides, and free fatty acids are reduced, and lipid absorption is impaired and delayed in the distal small intestine. Genes controlling lipogenesis, carbohydrate metabolism, and neoglucogenesis are upregulated. Microbiota composition is rapidly altered after EECs depletion and is characterized by decreased a-diversity. Bacteroides and Lactobacillus were progressively enriched, whereas Lachnospiraceae declined without impacting fecal short-chain fatty acid concentrations.

CONCLUSIONS

EECs are dispensable for survival in adult male mice under a standard chow diet. The absence of EECs impairs intestinal lipid absorption, leading to transcriptomic and metabolic adaptations and remodeling of the gut microbiota.

Transcriptome of GH-producing pituitary neuroendocrine tumours and models are significantly affected by somatostatin analogues

Pituitary neuroendocrine tumours (PitNETs) are neoplasms of the pituitary that overproduce hormones or cause unspecific symptoms due to mass effect. Growth hormone overproducing GH-producing PitNETs cause acromegaly leading to connective tissue, metabolic or oncologic disorders. The medical treatment of acromegaly is somatostatin analogues (SSA) in specific cases combined with dopamine agonists (DA), but almost half of patients display partial or full SSA resistance and potential causes of this are unknown. In this study we investigated transcriptomic landscape of GH-producing PitNETs on several levels and functional models-tumour tissue of patients with and without SSA preoperative treatment, tumour derived pituispheres and GH3 cell line incubated with SSA to study effect of medication on gene expression. MGI sequencing platform was used to sequence total RNA from PitNET tissue, pituispheres, mesenchymal stromal stem-like cells (MSC), and GH3 cell cultures, and data were analysed with Salmon-DeSeq2 pipeline. We observed that the GH-producing PitNETs have distinct changes in growth hormone related pathways related to its functional status alongside inner cell signalling, ion transport, cell adhesion and extracellular matrix characteristic patterns. In pituispheres model, treatment regimens (octreotide and cabergoline) affect specific cell proliferation (MKI67) and core functionality pathways (RYR2, COL8A2, HLA-G, ARFGAP1, TGFBR2). In GH3 cells we observed that medication did not have transcriptomic effects similar to preoperative treatment in PitNET tissue or pituisphere model. This study highlights the importance of correct model system selection for cell transcriptomic profiling and data interpretation that could be achieved in future by incorporating NGS methods and detailed cell omics profiling in PitNET model research.

A systematic computational analysis of Mycobacterium tuberculosis H37Rv and human CD34+ genomic expression reveals crucial molecular entities involved in infection progression

The co-evolution of Mycobacterium tuberculosis H37Rv along with its host systems enables the pathogenic bacterium to emerge as a multi-drug resistant form. This creates challenges for a more efficacious treatment strategy that can mitigate the infection. Working towards the same, our study followed a mathematical and statistical approach proposing that mycobacterial transcription factors regulating virulence and adaptation, host cell cytoplasmic component metabolism, oxidoreductase activity and respiratory ETC would be targets for antibiotics against Mycobacterium tuberculosis. Simultaneously, extending the statistical study on Mycobacterium-infected human cord blood CD34+ cells revealed that the human CD34+ genes, S100A8 and FGR (tyrosine-protein kinase, Src2), might be affected in the infection pathogenesis by Mycobacterium. Further, the deduced Mycobacterium-human gene interaction network proposed that mycobacterial coregulators Rv0452 (MarR family regulator) and Rv3862c (WhiB6) triggered genes controlling bacterial metabolism, which influences human immunological pathways involving TLR2 and CXCL8/MAPK8.Communicated by Ramaswamy H. Sarma.

Whole genome resequencing identifies local adaptation associated with environmental variation for redband trout

Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low-coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L, ROCK1, SIX6). Genotype-environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high-elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high-elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

Differences by origin in methylome suggest eco-phenotypes in the kelp Saccharina latissima

Most kelp species are of high ecological and economic importance worldwide, but are highly susceptible to rising ocean temperatures due to their sessile lifestyle. Due to interference with reproduction, development and growth, natural kelp forests have vanished in multiple regions after extreme summer heat waves. Furthermore, increasing temperatures are likely to decrease biomass production and, thus, reduce production security of farmed kelp. Epigenetic variation, and cytosine methylation as a heritable epigenetic trait, is a rapid means of acclimation and adaptation to environmental conditions, including temperature. While the first methylome of brown macroalgae has been recently described in the kelp Saccharina japonica, its functional relevance and contribution to environmental acclimation is currently unknown. The main objective of our study was to identify the importance of the methylome in the congener kelp species Saccharina latissima for temperature acclimation. Our study is the first to compare DNA methylation in kelp between wild populations of different latitudinal origin, and the first to investigate the effect of cultivation and rearing temperature on genome-wide cytosine methylation. Origin appears to determine many traits in kelp, but it is unknown to what extent the effects of thermal acclimation may be overruled by lab-related acclimation. Our results suggest that seaweed hatchery conditions have strong effects on the methylome and, thus, putatively on the epigenetically controlled characteristics of young kelp sporophytes. However, culture origin could best explain epigenetic differences in our samples suggesting that epigenetic mechanisms contribute to local adaptation of eco-phenotypes. Our study is a first step to understand whether DNA methylation marks (via their effect on gene regulation) may be used as biological regulators to enhance production security and kelp restoration success under rising temperatures, and highlights the importance to match hatchery conditions to origin.

Chromatin complex dependencies reveal targeting opportunities in leukemia

Chromatin regulators are frequently mutated in human cancer and are attractive drug targets. They include diverse proteins that share functional domains and assemble into related multi-subunit complexes. To investigate functional relationships among these regulators, here we apply combinatorial CRISPR knockouts (KOs) to test over 35,000 gene-gene pairings in leukemia cells, using a library of over 300,000 constructs. Top pairs that demonstrate either compensatory non-lethal interactions or synergistic lethality enrich for paralogs and targets that occupy the same protein complex. The screen highlights protein complex dependencies not apparent in single KO screens, for example MCM histone exchange, the nucleosome remodeling and deacetylase (NuRD) complex, and HBO1 (KAT7) complex. We explore two approaches to NuRD complex inactivation. Paralog and non-paralog combinations of the KAT7 complex emerge as synergistic lethal and specifically nominate the ING5 PHD domain as a potential therapeutic target when paired with other KAT7 complex member losses. These findings highlight the power of combinatorial screening to provide mechanistic insight and identify therapeutic targets within redundant networks.

Macrophage differentiation is marked by increased abundance of the mRNA 3' end processing machinery, altered poly(A) site usage, and sensitivity to the level of CstF64

Regulation of mRNA polyadenylation is important for response to external signals and differentiation in several cell types, and results in mRNA isoforms that vary in the amount of coding sequence or 3' UTR regulatory elements. However, its role in differentiation of monocytes to macrophages has not been investigated. Macrophages are key effectors of the innate immune system that help control infection and promote tissue-repair. However, overactivity of macrophages contributes to pathogenesis of many diseases. In this study, we show that macrophage differentiation is characterized by shortening and lengthening of mRNAs in relevant cellular pathways. The cleavage/polyadenylation (C/P) proteins increase during differentiation, suggesting a possible mechanism for the observed changes in poly(A) site usage. This was surprising since higher C/P protein levels correlate with higher proliferation rates in other systems, but monocytes stop dividing after induction of differentiation. Depletion of CstF64, a C/P protein and known regulator of polyadenylation efficiency, delayed macrophage marker expression, cell cycle exit, attachment, and acquisition of structural complexity, and impeded shortening of mRNAs with functions relevant to macrophage biology. Conversely, CstF64 overexpression increased use of promoter-proximal poly(A) sites and caused the appearance of differentiated phenotypes in the absence of induction. Our findings indicate that regulation of polyadenylation plays an important role in macrophage differentiation.

Bromodomain-containing Protein 4 Regulates Innate Inflammation in Airway Epithelial Cells via Modulation of Alternative Splicing

Bromodomain-containing Protein 4 (BRD4) is a transcriptional regulator which coordinates gene expression programs controlling cancer biology, inflammation, and fibrosis. In airway viral infection, non-toxic BRD4-specific inhibitors (BRD4i) block the release of pro-inflammatory cytokines and prevent downstream remodeling. Although the chromatin modifying functions of BRD4 in inducible gene expression have been extensively investigated, its roles in post-transcriptional regulation are not as well understood. Based on its interaction with the transcriptional elongation complex and spliceosome, we hypothesize that BRD4 is a functional regulator of mRNA processing. To address this question, we combine data-independent analysis - parallel accumulation-serial fragmentation (diaPASEF) with RNA-sequencing to achieve deep and integrated coverage of the proteomic and transcriptomic landscapes of human small airway epithelial cells exposed to viral challenge and treated with BRD4i. The transcript-level data was further interrogated for alternative splicing analysis, and the resulting data sets were correlated to identify pathways subject to post-transcriptional regulation. We discover that BRD4 regulates alternative splicing of key genes, including Interferon-related Developmental Regulator 1 ( IFRD1 ) and X-Box Binding Protein 1 ( XBP1 ), related to the innate immune response and the unfolded protein response, respectively. These findings extend the transcriptional elongation-facilitating actions of BRD4 in control of post-transcriptional RNA processing in innate signaling.

The mitochondrial multi-omic response to exercise training across tissues

Mitochondria are adaptable organelles with diverse cellular functions critical to whole-body metabolic homeostasis. While chronic endurance exercise training is known to alter mitochondrial activity, these adaptations have not yet been systematically characterized. Here, the Molecular Transducers of Physical Activity Consortium (MoTrPAC) mapped the longitudinal, multi-omic changes in mitochondrial analytes across 19 tissues in male and female rats endurance trained for 1, 2, 4 or 8 weeks. Training elicited substantial changes in the adrenal gland, brown adipose, colon, heart and skeletal muscle, while we detected mild responses in the brain, lung, small intestine and testes. The colon response was characterized by non-linear dynamics that resulted in upregulation of mitochondrial function that was more prominent in females. Brown adipose and adrenal tissues were characterized by substantial downregulation of mitochondrial pathways. Training induced a previously unrecognized robust upregulation of mitochondrial protein abundance and acetylation in the liver, and a concomitant shift in lipid metabolism. The striated muscles demonstrated a highly coordinated response to increase oxidative capacity, with the majority of changes occurring in protein abundance and post-translational modifications. We identified exercise upregulated networks that are downregulated in human type 2 diabetes and liver cirrhosis. In both cases HSD17B10, a central dehydrogenase in multiple metabolic pathways and mitochondrial tRNA maturation, was the main hub. In summary, we provide a multi-omic, cross-tissue atlas of the mitochondrial response to training and identify candidates for prevention of disease-associated mitochondrial dysfunction.

Modulation of RNA splicing enhances response to BCL2 inhibition in leukemia

Therapy resistance is a major challenge in the treatment of cancer. Here, we performed CRISPR-Cas9 screens across a broad range of therapies used in acute myeloid leukemia to identify genomic determinants of drug response. Our screens uncover a selective dependency on RNA splicing factors whose loss preferentially enhances response to the BCL2 inhibitor venetoclax. Loss of the splicing factor RBM10 augments response to venetoclax in leukemia yet is completely dispensable for normal hematopoiesis. Combined RBM10 and BCL2 inhibition leads to mis-splicing and inactivation of the inhibitor of apoptosis XIAP and downregulation of BCL2A1, an anti-apoptotic protein implicated in venetoclax resistance. Inhibition of splicing kinase families CLKs (CDC-like kinases) and DYRKs (dual-specificity tyrosine-regulated kinases) leads to aberrant splicing of key splicing and apoptotic factors that synergize with venetoclax, and overcomes resistance to BCL2 inhibition. Our findings underscore the importance of splicing in modulating response to therapies and provide a strategy to improve venetoclax-based treatments.

Multi-OMICs landscape of SARS-CoV-2-induced host responses in human lung epithelial cells

COVID-19 pandemic continues to remain a global health concern owing to the emergence of newer variants. Several multi-Omics studies have produced extensive evidence on host-pathogen interactions and potential therapeutic targets. Nonetheless, an increased understanding of host signaling networks regulated by post-translational modifications and their ensuing effect on the cellular dynamics is critical to expanding the current knowledge on SARS-CoV-2 infections. Through an unbiased transcriptomics, proteomics, acetylomics, phosphoproteomics, and exometabolome analysis of a lung-derived human cell line, we show that SARS-CoV-2 Norway/Trondheim-S15 strain induces time-dependent alterations in the induction of type I IFN response, activation of DNA damage response, dysregulated Hippo signaling, among others. We identified interplay of phosphorylation and acetylation dynamics on host proteins and its effect on the altered release of metabolites, especially organic acids and ketone bodies. Together, our findings serve as a resource of potential targets that can aid in designing novel host-directed therapeutic strategies.

A KO mouse model for the lncRNA Lhx1os produces motor neuron alterations and locomotor impairment

Here, we describe a conserved motor neuron-specific long non-coding RNA, Lhx1os, whose knockout in mice produces motor impairment and postnatal reduction of mature motor neurons (MNs). The ER stress-response pathway result specifically altered with the downregulation of factors involved in the unfolded protein response (UPR). Lhx1os was found to bind the ER-associated PDIA3 disulfide isomerase and to affect the expression of the same set of genes controlled by this protein, indicating that the two factors act in conjunction to modulate the UPR. Altogether, the observed phenotype and function of Lhx1os indicate its important role in the control of MN homeostasis and function.

High-resolution single-cell atlas reveals diversity and plasticity of tissue-resident neutrophils in non-small cell lung cancer

Non-small cell lung cancer (NSCLC) is characterized by molecular heterogeneity with diverse immune cell infiltration patterns, which has been linked to therapy sensitivity and resistance. However, full understanding of how immune cell phenotypes vary across different patient subgroups is lacking. Here, we dissect the NSCLC tumor microenvironment at high resolution by integrating 1,283,972 single cells from 556 samples and 318 patients across 29 datasets, including our dataset capturing cells with low mRNA content. We stratify patients into immune-deserted, B cell, T cell, and myeloid cell subtypes. Using bulk samples with genomic and clinical information, we identify cellular components associated with tumor histology and genotypes. We then focus on the analysis of tissue-resident neutrophils (TRNs) and uncover distinct subpopulations that acquire new functional properties in the tissue microenvironment, providing evidence for the plasticity of TRNs. Finally, we show that a TRN-derived gene signature is associated with anti-programmed cell death ligand 1 (PD-L1) treatment failure.