gprofiler2 -- an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e., cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data (https://viz.datascience.arizona.edu/3DEpiEx/).

Comparative multi-omic analyses of cardiac mitochondrial stress in three mouse models of frataxin deficiency

Cardiomyopathy is often fatal in Friedreich ataxia (FA). However, FA hearts maintain adequate function until advanced disease stages, suggesting initial adaptation to the loss of frataxin (FXN). Conditional cardiac knockout mouse models of FXN show transcriptional and metabolic profiles of the mitochondrial integrated stress response (ISRmt), which could play an adaptive role. However, the ISRmt has not been investigated in models with disease-relevant, partial decrease in FXN. We characterized the heart transcriptomes and metabolomes of three mouse models with varying degrees of FXN depletion: YG8-800, KIKO-700 and FXNG127V. Few metabolites were changed in YG8-800 mice, which did not provide a signature of cardiomyopathy or ISRmt; several metabolites were altered in FXNG127V and KIKO-700 hearts. Transcriptional changes were found in all models, but differentially expressed genes consistent with cardiomyopathy and ISRmt were only identified in FXNG127V hearts. However, these changes were surprisingly mild even at advanced age (18 months), despite a severe decrease in FXN levels to 1% of those of wild type. These findings indicate that the mouse heart has low reliance on FXN, highlighting the difficulty in modeling genetically relevant FA cardiomyopathy.

Molecular classification and identification of an aggressive signature in low-grade B-cell lymphomas

Non-follicular low-grade B-cell lymphomas (LGBCL) are biologically diverse entities that share clinical and histologic features that make definitive pathologic categorization challenging. While most patients with LGBCL have an indolent course, some experience aggressive disease, highlighting additional heterogeneity across these subtypes. To investigate the potential for shared biology across subtypes, we performed RNA sequencing and applied machine learning approaches that identified five clusters of patients that grouped independently of subtype. One cluster was characterized by inferior outcome, upregulation of cell cycle genes, and increased tumor immune cell content. Integration of whole exome sequencing identified novel LGBCL mutations and enrichment of TNFAIP3 and BCL2 alterations in the poor survival cluster. Building on this, we further refined a transcriptomic signature associated with early clinical failure in two independent cohorts. Taken together, this study identifies unique clusters of LGBCL defined by novel gene expression signatures and immune profiles associated with outcome across diagnostic subtypes.

Differential voluntary feed intake and whole transcriptome profiling in the hypothalamus of young sheep offered CP and phosphorus-deficient diets

A reduction in voluntary feed intake is observed in ruminants consuming nutrient-deficient diets, such as those with a low CP or P content, and has been attributed to active metabolic regulation, rather than a physical constraint. The hypothalamus is the key integrator of feed intake regulation in mammals. The objectives of this experiment were to (1) establish a model of metabolic feed intake regulation in ruminants consuming diets of variable CP and P content, and (2) determine key biochemical pathways and influential points of regulation within the hypothalamus. Merino wethers [n = 40; 23.7 ± 1.4 kg liveweight (mean ± SD)] were fed one of five dietary treatments (n = 8/treatment) for 63 days in individual pens. The treatments included targeted combinations of high (H) and low (L) CP (110 and 55 g/kg DM) and high and low P (2.5 and 0.7 g/kg DM) with 9 MJ metabolisable energy (ME) per kg DM which were fed ad libitum (UMEI; unrestricted ME intake) resulting in four experimental diets (HCP-HP-UMEI, LCP-HP-UMEI, HCP-LP-UMEI and LCP-LP-UMEI). An additional nutritional treatment (HCP-HP-RMEI) restricted intake of the HCP-HP diet to an equivalent ME intake of wethers consuming the LCP-LP-UMEI treatment. Wethers offered the LCP-HP-UMEI, HCP-LP-UMEI and LCP-LP-UMEI treatments consumed 42, 32 and 49% less total DM (P ≤ 0.05), respectively than the HCP-HP-UMEI treatment, and this was not attributable to any physical limitation of the rumen. Plasma concentrations of urea nitrogen and inorganic phosphate indicated that these nutrient deficiencies were successfully established. To assess potential mechanisms, RNA-seq was conducted on samples from the arcuate nucleus (ARC), ventromedial hypothalamus and lateral hypothalamus of the wethers, yielding a total of 301, 8 and 148 differentially expressed genes across all pairwise comparisons, respectively. The expression of NPY, AGRP and CARTPT, known for their regulatory role in mammalian feed intake regulation, had a similar transcriptional response in the ARC of wethers consuming nutrient-deficient treatments and those consuming a ME-restricted treatment, despite these wethers expressing behaviours indicative of satiated and hungry states, respectively. In addition, genes involved with glycolysis (TPI1), the citric acid cycle (CS, OGDH, GLUD1, GOT1) and oxidative phosphorylation (COX5A, ATP5MC1, ATP5F1B, ATP5MC3) were downregulated in the ARC of wethers fed a nutrient deficient (LCP-LP-UMEI) relative to the non-deficient (HCP-HP-UMEI) treatment. In summary, a model for voluntary feed intake restriction was established to determine genome-wide molecular changes in the hypothalamus of young ruminants.

Distinct impacts of alpha-synuclein overexpression on the hippocampal epigenome of mice in standard and enriched environments

Elevated alpha-synuclein (SNCA) gene expression is associated with transcriptional deregulation and increased risk of Parkinson's disease, which may be partially ameliorated by environmental enrichment. At the molecular level, there is emerging evidence that excess alpha-synuclein protein (aSyn) impacts the epigenome through direct and/or indirect mechanisms. However, the extents to which the effects of both aSyn and the environment converge at the epigenome and whether epigenetic alterations underpin the preventive effects of environmental factors on transcription remain to be elucidated. Here, we profiled five DNA and histone modifications in the hippocampus of wild-type and transgenic mice overexpressing human SNCA. Mice of each genotype were housed under either standard conditions or in an enriched environment (EE) for 12 months. SNCA overexpression induced hippocampal CpG hydroxymethylation and histone H3K27 acetylation changes that associated with genotype more than environment. Excess aSyn was also associated with genotype- and environment-dependent changes in non-CpG (CpH) DNA methylation and H3K4 methylation. These H3K4 methylation changes included loci where the EE ameliorated the impacts of the transgene as well as loci resistant to the effects of environmental enrichment in transgenic mice. In addition, select H3K4 monomethylation alterations were associated with changes in mRNA expression. Our results suggested an environment-dependent impact of excess aSyn on some functionally relevant parts of the epigenome, and will ultimately enhance our understanding of the molecular etiology of Parkinson's disease and other synucleinopathies.

DISTINCT METABOLIC STATES DIRECT RETINAL PIGMENT EPITHELIUM CELL FATE DECISIONS

During tissue regeneration, proliferation, dedifferentiation, and reprogramming are necessary to restore lost structures. However, it is not fully understood how metabolism intersects with these processes. Chicken embryos can regenerate their retina through retinal pigment epithelium (RPE) reprogramming when treated with fibroblast factor 2 (FGF2). Using transcriptome profiling, we uncovered extensive regulation of gene sets pertaining to proliferation, neurogenesis, and glycolysis throughout RPE-to-neural retina reprogramming. By manipulating cell media composition, we determined that glucose, glutamine, or pyruvate are sufficient to support RPE reprogramming identifying glycolysis as a requisite. Conversely, the induction of oxidative metabolism by activation of pyruvate dehydrogenase induces Epithelial-to-mesenchymal transition (EMT), while simultaneously blocking the activation of neural retina fate. We also identify that EMT is partially driven by an oxidative environment. Our findings provide evidence that metabolism controls RPE cell fate decisions and provide insights into the metabolic state of RPE cells, which are prone to fate changes in regeneration and pathologies, such as proliferative vitreoretinopathy.

Single-cell RNA sequencing and analysis of rodent blood stage Plasmodium

Bulk RNA sequencing of Plasmodium spp., the causative parasite of malaria, fails to discriminate developmental-stage-specific gene regulation. Here, we provide a protocol that uses single-cell RNA sequencing of FACS-sorted Plasmodium-chabaudi-chabaudi-AS-infected red blood cells (iRBCs) to characterize developmental-stage-specific modulation of gene expression during malaria blood stage. We describe steps for infecting mice, monitoring disease progression, preparing iRBCs, and single-cell sequencing iRBCs. We then detail procedures for analyzing scRNA-seq data. For complete details on the use and execution of this protocol, please refer to Ramos et al.1.

Multi-omic approach identifies hypoxic tumor-associated myeloid cells that drive immunobiology of high-risk pediatric ependymoma

Ependymoma (EPN) is a devastating childhood brain tumor. Single-cell analyses have illustrated the cellular heterogeneity of EPN tumors, identifying multiple neoplastic cell states including a mesenchymal-differentiated subpopulation which characterizes the PFA1 subtype. Here, we characterize the EPN immune environment, in the context of both tumor subtypes and tumor cell subpopulations using single-cell sequencing (scRNAseq, n = 27), deconvolution of bulk tumor gene expression (n = 299), spatial proteomics (n = 54), and single-cell cytokine release assays (n = 12). We identify eight distinct myeloid-derived subpopulations from which a group of cells, termed hypoxia myeloid cells, demonstrate features of myeloid-derived suppressor cells, including IL6/STAT3 pathway activation and wound healing ontologies. In PFA tumors, hypoxia myeloid cells colocalize with mesenchymal-differentiated cells in necrotic and perivascular niches and secrete IL-8, which we hypothesize amplifies the EPN immunosuppressive microenvironment. This myeloid cell-driven immunosuppression will need to be targeted for immunotherapy to be effective in this difficult-to-cure childhood brain tumor.

MBMethPred: a computational framework for the accurate classification of childhood medulloblastoma subgroups using data integration and AI-based approaches

Childhood medulloblastoma is a malignant form of brain tumor that is widely classified into four subgroups based on molecular and genetic characteristics. Accurate classification of these subgroups is crucial for appropriate treatment, monitoring plans, and targeted therapies. However, misclassification between groups 3 and 4 is common. To address this issue, an AI-based R package called MBMethPred was developed based on DNA methylation and gene expression profiles of 763 medulloblastoma samples to classify subgroups using machine learning and neural network models. The developed prediction models achieved a classification accuracy of over 96% for subgroup classification by using 399 CpGs as prediction biomarkers. We also assessed the prognostic relevance of prediction biomarkers using survival analysis. Furthermore, we identified subgroup-specific drivers of medulloblastoma using functional enrichment analysis, Shapley values, and gene network analysis. In particular, the genes involved in the nervous system development process have the potential to separate medulloblastoma subgroups with 99% accuracy. Notably, our analysis identified 16 genes that were specifically significant for subgroup classification, including EP300, CXCR4, WNT4, ZIC4, MEIS1, SLC8A1, NFASC, ASCL2, KIF5C, SYNGAP1, SEMA4F, ROR1, DPYSL4, ARTN, RTN4RL1, and TLX2. Our findings contribute to enhanced survival outcomes for patients with medulloblastoma. Continued research and validation efforts are needed to further refine and expand the utility of our approach in other cancer types, advancing personalized medicine in pediatric oncology.

Deciphering the decline of metabolic elasticity in aging and obesity

Organisms must adapt to fluctuating nutrient availability to maintain energy homeostasis. Here, we term the capacity for such adaptation and restoration "metabolic elasticity" and model it through ad libitum-fasting-refeeding cycles. Metabolic elasticity is achieved by coordinate versatility in gene expression, which we call "gene elasticity." We have developed the gene elasticity score as a systematic method to quantify the elasticity of the transcriptome across metabolically active tissues in mice and non-human primates. Genes involved in lipid and carbohydrate metabolism show high gene elasticity, and their elasticity declines with age, particularly with PPARγ dysregulation in adipose tissue. Synchronizing PPARγ activity with nutrient conditions through feeding-timed agonism optimizes their metabolic benefits and safety. We further broaden the conceptual scope of metabolic and gene elasticity to dietary challenges, revealing declines in diet-induced obesity similar to those in aging. Altogether, our findings provide a dynamic perspective on the dysmetabolic consequences of aging and obesity.

Self-formation of concentric zones of telencephalic and ocular tissues and directional retinal ganglion cell axons

The telencephalon and eye in mammals are originated from adjacent fields at the anterior neural plate. Morphogenesis of these fields generates telencephalon, optic-stalk, optic-disc, and neuroretina along a spatial axis. How these telencephalic and ocular tissues are specified coordinately to ensure directional retinal ganglion cell (RGC) axon growth is unclear. Here, we report self-formation of human telencephalon-eye organoids comprising concentric zones of telencephalic, optic-stalk, optic-disc, and neuroretinal tissues along the center-periphery axis. Initially-differentiated RGCs grew axons towards and then along a path defined by adjacent PAX2+ VSX2+ optic-disc cells. Single-cell RNA sequencing of these organoids not only confirmed telencephalic and ocular identities but also identified expression signatures of early optic-disc, optic-stalk, and RGCs. These signatures were similar to those in human fetal retinas. Optic-disc cells in these organoids differentially expressed FGF8 and FGF9; FGFR inhibitions drastically decreased early RGC differentiation and directional axon growth. Through the RGC-specific cell-surface marker CNTN2 identified here, electrophysiologically excitable RGCs were isolated under a native condition. Our findings provide insight into the coordinated specification of early telencephalic and ocular tissues in humans and establish resources for studying RGC-related diseases such as glaucoma.

Epigenome-wide DNA methylation association study of circulating IgE levels identifies novel targets for asthma

BACKGROUND

Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases.

METHODS

We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables.

FINDINGS

We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P < 7.94E-04) that can be explored as potential therapeutic targets.

INTERPRETATION

By integrating eQTM and MR analyses in general and clinical asthma populations, our findings provide a deeper understanding of the multidimensional inter-relations of DNA methylation, gene expression, and IgE levels.

FUNDING

US NIH/NHLBI grants: P01HL132825, K99HL159234. N01-HC-25195 and HHSN268201500001I.

Evaluating cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico

BACKGROUND

Preclinical models like cancer cell lines and patient-derived xenografts (PDXs) are vital for studying disease mechanisms and evaluating treatment options. It is essential that they accurately recapitulate the disease state of interest to generate results that will translate in the clinic. Prior studies have demonstrated that preclinical models do not recapitulate all biological aspects of human tissues, particularly with respect to the tissue of origin gene expression signatures. Therefore, it is critical to assess how well preclinical model gene expression profiles correlate with human cancer tissues to inform preclinical model selection and data analysis decisions.

AIMS

Here we evaluated how well preclinical models recapitulate human cancer and non-diseased tissue gene expression patterns in silico with respect to the full gene expression profile as well as subsetting by the most variable genes, genes significantly correlated with tumor purity, and tissue-specific genes.

METHODS

By using publicly available gene expression profiles across multiple sources, we evaluated cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico.

RESULTS

We found that using the full gene set improves correlations between preclinical model and tissue global gene expression profiles, confirmed that glioblastoma (GBM) PDX global gene expression correlation to GBM tumor global gene expression outperforms GBM cell line to GBM tumor global gene expression correlations, and demonstrated that preclinical models in our study often failed to reproduce tissue-specific expression. While including additional genes for global gene expression comparison between cell lines and tissues decreases the overall correlation, it improves the relative rank between a cell line and its tissue of origin compared to other tissues. Our findings underscore the importance of using the full gene expression set measured when comparing preclinical models and tissues and confirm that tissue-specific patterns are better preserved in GBM PDX models than in GBM cell lines.

CONCLUSION

Future studies can build on these findings to determine the specific pathways and gene sets recapitulated by particular preclinical models to facilitate model selection for a given study design or goal.

Transcriptomic Profiling of Subcutaneous Backfat in Castrated and Intact Alentejano Pigs Finished Outdoors with Commercial and Fiber-Rich Diets

In this work, we studied the backfat transcriptome of surgically castrated (C), intact (I) and intact fed an experimental diet (IE) outdoor-reared male Alentejano (AL) pigs. The experimental diet was a high-fiber diet with locally produced legumes and by-products associated with a boar taint reduction effect. At slaughter (~160 kg), backfat samples were collected for total RNA sequencing. Intact pigs presented leaner carcasses, more total collagen, and more unsaturated intramuscular fat content than C animals. A total of 2726 differentially expressed genes (DEGs, |log2 FC|> 0.58, q < 0.05) were identified between C and I with overexpressed genes related to muscular activity (MYH1, ACTA1) or collagen metabolism (COL1A1, COL1A2) in I pigs. Between C and IE, 1639 DEGs of genes involved in lipidic metabolism (LEP, ME1, FABP4, ELOVL6) were overexpressed in C. Finally, only 28 DEGs were determined between I and IE. Clustering results indicated a drastic influence of the testis in the transcriptome of subcutaneous fat of AL pigs, while the diet had a marginal effect. Diet can reduce stress by increasing satiety in animals, and could have induced an increase of skatole degradation due to the higher expression of the CYP2A19 gene in the IE group.

MyD88-dependent Toll-like receptor 2 signaling modulates macrophage activation on lysate-adsorbed Teflon™ AF surfaces in an in vitro biomaterial host response model

The adsorbed protein layer on an implanted biomaterial surface is known to mediate downstream cell-material interactions that drive the host response. While the adsorption of plasma-derived proteins has been studied extensively, the adsorption of damage-associated molecular patterns (DAMPs) derived from damaged cells and matrix surrounding the implant remains poorly understood. Previously, our group developed a DAMP-adsorption model in which 3T3 fibroblast lysates were used as a complex source of cell-derived DAMPs and we demonstrated that biomaterials with adsorbed lysate potently activated RAW-Blue macrophages via Toll-like receptor 2 (TLR2). In the present study, we characterized the response of mouse bone marrow derived macrophages (BMDM) from wildtype (WT), TLR2-/- and MyD88-/- mice on Teflon™ AF surfaces pre-adsorbed with 10% plasma or lysate-spiked plasma (10% w/w total protein from 3T3 fibroblast lysate) for 24 hours. WT BMDM cultured on adsorbates derived from 10% lysate in plasma had significantly higher gene and protein expression of IL-1β, IL-6, TNF-α, IL-10, RANTES/CCL5 and CXCL1/KC, compared to 10% plasma-adsorbed surfaces. Furthermore, the upregulation of pro-inflammatory cytokine and chemokine expression in the 10% lysate in plasma condition was attenuated in TLR2-/- and MyD88-/- BMDM. Proteomic analysis of the adsorbed protein layers showed that even this relatively small addition of lysate-derived proteins within plasma (10% w/w) caused a significant change to the adsorbed protein profile. The 10% plasma condition had fibrinogen, albumin, apolipoproteins, complement, and fibronectin among the top 25 most abundant proteins. While proteins layers generated from 10% lysate in plasma retained fibrinogen and fibronectin among the top 25 proteins, there was a disproportionate increase in intracellular proteins, including histones, tubulins, actins, and vimentin. Furthermore, we identified 7 DAMPs or DAMP-related proteins enriched in the 10% plasma condition (fibrinogen, apolipoproteins), compared to 39 DAMPs enriched in the 10% lysate in plasma condition, including high mobility group box 1 and histones. Together, these findings indicate that DAMPs and other intracellular proteins readily adsorb to biomaterial surfaces in competition with plasma proteins, and that adsorbed DAMPs induce an inflammatory response in adherent macrophages that is mediated by the MyD88-dependent TLR2 signaling pathway.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

α-parvin controls chondrocyte column formation and regulates long bone development

Endochondral ossification requires proper control of chondrocyte proliferation, differentiation, survival, and organization. Here we show that knockout of α-parvin, an integrin-associated focal adhesion protein, from murine limbs causes defects in endochondral ossification and dwarfism. The mutant long bones were shorter but wider, and the growth plates became disorganized, especially in the proliferative zone. With two-photon time-lapse imaging of bone explant culture, we provide direct evidence showing that α-parvin regulates chondrocyte rotation, a process essential for chondrocytes to form columnar structure. Furthermore, loss of α-parvin increased binucleation, elevated cell death, and caused dilation of the resting zones of mature growth plates. Single-cell RNA-seq analyses revealed alterations of transcriptome in all three zones (i.e., resting, proliferative, and hypertrophic zones) of the growth plates. Our results demonstrate a crucial role of α-parvin in long bone development and shed light on the cellular mechanism through which α-parvin regulates the longitudinal growth of long bones.

LINE-1 regulates cortical development by acting as long non-coding RNAs

Long Interspersed Nuclear Elements-1s (L1s) are transposable elements that constitute most of the genome's transcriptional output yet have still largely unknown functions. Here we show that L1s are required for proper mouse brain corticogenesis operating as regulatory long non-coding RNAs. They contribute to the regulation of the balance between neuronal progenitors and differentiation, the migration of post-mitotic neurons and the proportions of different cell types. In cortical cultured neurons, L1 RNAs are mainly associated to chromatin and interact with the Polycomb Repressive Complex 2 (PRC2) protein subunits enhancer of Zeste homolog 2 (Ezh2) and suppressor of zeste 12 (Suz12). L1 RNA silencing influences PRC2's ability to bind a portion of its targets and the deposition of tri-methylated histone H3 (H3K27me3) marks. Our results position L1 RNAs as crucial signalling hubs for genome-wide chromatin remodelling, enabling the fine-tuning of gene expression during brain development and evolution.

Inhibition of Toll-like receptor 4 and Interleukin-1 receptor prevent SARS-CoV-2 mediated kidney injury

Acute kidney injury (AKI) is a common and severe complication of the coronavirus disease 2019 (COVID-19). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) directly affects the glomerular and tubular epithelial cells to induce AKI; however, its pathophysiology remains unclear. Here, we explored the underlying mechanisms and therapeutic targets of renal involvement in COVID-19. We developed an in vitro human kidney cellular model, including immortalized tubular epithelial and endothelial cell lines, demonstrating that SARS-CoV-2 directly triggers cell death. To identify the molecular targets in the process of SARS-CoV-2-mediated cell injury, we performed transcriptional analysis using RNA sequencing. Tubular epithelial cells were more prone to dying by SARS-CoV-2 than endothelial cells; however, SARS-CoV-2 did not replicate in renal cells, distinct from VeroE6/transmembrane protease serine 2 cells. Transcriptomic analysis revealed increased inflammatory and immune-related gene expression levels in renal cells incubated with SARS-CoV-2. Toll-like receptor (TLR) 3 in renal cells recognized viral RNA and underwent cell death. Furthermore, analysis of upstream regulators identified several key transcriptional regulators. Among them, inhibition of the interleukin-1 receptor (IL-1R) and TLR4 pathways protects tubular epithelial and endothelial cells from injury via regulation of the signal transducer and activator of transcription protein-3/nuclear factor-kB pathway. Our results reveal that SARS-CoV-2 directly injures renal cells via the proinflammatory response without viral replication, and that IL-1R and TLR4 may be used as therapeutic targets for SARS-CoV-2 mediated kidney injury.

Large scale microfluidic CRISPR screening for increased amylase secretion in yeast

Key to our ability to increase recombinant protein production through secretion is a better understanding of the pathways that interact to translate, process and export mature proteins to the surrounding environment, including the supporting cellular machinery that supplies necessary energy and building blocks. By combining droplet microfluidic screening with large-scale CRISPR libraries that perturb the expression of the majority of coding and non-coding genes in S. cerevisiae, we identified 345 genes for which an increase or decrease in gene expression resulted in increased secretion of α-amylase. Our results show that modulating the expression of genes involved in the trafficking of vesicles, endosome to Golgi transport, the phagophore assembly site, the cell cycle and energy supply improve α-amylase secretion. Besides protein-coding genes, we also find multiple long non-coding RNAs enriched in the vicinity of genes associated with endosomal, Golgi and vacuolar processes. We validated our results by overexpressing or deleting selected genes, which resulted in significant improvements in α-amylase secretion. The advantages, in terms of precision and speed, inherent to CRISPR based perturbations, enables iterative testing of new strains for increased protein secretion.

Mitochondrial sulfide promotes life span and health span through distinct mechanisms in developing versus adult treated Caenorhabditis elegans

Living longer without simultaneously extending years spent in good health ("health span") is an increasing societal burden, demanding new therapeutic strategies. Hydrogen sulfide (H2S) can correct disease-related mitochondrial metabolic deficiencies, and supraphysiological H2S concentrations can pro health span. However, the efficacy and mechanisms of mitochondrion-targeted sulfide delivery molecules (mtH2S) administered across the adult life course are unknown. Using a Caenorhabditis elegans aging model, we compared untargeted H2S (NaGYY4137, 100 µM and 100 nM) and mtH2S (AP39, 100 nM) donor effects on life span, neuromuscular health span, and mitochondrial integrity. H2S donors were administered from birth or in young/middle-aged animals (day 0, 2, or 4 postadulthood). RNAi pharmacogenetic interventions and transcriptomics/network analysis explored molecular events governing mtH2S donor-mediated health span. Developmentally administered mtH2S (100 nM) improved life/health span vs. equivalent untargeted H2S doses. mtH2S preserved aging mitochondrial structure, content (citrate synthase activity) and neuromuscular strength. Knockdown of H2S metabolism enzymes and FoxO/daf-16 prevented the positive health span effects of mtH2S, whereas DCAF11/wdr-23 - Nrf2/skn-1 oxidative stress protection pathways were dispensable. Health span, but not life span, increased with all adult-onset mtH2S treatments. Adult mtH2S treatment also rejuvenated aging transcriptomes by minimizing expression declines of mitochondria and cytoskeletal components, and peroxisome metabolism hub components, under mechanistic control by the elt-6/elt-3 transcription factor circuit. H2S health span extension likely acts at the mitochondrial level, the mechanisms of which dissociate from life span across adult vs. developmental treatment timings. The small mtH2S doses required for health span extension, combined with efficacy in adult animals, suggest mtH2S is a potential healthy aging therapeutic.

Repurposing of the multiciliation gene regulatory network in fate specification of Cajal-Retzius neurons

Cajal-Retzius cells (CRs) are key players in cerebral cortex development, and they display a unique transcriptomic identity. Here, we use scRNA-seq to reconstruct the differentiation trajectory of mouse hem-derived CRs, and we unravel the transient expression of a complete gene module previously known to control multiciliogenesis. However, CRs do not undergo centriole amplification or multiciliation. Upon deletion of Gmnc, the master regulator of multiciliogenesis, CRs are initially produced but fail to reach their normal identity resulting in their massive apoptosis. We further dissect the contribution of multiciliation effector genes and identify Trp73 as a key determinant. Finally, we use in utero electroporation to demonstrate that the intrinsic competence of hem progenitors as well as the heterochronic expression of Gmnc prevent centriole amplification in the CR lineage. Our work exemplifies how the co-option of a complete gene module, repurposed to control a distinct process, may contribute to the emergence of novel cell identities.

Integrated Transcriptome Analysis Identified Key Expansin Genes Associated with Wheat Cell Wall, Grain Weight and Yield

This research elucidates the dynamic expression of expansin genes during the wheat grain (Triticum aestivum L.) development process using comprehensive meta-analysis and experimental validation. We leveraged RNA-seq data from multiple public databases, applying stringent criteria for selection, and identified 60,852 differentially expressed genes across developmental stages. From this pool, 28,558 DEGs were found to exhibit significant temporal regulation in at least two different datasets and were enriched for processes integral to grain development such as carbohydrate metabolism and cell wall organization. Notably, 30% of the 241 known expansin genes showed differential expression during grain growth. Hierarchical clustering and expression level analysis revealed temporal regulation and distinct contributions of expansin subfamilies during the early stages of grain development. Further analysis using co-expression networks underscored the significance of expansin genes, revealing their substantial co-expression with genes involved in cell wall modification. Finally, qPCR validation and grain morphological analysis under field conditions indicated a significant negative correlation between the expression of select expansin genes, and grain size and weight. This study illuminates the potential role of expansin genes in wheat grain development and provides new avenues for targeted genetic improvements in wheat.

Stimulation of Immune Checkpoint Molecule B and T-Lymphocyte Attenuator Alleviates Experimental Crescentic Glomerulonephritis

SIGNIFICANCE STATEMENT

Treatment of acute, crescentic glomerulonephritis (GN) consists of unspecific and potentially toxic immunosuppression. T cells are central in the pathogenesis of GN, and various checkpoint molecules control their activation. The immune checkpoint molecule B and T-lymphocyte attenuator (BTLA) has shown potential for restraining inflammation in other T-cell-mediated disease models. To investigate its role in GN in a murine model of crescentic nephritis, the authors induced nephrotoxic nephritis in BTLA-deficient mice and wild-type mice. They found that BTLA has a renoprotective role through suppression of local Th1-driven inflammation and expansion of T regulatory cells and that administration of an agonistic anti-BTLA antibody attenuated experimental GN. These findings suggest that antibody-based modulation of BTLA may represent a treatment strategy in human glomerular disease.

BACKGROUND

Modulating T-lymphocytes represents a promising targeted therapeutic option for glomerulonephritis (GN) because these cells mediate damage in various experimental and human GN types. The immune checkpoint molecule B and T-lymphocyte attenuator (BTLA) has shown its potential to restrain inflammation in other T-cell-mediated disease models. Its role in GN, however, has not been investigated.

METHODS

We induced nephrotoxic nephritis (NTN), a mouse model of crescentic GN, in Btla -deficient ( BtlaKO ) mice and wild-type littermate controls and assessed disease severity using functional and histologic parameters at different time points after disease induction. Immunologic changes were comprehensively evaluated by flow cytometry, RNA sequencing, and in vitro assays for dendritic cell and T-cell function. Transfer experiments into Rag1KO mice confirmed the observed in vitro findings. In addition, we evaluated the potential of an agonistic anti-BTLA antibody to treat NTN in vivo .

RESULTS

The BtlaKO mice developed aggravated NTN, driven by an increase of infiltrating renal Th1 cells. Single-cell RNA sequencing showed increased renal T-cell activation and positive regulation of the immune response. Although BTLA-deficient regulatory T cells (Tregs) exhibited preserved suppressive function in vitro and in vivo , BtlaKO T effector cells evaded Treg suppression. Administration of an agonistic anti-BTLA antibody robustly attenuated NTN by suppressing nephritogenic T effector cells and promoting Treg expansion.

CONCLUSIONS

In a model of crescentic GN, BTLA signaling effectively restrained nephritogenic Th1 cells and promoted regulatory T cells. Suppression of T-cell-mediated inflammation by BTLA stimulation may prove relevant for a broad range of conditions involving acute GN.

Plasma proteome of growing tumors

Early detection of cancer is vital for the best chance of successful treatment, but half of all cancers are diagnosed at an advanced stage. A simple and reliable blood screening test applied routinely would therefore address a major unmet medical need. To gain insight into the value of protein biomarkers in early detection and stratification of cancer we determined the time course of changes in the plasma proteome of mice carrying transplanted human lung, breast, colon, or ovarian tumors. For protein measurements we used an aptamer-based assay which simultaneously measures ~ 5000 proteins. Along with tumor lineage-specific biomarkers, we also found 15 markers shared among all cancer types that included the energy metabolism enzymes glyceraldehyde-3-phosphate dehydrogenase, glucose-6-phophate isomerase and dihydrolipoyl dehydrogenase as well as several important biomarkers for maintaining protein, lipid, nucleotide, or carbohydrate balance such as tryptophanyl t-RNA synthetase and nucleoside diphosphate kinase. Using significantly altered proteins in the tumor bearing mice, we developed models to stratify tumor types and to estimate the minimum detectable tumor volume. Finally, we identified significantly enriched common and unique biological pathways among the eight tumor cell lines tested.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Functional Landscape of African Swine Fever Virus-Host and Virus-Virus Protein Interactions

Viral replication fully relies on the host cell machinery, and physical interactions between viral and host proteins mediate key steps of the viral life cycle. Therefore, identifying virus-host protein-protein interactions (PPIs) provides insights into the molecular mechanisms governing virus infection and is crucial for designing novel antiviral strategies. In the case of the African swine fever virus (ASFV), a large DNA virus that causes a deadly panzootic disease in pigs, the limited understanding of host and viral targets hinders the development of effective vaccines and treatments. This review summarizes the current knowledge of virus-host and virus-virus PPIs by collecting and analyzing studies of individual viral proteins. We have compiled a dataset of experimentally determined host and virus protein targets, the molecular mechanisms involved, and the biological functions of the identified virus-host and virus-virus protein interactions during infection. Ultimately, this work provides a comprehensive and systematic overview of ASFV interactome, identifies knowledge gaps, and proposes future research directions.

The human placenta exhibits a unique transcriptomic void

The human placenta exhibits a unique genomic architecture with an unexpectedly high mutation burden and many uniquely expressed genes. The aim of this study is to identify transcripts that are uniquely absent or depleted in the placenta. Here, we show that 40 of 46 of the other organs have no selectively depleted transcripts and that, of the remaining six, the liver has the largest number, with 26. In contrast, the term placenta has 762 depleted transcripts. Gene Ontology analysis of this depleted set highlighted multiple pathways reflecting known unique elements of placental physiology. For example, transcripts associated with neuronal function are in the depleted set-as expected given the lack of placental innervation. However, this demonstrated overrepresentation of genes involved in mitochondrial function (p = 5.8 × 10-10), including PGC-1α, the master regulator of mitochondrial biogenesis, and genes involved in polyamine metabolism (p = 2.1 × 10-4).

Autism-linked UBE3A gain-of-function mutation causes interneuron and behavioral phenotypes when inherited maternally or paternally in mice

The E3 ubiquitin ligase Ube3a is biallelically expressed in neural progenitors and glial cells, suggesting that UBE3A gain-of-function mutations might cause neurodevelopmental disorders irrespective of parent of origin. Here, we engineered a mouse line that harbors an autism-linked UBE3AT485A (T503A in mouse) gain-of-function mutation and evaluated phenotypes in animals that inherited the mutant allele paternally, maternally, or from both parents. We find that paternally and maternally expressed UBE3AT503A results in elevated UBE3A activity in neural progenitors and glial cells. Expression of UBE3AT503A from the maternal allele, but not the paternal one, leads to a persistent elevation of UBE3A activity in neurons. Mutant mice display behavioral phenotypes that differ by parent of origin. Expression of UBE3AT503A, irrespective of its parent of origin, promotes transient embryonic expansion of Zcchc12 lineage interneurons. Phenotypes of Ube3aT503A mice are distinct from Angelman syndrome model mice. Our study has clinical implications for a growing number of disease-linked UBE3A gain-of-function mutations.

Age-related gene expression signatures from limb skeletal muscles and the diaphragm in mice and rats reveal common and species-specific changes

BACKGROUND

As a result of aging, skeletal muscle undergoes atrophy and a decrease in function. This age-related skeletal muscle weakness is known as "sarcopenia". Sarcopenia is part of the frailty observed in humans. In order to discover treatments for sarcopenia, it is necessary to determine appropriate preclinical models and the genes and signaling pathways that change with age in these models.

METHODS AND RESULTS

To understand the changes in gene expression that occur as a result of aging in skeletal muscles, we generated a multi-time-point gene expression signature throughout the lifespan of mice and rats, as these are the most commonly used species in preclinical research and intervention testing. Gastrocnemius, tibialis anterior, soleus, and diaphragm muscles from male and female C57Bl/6J mice and male Sprague Dawley rats were analyzed at ages 6, 12, 18, 21, 24, and 27 months, plus an additional 9-month group was used for rats. More age-related genes were identified in rat skeletal muscles compared with mice; this was consistent with the finding that rat muscles undergo more robust age-related decline in mass. In both species, pathways associated with innate immunity and inflammation linearly increased with age. Pathways linked with extracellular matrix remodeling were also universally downregulated. Interestingly, late downregulated pathways were exclusively found in the rat limb muscles and these were linked to metabolism and mitochondrial respiration; this was not seen in the mouse.

CONCLUSIONS

This extensive, side-by-side transcriptomic profiling shows that the skeletal muscle in rats is impacted more by aging compared with mice, and the pattern of decline in the rat may be more representative of the human. The observed changes point to potential therapeutic interventions to avoid age-related decline in skeletal muscle function.

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataract. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq, and CUT&RUN-seq to discover novel mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Further, we divulge a conserved epigenetic paradigm of cellular differentiation, defined by progressive loss of H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

High-density linkage mapping and genetic dissection of resistance to broomrape (Orobanche crenata Forsk.) in pea (Pisum sativum L.)

Pea (Pisum sativum L.) is a widely cultivated legume of major importance for global food security and agricultural sustainability. Crenate broomrape (Orobanche crenata Forsk.) (Oc) is a parasitic weed severely affecting legumes, including pea, in the Mediterranean Basin and the Middle East. Previously, the identification of the pea line "ROR12", displaying resistance to Oc, was reported. Two-year field trials on a segregant population of 148 F₇ recombinant inbred lines (RILs), originating from a cross between "ROR12" and the susceptible cultivar "Sprinter", revealed high heritability (0.84) of the "ROR12" resistance source. Genotyping-by-sequencing (GBS) on the same RIL population allowed the construction of a high-density pea linkage map, which was compared with the pea reference genome and used for quantitative trait locus (QTL) mapping. Three QTLs associated with the response to Oc infection, named PsOcr-1, PsOcr-2, and PsOcr-3, were identified, with PsOcr-1 explaining 69.3% of the genotypic variance. Evaluation of the effects of different genotypic combinations indicated additivity between PsOcr-1 and PsOcr-2, and between PsOcr-1 and PsOcr-3, and epistasis between PsOcr-2 and PsOcr-3. Finally, three Kompetitive Allele Specific PCR (KASP) marker assays were designed on the single-nucleotide polymorphisms (SNPs) associated with the QTL significance peaks. Besides contributing to the development of pea genomic resources, this work lays the foundation for the obtainment of pea cultivars resistant to Oc and the identification of genes involved in resistance to parasitic Orobanchaceae.

A high-fat diet increases hepatic mitochondrial turnover through restricted acetylation in a NAFLD mouse model

Lysine acetylation of proteins has emerged as a key posttranslational modification (PTM) that regulates mitochondrial metabolism. Acetylation may regulate energy metabolism by inhibiting and affecting the stability of metabolic enzymes and oxidative phosphorylation (OxPhos) subunits. Although protein turnover can be easily measured, due to the low abundance of modified proteins, it has been difficult to evaluate the effect of acetylation on the stability of proteins in vivo. We applied 2H2O-metabolic labeling coupled with immunoaffinity and high-resolution mass spectrometry method to measure the stability of acetylated proteins in mouse liver based on their turnover rates. As a proof-of-concept, we assessed the consequence of high-fat diet (HFD)-induced altered acetylation in protein turnover in LDL receptor-deficient (LDLR-/-) mice susceptible to diet-induced nonalcoholic fatty liver disease (NAFLD). HFD feeding for 12 wk led to steatosis, the early stage of NAFLD. A significant reduction in acetylation of hepatic proteins was observed in NAFLD mice, based on immunoblot analysis and label-free quantification with mass spectrometry. Compared with control mice on a normal diet, NAFLD mice had overall increased turnover rates of hepatic proteins, including mitochondrial metabolic enzymes (0.159 ± 0.079 vs. 0.132 ± 0.068 day-1), suggesting their reduced stability. Also, acetylated proteins had slower turnover rates (increased stability) than native proteins in both groups (0.096 ± 0.056 vs. 0.170 ± 0.059 day-1 in control, and 0.111 ± 0.050 vs. 0.208 ± 0.074 day-1 in NAFLD). Furthermore, association analysis revealed a relationship between the HFD-induced decrease in acetylation and increased turnover rates for hepatic proteins in NAFLD mice. These changes were associated with increased expressions of the hepatic mitochondrial transcriptional factor (TFAM) and complex II subunit without any changes to other OxPhos proteins, suggesting that enhanced mitochondrial biogenesis prevented restricted acetylation-mediated depletion of mitochondrial proteins. We conclude that decreased acetylation of mitochondrial proteins may contribute to adaptive improved hepatic mitochondrial function in the early stages of NAFLD.NEW & NOTEWORTHY This is the first method to quantify acetylome dynamics in vivo. This method revealed acetylation-mediated altered hepatic mitochondrial protein turnover in response to a high-fat diet in a mouse model of NAFLD.

Wound microbiota-mediated correction of matrix metalloproteinase expression promotes re-epithelialization of diabetic wounds

Chronic wounds are a common and costly complication of diabetes, where multifactorial defects contribute to dysregulated skin repair, inflammation, tissue damage, and infection. We previously showed that aspects of the diabetic foot ulcer microbiota were correlated with poor healing outcomes, but many microbial species recovered remain uninvestigated with respect to wound healing. Here we focused on Alcaligenes faecalis , a Gram-negative bacterium that is frequently recovered from chronic wounds but rarely causes infection. Treatment of diabetic wounds with A. faecalis accelerated healing during early stages. We investigated the underlying mechanisms and found that A. faecalis treatment promotes re-epithelialization of diabetic keratinocytes, a process which is necessary for healing but deficient in chronic wounds. Overexpression of matrix metalloproteinases in diabetes contributes to failed epithelialization, and we found that A. faecalis treatment balances this overexpression to allow proper healing. This work uncovers a mechanism of bacterial-driven wound repair and provides a foundation for the development of microbiota-based wound interventions.

Distinct longevity mechanisms across and within species and their association with aging

Lifespan varies within and across species, but the general principles of its control remain unclear. Here, we conducted multi-tissue RNA-seq analyses across 41 mammalian species, identifying longevity signatures and examining their relationship with transcriptomic biomarkers of aging and established lifespan-extending interventions. An integrative analysis uncovered shared longevity mechanisms within and across species, including downregulated Igf1 and upregulated mitochondrial translation genes, and unique features, such as distinct regulation of the innate immune response and cellular respiration. Signatures of long-lived species were positively correlated with age-related changes and enriched for evolutionarily ancient essential genes, involved in proteolysis and PI3K-Akt signaling. Conversely, lifespan-extending interventions counteracted aging patterns and affected younger, mutable genes enriched for energy metabolism. The identified biomarkers revealed longevity interventions, including KU0063794, which extended mouse lifespan and healthspan. Overall, this study uncovers universal and distinct strategies of lifespan regulation within and across species and provides tools for discovering longevity interventions.

Experimental gene expression of developmentally downregulated Crmp1, Crmp4, and Crmp5 promotes axon regeneration and retinal ganglion cell survival after optic nerve injury

Collapsin response mediator proteins (Crmps) play roles in neuronal development and axon growth. However, neuronal-specific roles of Crmp1, Crmp4, and Crmp5 in regeneration of injured central nervous system (CNS) axons in vivo are unclear. Here, we analyzed developmental and subtype-specific expression of Crmp genes in retinal ganglion cells (RGCs), tested whether overexpressing Crmp1, Crmp4, or Crmp5 in RGCs through localized intralocular AAV2 delivery promotes axon regeneration after optic nerve injury in vivo, and characterized developmental co-regulation of gene-concept networks associated with Crmps. We found that all Crmp genes are developmentally downregulated in RGCs during maturation. However, while Crmp1, Crmp2, and Crmp4 were expressed to a varying degree in most RGC subtypes, Crmp3 and Crmp5 were expressed only in a small subset of RGC subtypes. We then found that after optic nerve injury, Crmp1, Crmp4, and Crmp5 promote RGC axon regeneration to varying extents, with Crmp4 promoting the most axon regeneration and also localizing to axons. We also found that Crmp1 and Crmp4, but not Crmp5, promote RGC survival. Finally, we found that Crmp1, Crmp2, Crmp4, and Crmp5's ability to promote axon regeneration is associated with neurodevelopmental mechanisms, which control RGC's intrinsic axon growth capacity.

Evaluating cancer cell line and patient-derived xenograft recapitulation of tumor and non-diseased tissue gene expression profiles in silico

Preclinical models like cancer cell lines and patient-derived xenografts (PDXs) are vital for studying disease mechanisms and evaluating treatment options. It is essential that they accurately recapitulate the disease state of interest to generate results that will translate in the clinic. Prior studies have demonstrated that preclinical models do not recapitulate all biological aspects of human tissues, particularly with respect to the tissue of origin gene expression signatures. Therefore, it is critical to assess how well preclinical model gene expression profiles correlate with human cancer tissues to inform preclinical model selection and data analysis decisions. Here we evaluated how well preclinical models recapitulate human cancer and non-diseased tissue gene expression patterns in silico with respect to the full gene expression profile as well as subsetting by the most variable genes, genes significantly correlated with tumor purity, and tissue-specific genes by using publicly available gene expression profiles across multiple sources. We found that using the full gene set improves correlations between preclinical model and tissue global gene expression profiles, confirmed that GBM PDX global gene expression correlation to GBM tumor global gene expression outperforms GBM cell line to GBM tumor global gene expression correlations, and demonstrated that preclinical models in our study often failed to reproduce tissue-specific expression. While including additional genes for global gene expression comparison between cell lines and tissues decreases the overall correlation, it improves the relative rank between a cell line and its tissue of origin compared to other tissues. Our findings underscore the importance of using the full gene expression set measured when comparing preclinical models and tissues and confirm that tissue-specific patterns are better preserved in GBM PDX models than in GBM cell lines. Future studies can build on these findings to determine the specific pathways and gene sets recapitulated by particular preclinical models to facilitate model selection for a given study design or goal.

Single-cell transcriptomics combined with proteomics of intrathecal IgG reveal transcriptional heterogeneity of oligoclonal IgG-secreting cells in multiple sclerosis

The phenotypes of B lineage cells that produce oligoclonal IgG in multiple sclerosis have not been unequivocally determined. Here, we utilized single-cell RNA-seq data of intrathecal B lineage cells in combination with mass spectrometry of intrathecally synthesized IgG to identify its cellular source. We found that the intrathecally produced IgG matched a larger fraction of clonally expanded antibody-secreting cells compared to singletons. The IgG was traced back to two clonally related clusters of antibody-secreting cells, one comprising highly proliferating cells, and the other consisting of more differentiated cells expressing genes associated with immunoglobulin synthesis. These findings suggest some degree of heterogeneity among cells that produce oligoclonal IgG in multiple sclerosis.

Altered hair root gene expression profiles highlight calcium signaling and lipid metabolism pathways to be associated with curly hair initiation and maintenance in Mangalitza pigs

Hair types have been under strong targeted selection in domestic animals for their impact on skin protection, thermoregulation and exterior morphology, and subsequent economic importance. In pigs, a very special hair phenotype was observed in Mangalitza, who expresses a thick coat of curly bristles and downy hair. Two breed-specific missense variants in TRPM2 and CYP4F3 were suggested to be associated with the Mangalitza pig's hair shape due to their role in hair follicle morphogenesis reported for human and mice. However, the mechanism behind this expression of a curly hair type is still unclear and needs to be explored. In our study, hair shafts were measured and investigated for the curvature of the hair in Mangalitza and crossbreeds in comparison to straight-coated pigs. For molecular studies, hair roots underwent RNA sequencing for a differential gene expression analysis using DESeq2. The output matrix of normalized counts was then used to construct weighted gene co-expression networks. The resulting hair root gene expression profiles highlighted 454 genes to be significantly differentially expressed for initiation of curly hair phenotype in newborn Mangalitza piglets versus post-initiation in later development. Furthermore, 2,554 genes showed a significant differential gene expression in curly hair in comparison to straight hair. Neither TRPM2 nor CYP4F3 were identified as differentially expressed. Incidence of the genes in weighted co-expression networks associated with TRPM2 and CYP4F3, and prominent interactions of subsequent proteins with lipids and calcium-related pathways suggested calcium signaling and/or lipid metabolism as essential players in the induction of the curly hair as well as an ionic calcium-dependency to be a prominent factor for the maintenance of this phenotype. Subsequently, our study highlights the complex interrelations and dependencies of mutant genes TRPM2 and CYP4F3 and associated gene expression patterns, allowing the initiation of curly hair type during the development of a piglet as well as the maintenance in adult individuals.

Direct and indirect consequences of PAB1 deletion in the regulation of translation initiation, translation termination, and mRNA decay

Cytoplasmic poly(A)-binding protein (PABPC; Pab1 in yeast) is thought to be involved in multiple steps of post-transcriptional control, including translation initiation, translation termination, and mRNA decay. To understand these roles of PABPC in more detail for endogenous mRNAs, and to distinguish its direct effects from indirect effects, we have employed RNA-Seq and Ribo-Seq to analyze changes in the abundance and translation of the yeast transcriptome, as well as mass spectrometry to assess the abundance of the components of the yeast proteome, in cells lacking the PAB1 gene. We observed drastic changes in the transcriptome and proteome, as well as defects in translation initiation and termination, in pab1Δ cells. Defects in translation initiation and the stabilization of specific classes of mRNAs in pab1Δ cells appear to be partly indirect consequences of reduced levels of specific initiation factors, decapping activators, and components of the deadenylation complex in addition to the general loss of Pab1's direct role in these processes. Cells devoid of Pab1 also manifested a nonsense codon readthrough phenotype indicative of a defect in translation termination, but this defect may be a direct effect of the loss of Pab1 as it could not be attributed to significant reductions in the levels of release factors.

Myelin dysfunction drives amyloid-β deposition in models of Alzheimer's disease

The incidence of Alzheimer's disease (AD), the leading cause of dementia, increases rapidly with age, but why age constitutes the main risk factor is still poorly understood. Brain ageing affects oligodendrocytes and the structural integrity of myelin sheaths1, the latter of which is associated with secondary neuroinflammation2,3. As oligodendrocytes support axonal energy metabolism and neuronal health4-7, we hypothesized that loss of myelin integrity could be an upstream risk factor for neuronal amyloid-β (Aβ) deposition, the central neuropathological hallmark of AD. Here we identify genetic pathways of myelin dysfunction and demyelinating injuries as potent drivers of amyloid deposition in mouse models of AD. Mechanistically, myelin dysfunction causes the accumulation of the Aβ-producing machinery within axonal swellings and increases the cleavage of cortical amyloid precursor protein. Suprisingly, AD mice with dysfunctional myelin lack plaque-corralling microglia despite an overall increase in their numbers. Bulk and single-cell transcriptomics of AD mouse models with myelin defects show that there is a concomitant induction of highly similar but distinct disease-associated microglia signatures specific to myelin damage and amyloid plaques, respectively. Despite successful induction, amyloid disease-associated microglia (DAM) that usually clear amyloid plaques are apparently distracted to nearby myelin damage. Our data suggest a working model whereby age-dependent structural defects of myelin promote Aβ plaque formation directly and indirectly and are therefore an upstream AD risk factor. Improving oligodendrocyte health and myelin integrity could be a promising target to delay development and slow progression of AD.

Vitamin A deficiency compromises the barrier function of the retinal pigment epithelium

A major cause for childhood blindness worldwide is attributed to nutritional vitamin A deficiency. Surprisingly, the molecular basis of the ensuing retinal degeneration has not been well defined. Abundant expression of the retinoid transporter STRA6 in the retinal pigment epithelium (RPE) and homeostatic blood levels of retinol-binding protein delay vitamin A deprivation of the mouse eyes. Hence, genetic dissection of STRA6 makes mice susceptible to nutritional manipulation of ocular retinoid status. We performed RNA-seq analyses and complemented the data with tests of visual physiology, ocular morphology, and retinoid biochemistry to compare eyes with different vitamin A status. Mild ocular vitamin A deficiency decreased transcripts of photoreceptor transduction pathway-related genes and increased transcripts of oxidative stress pathways. The response was associated with impaired visual sensitivity and an accumulation of fluorescent debris in the retina. Severe vitamin A deficiency did not only impair visual perception but also decreased transcripts of genes encoding cell adhesion and cellular junction proteins. This response altered cell morphology, resulted in significant changes in transport pathways of small molecules, and compromised the barrier function of the RPE. Together, our analyses characterize the molecular events underlying nutritional blindness in a novel mouse model and indicate that breakdown of the outer blood-retinal barrier contributes to retinal degeneration and photoreceptor cell death in severe vitamin A deficiency.

Mass-Spectrometric Evaluation of the African Swine Fever Virus-Induced Host Shutoff Using Dynamic Stable Isotope Labeling with Amino Acids in Cell Culture (SILAC)

African swine fever is a viral disease of swine caused by the African swine fever virus (ASFV). Currently, ASFV is spreading over the Eurasian continent and threatening global pig husbandry. One viral strategy to undermine an efficient host cell response is to establish a global shutoff of host protein synthesis. This shutoff has been observed in ASFV-infected cultured cells using two-dimensional electrophoresis combined with metabolic radioactive labeling. However, it remained unclear if this shutoff was selective for certain host proteins. Here, we characterized ASFV-induced shutoff in porcine macrophages by measurement of relative protein synthesis rates using a mass spectrometric approach based on stable isotope labeling with amino acids in cell culture (SILAC). The impact of ASFV infection on the synthesis of >2000 individual host proteins showed a high degree of variability, ranging from complete shutoff to a strong induction of proteins that are absent from naïve cells. GO-term enrichment analysis revealed that the most effective shutoff was observed for proteins related to RNA metabolism, while typical representatives of the innate immune system were strongly induced after infection. This experimental setup is suitable to quantify a virion-induced host shutoff (vhs) after infection with different viruses.

Cancer-specific functional profiling in microsatellite-unstable (MSI) colon and endometrial cancers using combined differentially expressed genes and biclustering analysis

Microsatellite-unstable (MSI) cancers have distinct genetic and clinical features from microsatellite-stable cancers, but the molecular functional differences between MSI cancers originating from different tissues or organs have not been well studied because the application of usual differentially expressed gene (DEG) analysis is error-prone, producing too many noncancer-specific normally functioning genes. To maximize therapeutic efficacy, biomarkers reflecting cancer-specific differences between MSI cancers of different tissue origins should be identified. To identify functional differences between MSI colon and endometrial cancers, we combined DEG analysis and biclustering instead of DEG analysis alone and refined functionally relevant biclusters reflecting genuine functional differences between the 2 tumors. Specifically, using The Cancer Genome Atlas and genome-tissue expression as data sources, gene ontology (GO) enrichment tests were performed after routinely identifying DEGs between the 2 tumors with the exclusion of DEGs identified in their normal counterparts. Cancer-specific biclusters and associated enriched GO terms were obtained by biclustering with enrichment tests for the preferences for cancer type (either colon or endometrium) and GO enrichment tests for each cancer-specific bicluster, respectively. A novel childness score was developed to select functionally relevant biclusters among cancer-specific biclusters based on the extent to which the enriched GO terms of the biclusters tended to be child terms of the enriched GO terms in DEGs. The selected biclusters were tested using survival analysis to validate their clinical significance. We performed multiple sequential analyses to produce functionally relevant biclusters from the RNA sequencing data of MSI colon and endometrial cancer samples and their normal counterparts. We identified 3066 cancer-specific DEGs. Biclustering analysis revealed 153 biclusters and 41 cancer-specific biclusters were selected using Fisher exact test. A mean childness score over 0.6 was applied as the threshold and yielded 8 functionally relevant biclusters from cancer-specific biclusters. Functional differences appear to include gland cavitation and the TGF-β receptor, G protein, and cytokine pathways. In the survival analysis, 6 of the 8 functionally relevant biclusters were statistically significant. By attenuating noise and applying a synergistic contribution of DEG results, we refined candidate biomarkers to complement tissue-specific features of MSI tumors.

Cross-disorder genetic analysis of immune diseases reveals distinct gene associations that converge on common pathways

Genome-wide association studies (GWAS) have mapped thousands of susceptibility loci associated with immune-mediated diseases. To assess the extent of the genetic sharing across nine immune-mediated diseases we apply genomic structural equation modelling to GWAS data from European populations. We identify three disease groups: gastrointestinal tract diseases, rheumatic and systemic diseases, and allergic diseases. Although loci associated with the disease groups are highly specific, they converge on perturbing the same pathways. Finally, we test for colocalization between loci and single-cell eQTLs derived from peripheral blood mononuclear cells. We identify the causal route by which 46 loci predispose to three disease groups and find evidence for eight genes being candidates for drug repurposing. Taken together, here we show that different constellations of diseases have distinct patterns of genetic associations, but that associated loci converge on perturbing different nodes in T cell activation and signalling pathways.

Genome graphs detect human polymorphisms in active epigenomic state during influenza infection

Genetic variants, including mobile element insertions (MEIs), are known to impact the epigenome. We hypothesized that genome graphs, which encapsulate genetic diversity, could reveal missing epigenomic signals. To test this, we sequenced the epigenome of monocyte-derived macrophages from 35 ancestrally diverse individuals before and after influenza infection, allowing us to investigate the role of MEIs in immunity. We characterized genetic variants and MEIs using linked reads and built a genome graph. Mapping epigenetic data revealed 2.3%-3% novel peaks for H3K4me1, H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq), and ATAC-seq. Additionally, the use of a genome graph modified some quantitative trait loci estimates and revealed 375 polymorphic MEIs in an active epigenomic state. Among these is an AluYh3 polymorphism whose chromatin state changed after infection and was associated with the expression of TRIM25, a gene that restricts influenza RNA synthesis. Our results demonstrate that graph genomes can reveal regulatory regions that would have been overlooked by other approaches.

g:Profiler-interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update)

g:Profiler is a reliable and up-to-date functional enrichment analysis tool that supports various evidence types, identifier types and organisms. The toolset integrates many databases, including Gene Ontology, KEGG and TRANSFAC, to provide a comprehensive and in-depth analysis of gene lists. It also provides interactive and intuitive user interfaces and supports ordered queries and custom statistical backgrounds, among other settings. g:Profiler provides multiple programmatic interfaces to access its functionality. These can be easily integrated into custom workflows and external tools, making them valuable resources for researchers who want to develop their own solutions. g:Profiler has been available since 2007 and is used to analyse millions of queries. Research reproducibility and transparency are achieved by maintaining working versions of all past database releases since 2015. g:Profiler supports 849 species, including vertebrates, plants, fungi, insects and parasites, and can analyse any organism through user-uploaded custom annotation files. In this update article, we introduce a novel filtering method highlighting Gene Ontology driver terms, accompanied by new graph visualizations providing a broader context for significant Gene Ontology terms. As a leading enrichment analysis and gene list interoperability service, g:Profiler offers a valuable resource for genetics, biology and medical researchers. It is freely accessible at https://biit.cs.ut.ee/gprofiler.

GTAC enables parallel genotyping of multiple genomic loci with chromatin accessibility profiling in single cells

Understanding clonal evolution and cancer development requires experimental approaches for characterizing the consequences of somatic mutations on gene regulation. However, no methods currently exist that efficiently link high-content chromatin accessibility with high-confidence genotyping in single cells. To address this, we developed Genotyping with the Assay for Transposase-Accessible Chromatin (GTAC), enabling accurate mutation detection at multiple amplified loci, coupled with robust chromatin accessibility readout. We applied GTAC to primary acute myeloid leukemia, obtaining high-quality chromatin accessibility profiles and clonal identities for multiple mutations in 88% of cells. We traced chromatin variation throughout clonal evolution, showing the restriction of different clones to distinct differentiation stages. Furthermore, we identified switches in transcription factor motif accessibility associated with a specific combination of driver mutations, which biased transformed progenitors toward a leukemia stem cell-like chromatin state. GTAC is a powerful tool to study clonal heterogeneity across a wide spectrum of pre-malignant and neoplastic conditions.

The Fgf/Erf/NCoR1/2 repressive axis controls trophoblast cell fate

Placental development relies on coordinated cell fate decisions governed by signalling inputs. However, little is known about how signalling cues are transformed into repressive mechanisms triggering lineage-specific transcriptional signatures. Here, we demonstrate that upon inhibition of the Fgf/Erk pathway in mouse trophoblast stem cells (TSCs), the Ets2 repressor factor (Erf) interacts with the Nuclear Receptor Co-Repressor Complex 1 and 2 (NCoR1/2) and recruits it to key trophoblast genes. Genetic ablation of Erf or Tbl1x (a component of the NCoR1/2 complex) abrogates the Erf/NCoR1/2 interaction. This leads to mis-expression of Erf/NCoR1/2 target genes, resulting in a TSC differentiation defect. Mechanistically, Erf regulates expression of these genes by recruiting the NCoR1/2 complex and decommissioning their H3K27ac-dependent enhancers. Our findings uncover how the Fgf/Erf/NCoR1/2 repressive axis governs cell fate and placental development, providing a paradigm for Fgf-mediated transcriptional control.

Patterns of Oncogene Coexpression at Single-Cell Resolution Influence Survival in Lymphoma

Cancers often overexpress multiple clinically relevant oncogenes, but it is not known if combinations of oncogenes in cellular subpopulations within a cancer influence clinical outcomes. Using quantitative multispectral imaging of the prognostically relevant oncogenes MYC, BCL2, and BCL6 in diffuse large B-cell lymphoma (DLBCL), we show that the percentage of cells with a unique combination MYC+BCL2+BCL6- (M+2+6-) consistently predicts survival across four independent cohorts (n = 449), an effect not observed with other combinations including M+2+6+. We show that the M+2+6- percentage can be mathematically derived from quantitative measurements of the individual oncogenes and correlates with survival in IHC (n = 316) and gene expression (n = 2,521) datasets. Comparative bulk/single-cell transcriptomic analyses of DLBCL samples and MYC/BCL2/BCL6-transformed primary B cells identify molecular features, including cyclin D2 and PI3K/AKT as candidate regulators of M+2+6- unfavorable biology. Similar analyses evaluating oncogenic combinations at single-cell resolution in other cancers may facilitate an understanding of cancer evolution and therapy resistance.

SIGNIFICANCE

Using single-cell-resolved multiplexed imaging, we show that selected subpopulations of cells expressing specific combinations of oncogenes influence clinical outcomes in lymphoma. We describe a probabilistic metric for the estimation of cellular oncogenic coexpression from IHC or bulk transcriptomes, with possible implications for prognostication and therapeutic target discovery in cancer. This article is highlighted in the In This Issue feature, p. 1027.