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

Genetic background of walking ability and its relationship with leg defects, mortality, and performance traits in turkeys (Meleagris gallopavo)

This study aimed to explore the genetic basis of walking ability and potentially related performance traits in turkey purebred populations. Phenotypic, pedigree, and genomic datasets from 2 turkey lines hatched between 2010 and 2023 were included in the study. Walking ability data, defined based on a scoring system ranging from 1 (worst) to 6 (best), were collected on 192,019 animals of a female line and 235,461 animals of a male line. Genomic information was obtained for 46,427 turkeys (22,302 from a female line and 24,125 from a male line) using a 65K single nucleotide polymorphism (SNP) panel. Variance components and heritability for walking ability were estimated. Furthermore, genetic and phenotypic correlations among walking ability, mortality and disorders, and performance traits were calculated. A genome-wide association study (GWAS) was also conducted to identify SNPs associated with walking ability. Walking ability is moderately heritable (0.23 ± 0.01) in both turkey lines. The genetic correlations between walking ability and the other evaluated traits ranged from -0.02 to -0.78, with leg defects exhibiting the strongest negative correlation with walking ability. In the female line, 31 SNPs were associated with walking ability and overlapped with 116 genes. These positional genes are linked to 6 gene ontology (GO) terms. Notably, genes such as CSRP2, DDX1, RHBDL1, SEZ6L, and CTSK are involved in growth, development, locomotion, and bone disorders. GO terms, including fibronectin binding (GO:0001968), peptide cross-linking (GO:0018149), and catabolic process (GO:0009057), are directly linked with mobility. In the male line, 66 markers associated with walking ability were identified and overlapped with 281 genes. These genes are linked to 12 GO terms. Genes such as RB1CC1, TNNI1, MSTN, FN1, SIK3, PADI2, ERBB4, B3GNT2, and BACE1 are associated with cell growth, myostatin development, and disorders. GO terms in the male line are predominantly related to lipid metabolism. In conclusion, walking ability is moderately heritable in both populations. Furthermore, walking ability can be enhanced through targeted genetic selection, emphasizing its relevance to both animal welfare and productivity.

Alcaligenes faecalis corrects aberrant matrix metalloproteinase expression to promote reepithelialization 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 reepithelialization of diabetic keratinocytes, a process that 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.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

The ethylene-responsive transcription factor ERF024 is a novel regulator of climacteric fruit ripening in melon

Fruit ripening is an essential developmental stage in Angiosperms triggered by hormonal signals such as ethylene, a major player in climacteric ripening. Melon is a unique crop showing both climacteric and non-climacteric cultivars, offering an ideal model for dissecting the genetic mechanisms underpinning this process. The major quantitative trait locus ETHQV8.1 was previously identified as a key regulator of melon fruit ripening. Here, we narrowed down ETHQV8.1 to a precise genomic region containing a single gene, the transcription factor CmERF024. Functional validation using CRISPR/Cas9 knock-out plants unequivocally identified CmERF024 as the causal gene governing ETHQV8.1. The erf024 mutants exhibited suppression of ethylene production, leading to a significant delay and attenuation of fruit ripening. Integrative multi-omic analyses encompassing RNA-seq, DAP-seq, and DNase-seq revealed the association of CmERF024 with chromatin accessibility and gene expression dynamics throughout fruit ripening. Our data suggest CmERF024 as a novel regulator of climacteric fruit ripening in melon.

The hepatocyte epidermal growth factor receptor (EGFR) pathway regulates the cellular interactome within the liver fibrotic niche

Liver fibrosis is the consequence of chronic liver injury in the presence of an inflammatory component. Although the main executors of this activation are known, the mechanisms that lead to the inflammatory process that mediates the production of pro-fibrotic factors are not well characterized. Epidermal growth factor receptor (EGFR) signaling in hepatocytes is essential for the regenerative processes of the liver; however, its potential role in regulating the fibrotic niche is not yet clear. Our group generated a mouse model that expresses an inactive truncated form of the EGFR specifically in hepatocytes (ΔEGFR mice). Here, we have analyzed the response of WT and ΔEGFR mice to chronic treatment with carbon tetrachloride (CCl4), which induces a pro-inflammatory and fibrotic process in the liver. The results indicated that the hallmarks of liver fibrosis were attenuated in CCl4-treated ΔEGFR mice when compared with CCl4-treated WT mice, coinciding with a faster resolution of the fibrotic process and ameliorated damage. The absence of EGFR activity in hepatocytes induced changes in the pattern of immune cells in the liver, with a notable increase in the population of M2 macrophages, more related to fibrosis resolution, as well as in the population of lymphocytes related to eradication of the damage. Transcriptome analysis of hepatocytes, and secretome studies of extracellular media from in vitro experiments, allowed us to elucidate the specific molecular mechanisms regulated by EGFR that mediate hepatocyte production of both pro-fibrotic and pro-inflammatory mediators; these have consequences for the deposition of extracellular matrix proteins, as well as for the immune microenvironment. Overall, our study uncovered novel mechanistic insights regarding EGFR kinase-dependent actions in hepatocytes that reveal its key role in chronic liver damage. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Whole-genome SNP allele frequency differences between Tibetan and Large white pigs reveal genes associated with skeletal muscle growth

BACKGROUND

The skeletal muscle growth rate and body size of Tibetan pigs (TIB) are lower than Large white pigs (LW). However, the underlying genetic basis attributing to these differences remains uncertain. To address this knowledge gap, the present study employed whole-genome sequencing of TIB (slow growth) and LW (fast growth) individuals, and integrated with existing NCBI sequencing datasets of TIB and LW individuals, enabling the identification of a comprehensive set of genetic variations for each breed. The specific and predominant SNPs in the TIB and LW populations were detected by using a cutoff value of 0.50 for SNP allele frequency and absolute allele frequency differences (△AF) between the TIB and LW populations.

RESULTS

A total of 21,767,938 SNPs were retrieved from 44 TIB and 29 LW genomes. The analysis detected 2,893,106 (13.29%) and 813,310 (3.74%) specific and predominant SNPs in the TIB and LW populations, and annotated to 24,560 genes. Further GO analysis revealed 291 genes involved in biological processes related to striated and/or skeletal muscle differentiation, proliferation, hypertrophy, regulation of striated muscle cell differentiation and proliferation, and myoblast differentiation and fusion. These 291 genes included crucial regulators of muscle cell determination, proliferation, differentiation, and hypertrophy, such as members of the Myogenic regulatory factors (MRF) (MYOD, MYF5, MYOG, MYF6) and Myocyte enhancer factor 2 (MEF2) (MEF2A, MEF2C, MEF2D) families, as well as muscle growth inhibitors (MSTN, ACVR1, and SMAD1); KEGG pathway analysis revealed 106 and 20 genes were found in muscle growth related positive and negative regulatory signaling pathways. Notably, genes critical for protein synthesis, such as MTOR, IGF1, IGF1R, IRS1, INSR, and RPS6KA6, were implicated in these pathways.

CONCLUSION

This study employed an effective methodology to rigorously identify the potential genes associated with skeletal muscle development. A substantial number of SNPs and genes that potentially play roles in the divergence observed in skeletal muscle growth between the TIB and LW breeds were identified. These findings offer valuable insights into the genetic underpinnings of skeletal muscle development and present opportunities for enhancing meat production through pig breeding.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Minor Spliceosomal 65K/RNPC3 Interacts with ANKRD11 and Mediates HDAC3-Regulated Histone Deacetylation and Transcription

RNA splicing is crucial in the multilayer regulatory networks for gene expression, making functional interactions with DNA- and other RNA-processing machineries in the nucleus. However, these established couplings are all major spliceosome-related; whether the minor spliceosome is involved remains unclear. Here, through affinity purification using Drosophila lysates, an interaction is identified between the minor spliceosomal 65K/RNPC3 and ANKRD11, a cofactor of histone deacetylase 3 (HDAC3). Using a CRISPR/Cas9 system, Deletion strains are constructed and found that both Dm65KΔ/Δ and Dmankrd11Δ/Δ mutants have reduced histone deacetylation at Lys9 of histone H3 (H3K9) and Lys5 of histone H4 (H4K5) in their heads, exhibiting various neural-related defects. The 65K-ANKRD11 interaction is also conserved in human cells, and the HsANKRD11 middle-uncharacterized domain mediates Hs65K association with HDAC3. Cleavage under targets and tagmentation (CUT&Tag) assays revealed that HsANKRD11 is a bridging factor, which facilitates the synergistic common chromatin-binding of HDAC3 and Hs65K. Knockdown (KD) of HsANKRD11 simultaneously decreased their common binding, resulting in reduced deacetylation of nearby H3K9. Ultimately, this study demonstrates that expression changes of many genes caused by HsANKRD11-KD are due to the decreased common chromatin-binding of HDAC3 and Hs65K and subsequently reduced deacetylation of H3K9, illustrating a novel and conserved coupling mechanism that links the histone deacetylation with minor spliceosome for the regulation of gene expression.

MBNL1 Regulates Programmed Postnatal Switching Between Regenerative and Differentiated Cardiac States

BACKGROUND

Discovering determinants of cardiomyocyte maturity is critical for deeply understanding the maintenance of differentiated states and potentially reawakening endogenous regenerative programs in adult mammalian hearts as a therapeutic strategy. Forced dedifferentiation paired with oncogene expression is sufficient to drive cardiac regeneration, but elucidation of endogenous developmental regulators of the switch between regenerative and mature cardiomyocyte cell states is necessary for optimal design of regenerative approaches for heart disease. MBNL1 (muscleblind-like 1) regulates fibroblast, thymocyte, and erythroid differentiation and proliferation. Hence, we examined whether MBNL1 promotes and maintains mature cardiomyocyte states while antagonizing cardiomyocyte proliferation.

METHODS

MBNL1 gain- and loss-of-function mouse models were studied at several developmental time points and in surgical models of heart regeneration. Multi-omics approaches were combined with biochemical, histological, and in vitro assays to determine the mechanisms through which MBNL1 exerts its effects.

RESULTS

MBNL1 is coexpressed with a maturation-association genetic program in the heart and is regulated by the MEIS1/calcineurin signaling axis. Targeted MBNL1 overexpression early in development prematurely transitioned cardiomyocytes to hypertrophic growth, hypoplasia, and dysfunction, whereas loss of MBNL1 function increased cardiomyocyte cell cycle entry and proliferation through altered cell cycle inhibitor transcript stability. Moreover, MBNL1-dependent stabilization of estrogen-related receptor signaling was essential for maintaining cardiomyocyte maturity in adult myocytes. In accordance with these data, modulating MBNL1 dose tuned the temporal window of neonatal cardiac regeneration, where increased MBNL1 expression arrested myocyte proliferation and regeneration and MBNL1 deletion promoted regenerative states with prolonged myocyte proliferation. However, MBNL1 deficiency was insufficient to promote regeneration in the adult heart because of cell cycle checkpoint activation.

CONCLUSIONS

Here, MBNL1 was identified as an essential regulator of cardiomyocyte differentiated states, their developmental switch from hyperplastic to hypertrophic growth, and their regenerative potential through controlling an entire maturation program by stabilizing adult myocyte mRNAs during postnatal development and throughout adulthood. Targeting loss of cardiomyocyte maturity and downregulation of cell cycle inhibitors through MBNL1 deletion was not sufficient to promote adult regeneration.

Neuronal activity-driven O-GlcNAcylation promotes mitochondrial plasticity

Neuronal activity is an energy-intensive process that is largely sustained by instantaneous fuel utilization and ATP synthesis. However, how neurons couple ATP synthesis rate to fuel availability is largely unknown. Here, we demonstrate that the metabolic sensor enzyme O-linked N-acetyl glucosamine (O-GlcNAc) transferase regulates neuronal activity-driven mitochondrial bioenergetics in hippocampal and cortical neurons. We show that neuronal activity upregulates O-GlcNAcylation in mitochondria. Mitochondrial O-GlcNAcylation is promoted by activity-driven glucose consumption, which allows neurons to compensate for high energy expenditure based on fuel availability. To determine the proteins that are responsible for these adjustments, we mapped the mitochondrial O-GlcNAcome of neurons. Finally, we determine that neurons fail to meet activity-driven metabolic demand when O-GlcNAcylation dynamics are prevented. Our findings suggest that O-GlcNAcylation provides a fuel-dependent feedforward control mechanism in neurons to optimize mitochondrial performance based on neuronal activity. This mechanism thereby couples neuronal metabolism to mitochondrial bioenergetics and plays a key role in sustaining energy homeostasis.

Pulmonary maternal immune activation does not cross the placenta but leads to fetal metabolic adaptation

The fetal development of organs and functions is vulnerable to perturbation by maternal inflammation which may increase susceptibility to disorders after birth. Because it is not well understood how the placenta and fetus respond to acute lung- inflammation, we characterize the response to maternal pulmonary lipopolysaccharide exposure across 24 h in maternal and fetal organs using multi-omics, imaging and integrative analyses. Unlike maternal organs, which mount strong inflammatory immune responses, the placenta upregulates immuno-modulatory genes, in particular the IL-6 signaling suppressor Socs3. Similarly, we observe no immune response in the fetal liver, which instead displays metabolic changes, including increases in lipids containing docosahexaenoic acid, crucial for fetal brain development. The maternal liver and plasma display similar metabolic alterations, potentially increasing bioavailability of docosahexaenoic acid for the mother and fetus. Thus, our integrated temporal analysis shows that systemic inflammation in the mother leads to a metabolic perturbation in the fetus.

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

High fat diet ameliorates mitochondrial cardiomyopathy in CHCHD10 mutant mice

Mutations in CHCHD10, a mitochondrial protein with undefined functions, are associated with autosomal dominant mitochondrial diseases. Chchd10 knock-in mice harboring a heterozygous S55L mutation (equivalent to human pathogenic S59L) develop a fatal mitochondrial cardiomyopathy caused by CHCHD10 aggregation and proteotoxic mitochondrial integrated stress response (mtISR). In mutant hearts, mtISR is accompanied by a metabolic rewiring characterized by increased reliance on glycolysis rather than fatty acid oxidation. To counteract this metabolic rewiring, heterozygous S55L mice were subjected to chronic high-fat diet (HFD) to decrease insulin sensitivity and glucose uptake and enhance fatty acid utilization in the heart. HFD ameliorated the ventricular dysfunction of mutant hearts and significantly extended the survival of mutant female mice affected by severe pregnancy-induced cardiomyopathy. Gene expression profiles confirmed that HFD increased fatty acid utilization and ameliorated cardiomyopathy markers. Importantly, HFD also decreased accumulation of aggregated CHCHD10 in the S55L heart, suggesting activation of quality control mechanisms. Overall, our findings indicate that metabolic therapy can be effective in mitochondrial cardiomyopathies associated with proteotoxic stress.

Single nucleus and spatial transcriptomic profiling of healthy human hamstring tendon

The molecular and cellular basis of health in human tendons remains poorly understood. Among human tendons, hamstring tendon has markedly low pathology and can provide a prototypic healthy tendon reference. The aim of this study was to determine the transcriptomes and location of all cell types in healthy hamstring tendon. Using single nucleus RNA sequencing, we profiled the transcriptomes of 10 533 nuclei from four healthy donors and identified 12 distinct cell types. We confirmed the presence of two fibroblast cell types, endothelial cells, mural cells, and immune cells, and identified cell types previously unreported in tendons, including different skeletal muscle cell types, satellite cells, adipocytes, and undefined nervous system cells. The location of these cell types within tendon was defined using spatial transcriptomics and imaging, and potential transcriptional networks and cell-cell interactions were analyzed. We demonstrate that fibroblasts have the highest number of potential cell-cell interactions in our dataset, are present throughout the tendon, and play an important role in the production and organization of extracellular matrix, thus confirming their role as key regulators of hamstring tendon homeostasis. Overall, our findings underscore the complexity of the cellular networks that underpin healthy human tendon function and the central role of fibroblasts as key regulators of hamstring tendon tissue homeostasis.

Small-molecule inhibition of the METTL3/METTL14 complex suppresses neuroblastoma tumor growth and promotes differentiation

The N6-methyladenosine (m6A) RNA modification is an important regulator of gene expression. m6A is deposited by a methyltransferase complex that includes methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14). High levels of METTL3/METTL14 drive the growth of many types of adult cancer, and METTL3/METTL14 inhibitors are emerging as new anticancer agents. However, little is known about the m6A epitranscriptome or the role of the METTL3/METTL14 complex in neuroblastoma, a common pediatric cancer. Here, we show that METTL3 knockdown or pharmacologic inhibition with the small molecule STM2457 leads to reduced neuroblastoma cell proliferation and increased differentiation. These changes in neuroblastoma phenotype are associated with decreased m6A deposition on transcripts involved in nervous system development and neuronal differentiation, with increased stability of target mRNAs. In preclinical studies, STM2457 treatment suppresses the growth of neuroblastoma tumors in vivo. Together, these results support the potential of METTL3/METTL14 complex inhibition as a therapeutic strategy against neuroblastoma.

Meta-analysis towards FSHD reveals misregulation of neuromuscular junction, nuclear envelope, and spliceosome

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most common autosomal dominant muscle disorders, yet no cure or amelioration exists. The clinical presentation is diverse, making it difficult to identify the actual driving pathomechanism among many downstream events. To unravel this complexity, we performed a meta-analysis of 13 original omics datasets (in total 171 FSHD and 129 control samples). Our approach confirmed previous findings about the disease pathology and specified them further. We confirmed increased expression of former proposed DUX4 biomarkers, and furthermore impairment of the respiratory chain. Notably, the meta-analysis provides insights about so far not reported pathways, including misregulation of neuromuscular junction protein encoding genes, downregulation of the spliceosome, and extensive alterations of nuclear envelope protein expression. Finally, we developed a publicly available shiny app to provide a platform for researchers who want to search our analysis for genes of interest in the future.

Single-cell multi-cohort dissection of the schizophrenia transcriptome

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

Global, site-resolved analysis of ubiquitylation occupancy and turnover rate reveals systems properties

Ubiquitylation regulates most proteins and biological processes in a eukaryotic cell. However, the site-specific occupancy (stoichiometry) and turnover rate of ubiquitylation have not been quantified. Here we present an integrated picture of the global ubiquitylation site occupancy and half-life. Ubiquitylation site occupancy spans over four orders of magnitude, but the median ubiquitylation site occupancy is three orders of magnitude lower than that of phosphorylation. The occupancy, turnover rate, and regulation of sites by proteasome inhibitors are strongly interrelated, and these attributes distinguish sites involved in proteasomal degradation and cellular signaling. Sites in structured protein regions exhibit longer half-lives and stronger upregulation by proteasome inhibitors than sites in unstructured regions. Importantly, we discovered a surveillance mechanism that rapidly and site-indiscriminately deubiquitylates all ubiquitin-specific E1 and E2 enzymes, protecting them against accumulation of bystander ubiquitylation. The work provides a systems-scale, quantitative view of ubiquitylation properties and reveals general principles of ubiquitylation-dependent governance.

Spiral ganglion neuron degeneration in aminoglycoside-deafened rats involves innate and adaptive immune responses not requiring complement

Spiral ganglion neurons (SGNs) transmit auditory information from cochlear hair cells to the brain. SGNs are thus not only important for normal hearing, but also for effective functioning of cochlear implants, which stimulate SGNs when hair cells are missing. SGNs slowly degenerate following aminoglycoside-induced hair cell loss, a process thought to involve an immune response. However, the specific immune response pathways involved remain unknown. We used RNAseq to gain a deeper understanding immune-related and other transcriptomic changes that occur in the rat spiral ganglion after kanamycin-induced deafening. Among the immune and inflammatory genes that were selectively upregulated in deafened spiral ganglia, the complement cascade genes were prominent. We then assessed SGN survival, as well as immune cell numbers and activation, in the spiral ganglia of rats with a CRISPR-Cas9-mediated knockout of complement component 3 (C3). Similar to previous findings in our lab and other deafened rodent models, we observed an increase in macrophage number and increased expression of CD68, a marker of phagocytic activity and cell activation, in macrophages in the deafened ganglia. Moreover, we found an increase in MHCII expression on spiral ganglion macrophages and an increase in lymphocyte number in the deafened ganglia, suggestive of an adaptive immune response. However, C3 knockout did not affect SGN survival or increase in macrophage number/activation, implying that complement activation does not play a role in SGN death after deafening. Together, these data suggest that both innate and adaptive immune responses are activated in the deafened spiral ganglion, with the adaptive response directly contributing to cochlear neurodegeneration.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Syncytiotrophoblast Markers Are Downregulated in Placentas from Idiopathic Stillbirths

The trophoblast cells are responsible for the transfer of nutrients between the mother and the foetus and play a major role in placental endocrine function by producing and releasing large amounts of hormones and growth factors. Syncytiotrophoblast cells (STB), formed by the fusion of mononuclear cytotrophoblasts (CTB), constitute the interface between the foetus and the mother and are essential for all of these functions. We performed transcriptome analysis of human placental samples from two control groups-live births (LB), and stillbirths (SB) with a clinically recognised cause-and from our study group, idiopathic stillbirths (iSB). We identified 1172 DEGs in iSB, when comparing with the LB group; however, when we compared iSB with the SB group, only 15 and 12 genes were down- and upregulated in iSB, respectively. An assessment of these DEGs identified 15 commonly downregulated genes in iSB. Among these, several syncytiotrophoblast markers, like genes from the PSG and CSH families, as well as ALPP, KISS1, and CRH, were significantly downregulated in placental samples from iSB. The transcriptome analysis revealed underlying differences at a molecular level involving the syncytiotrophoblast. This suggests that defects in the syncytial layer may underlie unexplained stillbirths, therefore offering insights to improve clinical obstetrics practice.

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

The ch12q13 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 cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

MultiRNAflow: integrated analysis of temporal RNA-seq data with multiple biological conditions

MOTIVATION

The dynamic transcriptional mechanisms that govern eukaryotic cell function can now be analyzed by RNA sequencing. However, the packages currently available for the analysis of raw sequencing data do not provide automatic analysis of complex experimental designs with multiple biological conditions and multiple analysis time-points.

RESULTS

The MultiRNAflow suite combines several packages in a unified framework allowing exploratory and supervised statistical analyses of temporal data for multiple biological conditions.

AVAILABILITY AND IMPLEMENTATION

The R package MultiRNAflow is freely available on Bioconductor (https://bioconductor.org/packages/MultiRNAflow/), and the latest version of the source code is available on a GitHub repository (https://github.com/loubator/MultiRNAflow).

Like mother like son: Transgenerational memory and cross-tolerance from drought to heat stress are identified in chloroplast proteome and seed provisioning in Pinus radiata

How different stressors impact plant health and memory when they are imposed in different generations in wild ecosystems is still scarce. Here, we address how different environments shape heritable memory for the next generation in seeds and seedlings of Pinus radiata, a long-lived species with economic interest. The performance of the seedlings belonging to two wild clonal subpopulations (optimal fertirrigation vs. slightly stressful conditions) was tested under heat stress through physiological profiling and comparative time-series chloroplast proteomics. In addition, we explored the seeds conducting a physiological characterization and targeted transcriptomic profiling in both subpopulations. Our results showed differential responses between them, evidencing a cross-stress transgenerational memory. Seedlings belonging to the stressed subpopulation retained key proteins related to Photosystem II, chloroplast-to-nucleus signalling and osmoprotection which helped to overcome the applied heat stress. The seeds also showed a differential gene expression profile for targeted genes and microRNAs, as well as an increased content of starch and secondary metabolites, molecules which showed potential interest as biomarkers for early selection of primed plants. Thus, these finds not only delve into transgenerational cross-stress memory in trees, but also provide a new biotechnological tool for forest design.

Mapping the global interactome of the ARF family reveals spatial organization in cellular signaling pathways

The ADP-ribosylation factors (ARFs) and ARF-like (ARL) GTPases serve as essential molecular switches governing a wide array of cellular processes. In this study, we used proximity-dependent biotin identification (BioID) to comprehensively map the interactome of 28 out of 29 ARF and ARL proteins in two cellular models. Through this approach, we identified ∼3000 high-confidence proximal interactors, enabling us to assign subcellular localizations to the family members. Notably, we uncovered previously undefined localizations for ARL4D and ARL10. Clustering analyses further exposed the distinctiveness of the interactors identified with these two GTPases. We also reveal that the expression of the understudied member ARL14 is confined to the stomach and intestines. We identified phospholipase D1 (PLD1) and the ESCPE-1 complex, more precisely, SNX1, as proximity interactors. Functional assays demonstrated that ARL14 can activate PLD1 in cellulo and is involved in cargo trafficking via the ESCPE-1 complex. Overall, the BioID data generated in this study provide a valuable resource for dissecting the complexities of ARF and ARL spatial organization and signaling.

OmicsFootPrint: a framework to integrate and interpret multi-omics data using circular images and deep neural networks

The OmicsFootPrint framework addresses the need for advanced multi-omics data analysis methodologies by transforming data into intuitive two-dimensional circular images and facilitating the interpretation of complex diseases. Utilizing Deep Neural Networks and incorporating the SHapley Additive exPlanations (SHAP) algorithm, the framework enhances model interpretability. Tested with The Cancer Genome Atlas (TCGA) data, OmicsFootPrint effectively classified lung and breast cancer subtypes, achieving high Area Under Curve (AUC) scores - 0.98±0.02 for lung cancer subtype differentiation, 0.83±0.07 for breast cancer PAM50 subtypes, and successfully distinguishe between invasive lobular and ductal carcinomas in breast cancer, showcasing its robustness. It also demonstrated notable performance in predicting drug responses in cancer cell lines, with a median AUC of 0.74, surpassing existing algorithms. Furthermore, its effectiveness persists even with reduced training sample sizes. OmicsFootPrint marks an enhancement in multi-omics research, offering a novel, efficient, and interpretable approach that contributes to a deeper understanding of disease mechanisms.

RNA-Seq Data Analysis: A Practical Guide for Model and Non-Model Organisms

RNA sequencing (RNA-seq) has emerged as a powerful tool for assessing genome-wide gene expression, revolutionizing various fields of biology. However, analyzing large RNA-seq datasets can be challenging, especially for students or researchers lacking bioinformatics experience. To address these challenges, we present a comprehensive guide to provide step-by-step workflows for analyzing RNA-seq data, from raw reads to functional enrichment analysis, starting with considerations for experimental design. This is designed to aid students and researchers working with any organism, irrespective of whether an assembled genome is available. Within this guide, we employ various recognized bioinformatics tools to navigate the landscape of RNA-seq analysis and discuss the advantages and disadvantages of different tools for the same task. Our protocol focuses on clarity, reproducibility, and practicality to enable users to navigate the complexities of RNA-seq data analysis easily and gain valuable biological insights from the datasets. Additionally, all scripts and a sample dataset are available in a GitHub repository to facilitate the implementation of the analysis pipeline. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Analysis of data from a model plant with an available reference genome Basic Protocol 2: Gene ontology enrichment analysis Basic Protocol 3: De novo assembly of data from non-model plants.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

The post-septic peripheral myeloid compartment reveals unexpected diversity in myeloid-derived suppressor cells

Introduction

Sepsis engenders distinct host immunologic changes that include the expansion of myeloid-derived suppressor cells (MDSCs). These cells play a physiologic role in tempering acute inflammatory responses but can persist in patients who develop chronic critical illness.

Methods

Cellular Indexing of Transcriptomes and Epitopes by Sequencing and transcriptomic analysis are used to describe MDSC subpopulations based on differential gene expression, RNA velocities, and biologic process clustering.

Results

We identify a unique lineage and differentiation pathway for MDSCs after sepsis and describe a novel MDSC subpopulation. Additionally, we report that the heterogeneous response of the myeloid compartment of blood to sepsis is dependent on clinical outcome.

Discussion

The origins and lineage of these MDSC subpopulations were previously assumed to be discrete and unidirectional; however, these cells exhibit a dynamic phenotype with considerable plasticity.

The microbiota and T cells non-genetically modulate inherited phenotypes transgenerationally

The host-microbiota relationship has evolved to shape mammalian physiology, including immunity, metabolism, and development. Germ-free models are widely used to study microbial effects on host processes such as immunity. Here, we find that both germ-free and T cell-deficient mice exhibit a robust sebum secretion defect persisting across multiple generations despite microbial colonization and T cell repletion. These phenotypes are inherited by progeny conceived during in vitro fertilization using germ-free sperm and eggs, demonstrating that non-genetic information in the gametes is required for microbial-dependent phenotypic transmission. Accordingly, gene expression in early embryos derived from gametes from germ-free or T cell-deficient mice is strikingly and similarly altered. Our findings demonstrate that microbial- and immune-dependent regulation of non-genetic information in the gametes can transmit inherited phenotypes transgenerationally in mice. This mechanism could rapidly generate phenotypic diversity to enhance host adaptation to environmental perturbations.

Long Non-Coding RNA Levels Are Modulated in Schistosoma mansoni following In Vivo Praziquantel Exposure

Schistosomiasis is a disease caused by trematodes of the genus Schistosoma that affects over 200 million people worldwide. For decades, praziquantel (PZQ) has been the only available drug to treat the disease. Despite recent discoveries that identified a transient receptor ion channel as the target of PZQ, schistosome response to this drug remains incompletely understood, since effectiveness relies on other factors that may trigger a complex regulation of parasite gene expression. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential that play important roles in S. mansoni homeostasis, reproduction, and fertility. Here, we show that in vivo PZQ treatment modulates lncRNA levels in S. mansoni. We re-analyzed public RNA-Seq data from mature and immature S. mansoni worms treated in vivo with PZQ and detected hundreds of lncRNAs differentially expressed following drug exposure, many of which are shared among mature and immature worms. Through RT-qPCR, seven out of ten selected lncRNAs were validated as differentially expressed; interestingly, we show that these lncRNAs are not adult worm stage-specific and are co-expressed with PZQ-modulated protein-coding genes. By demonstrating that parasite lncRNA expression levels alter in response to PZQ, this study unravels an important step toward elucidating the complex mechanisms of S. mansoni response to PZQ.

Coprophagia in early life tunes expression of immune genes after weaning in rabbit ileum

Coprophagia by suckling rabbits, i.e. ingestion of feces from their mother, reduces mortality after weaning. We hypothesized that this beneficial effect of coprophagia is immune-mediated at the intestinal level. Therefore, this study investigated immune development after weaning by analyzing the ileal transcriptome at day 35 and 49 in rabbits with differential access to coprophagia in early life. Rabbit pups had access between day 1 and 15 to (i) no feces (NF) or (ii) feces from unrelated does (Foreign Feces, FF) or (iii) feces from unrelated does treated with antibiotics (FFab). 350 genes were differentially expressed between day 35 and day 49 in suckling rabbits with access to coprophagia. These genes coded for antimicrobial peptides, a mucin, cytokines and chemokines, pattern recognition receptors, proteins involved in immunoglobulin A secretion and in interferon signaling pathway. Strikingly, prevention of coprophagia or access to feces from antibiotic-treated does in early life blunted immune development between day 35 et 49 in the ileum of rabbits. Thus, coprophagia might be crucial for the maturation of intestinal immunity in rabbits and could explain why this behavior improves survival.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

NAC61 regulates late- and post-ripening osmotic, oxidative, and biotic stress responses in grapevine

During late- and post-ripening stages, grape berry undergoes profound biochemical and physiological changes whose molecular control is poorly understood. Here, we report the role of NAC61, a grapevine NAC transcription factor, in regulating different processes involved in berry ripening progression. NAC61 is highly expressed during post-harvest berry dehydration and its expression pattern is closely related to sugar concentration. The ectopic expression of NAC61 in Nicotiana benthamiana leaves resulted in low stomatal conductance, high leaf temperature, tissue collapse and a higher relative water content. Transcriptome analysis of grapevine leaves transiently overexpressing NAC61 and DNA affinity purification and sequencing analyses allowed us to narrow down a list of NAC61-regulated genes. Direct regulation of the stilbene synthase regulator MYB14, the osmotic stress-related gene DHN1b, the Botrytis cinerea susceptibility gene WRKY52, and NAC61 itself was validated. We also demonstrate that NAC61 interacts with NAC60, a proposed master regulator of grapevine organ maturation, in the activation of MYB14 and NAC61 expression. Overall, our findings establish NAC61 as a key player in a regulatory network that governs stilbenoid metabolism and osmotic, oxidative, and biotic stress responses that are the hallmark of late- and post-ripening grape stages.

Host Immune Response Modulation in Avian Coronavirus Infection: Tracheal Transcriptome Profiling In Vitro and In Vivo

Infectious bronchitis virus (IBV) is a highly contagious Gammacoronavirus causing moderate to severe respiratory infection in chickens. Understanding the initial antiviral response in the respiratory mucosa is crucial for controlling viral spread. We aimed to characterize the impact of IBV Delmarva (DMV)/1639 and IBV Massachusetts (Mass) 41 at the primary site of infection, namely, in chicken tracheal epithelial cells (cTECs) in vitro and the trachea in vivo. We hypothesized that some elements of the induced antiviral responses are distinct in both infection models. We inoculated cTECs and infected young specific pathogen-free (SPF) chickens with IBV DMV/1639 or IBV Mass41, along with mock-inoculated controls, and studied the transcriptome using RNA-sequencing (RNA-seq) at 3 and 18 h post-infection (hpi) for cTECs and at 4 and 11 days post-infection (dpi) in the trachea. We showed that IBV DMV/1639 and IBV Mass41 replicate in cTECs in vitro and the trachea in vivo, inducing host mRNA expression profiles that are strain- and time-dependent. We demonstrated the different gene expression patterns between in vitro and in vivo tracheal IBV infection. Ultimately, characterizing host-pathogen interactions with various IBV strains reveals potential mechanisms for inducing and modulating the immune response during IBV infection in the chicken trachea.

Genome-wide identification and functional profile analysis of long non-coding RNAs in Avicennia marina

Avicennia marina, known for its remarkable adaptability to the challenging coastal environment, including high salinity, tide, and anaerobic soils, holds pivotal functions in safeguarding the coastal ecosystem. Long non-coding RNAs (lncRNAs) have emerged as significant players in various natural processes of plants such as development. However, lncRNAs in A. marina remain largely unknown and uncharacterized. Here, we employed the transcriptome datasets from multiple tissues, such as root, leaf, and seed, to detect and characterize the lncRNAs of A. marina. Analyzing synthetically, we finally identified 6333 lncRNAs in the A. marina. These lncRNAs exhibited distinct features compared to messenger RNAs, including larger exons, lower guanine-cytosine contents, lower expression levels, and higher tissue specificities. Moreover, we identified thousands of tissue-specific lncRNAs across the examined tissues and further found that these tissue-specific lncRNAs were significantly enriched in biological processes related to the major functions of their corresponding tissues. For instance, leaf-specific lncRNAs showed prominent enrichment in photosynthesis, oxidation-reduction processes, and light harvesting. By providing a comprehensive dataset and functional annotations for A. marina lncRNAs, this study offers a valuable overview of lncRNAs in A. marina and lays the fundamental foundation for further functional exploring of them.

TMBcalc: a computational pipeline for identifying pan-cancer Tumor Mutational Burden gene signatures

Background

In the precision medicine era, identifying predictive factors to select patients most likely to benefit from treatment with immunological agents is a crucial and open challenge in oncology.

Methods

This paper presents a pan-cancer analysis of Tumor Mutational Burden (TMB). We developed a novel computational pipeline, TMBcalc, to calculate the TMB. Our methodology can identify small and reliable gene signatures to estimate TMB from custom targeted-sequencing panels. For this purpose, our pipeline has been trained on top of 17 cancer types data obtained from TCGA.

Results

Our results show that TMB, computed through the identified signature, strongly correlates with TMB obtained from whole-exome sequencing (WES).

Conclusion

We have rigorously analyzed the effectiveness of our methodology on top of several independent datasets. In particular we conducted a comprehensive testing on: (i) 126 samples sourced from the TCGA database; few independent whole-exome sequencing (WES) datasets linked to colon, breast, and liver cancers, all acquired from the EGA and the ICGC Data Portal. This rigorous evaluation clearly highlights the robustness and practicality of our approach, positioning it as a promising avenue for driving substantial progress within the realm of clinical practice.

Nucleosome reorganisation in breast cancer tissues

BACKGROUND

Nucleosome repositioning in cancer is believed to cause many changes in genome organisation and gene expression. Understanding these changes is important to elucidate fundamental aspects of cancer. It is also important for medical diagnostics based on cell-free DNA (cfDNA), which originates from genomic DNA regions protected from digestion by nucleosomes.

RESULTS

We have generated high-resolution nucleosome maps in paired tumour and normal tissues from the same breast cancer patients using MNase-assisted histone H3 ChIP-seq and compared them with the corresponding cfDNA from blood plasma. This analysis has detected single-nucleosome repositioning at key regulatory regions in a patient-specific manner and common cancer-specific patterns across patients. The nucleosomes gained in tumour versus normal tissue were particularly informative of cancer pathways, with ~ 20-fold enrichment at CpG islands, a large fraction of which marked promoters of genes encoding DNA-binding proteins. The tumour tissues were characterised by a 5-10 bp decrease in the average distance between nucleosomes (nucleosome repeat length, NRL), which is qualitatively similar to the differences between pluripotent and differentiated cells. This effect was correlated with gene activity, differential DNA methylation and changes in local occupancy of linker histone variants H1.4 and H1X.

CONCLUSIONS

Our study offers a novel resource of high-resolution nucleosome maps in breast cancer patients and reports for the first time the effect of systematic decrease of NRL in paired tumour versus normal breast tissues from the same patient. Our findings provide a new mechanistic understanding of nucleosome repositioning in tumour tissues that can be valuable for patient diagnostics, stratification and monitoring.

Dissecting the spatiotemporal diversity of adult neural stem cells

Adult stem cells are important for tissue turnover and regeneration. However, in most adult systems it remains elusive how stem cells assume different functional states and support spatially patterned tissue architecture. Here, we dissected the diversity of neural stem cells in the adult zebrafish brain, an organ that is characterized by pronounced zonation and high regenerative capacity. We combined single-cell transcriptomics of dissected brain regions with massively parallel lineage tracing and in vivo RNA metabolic labeling to analyze the regulation of neural stem cells in space and time. We detected a large diversity of neural stem cells, with some subtypes being restricted to a single brain region, while others were found globally across the brain. Global stem cell states are linked to neurogenic differentiation, with different states being involved in proliferative and non-proliferative differentiation. Our work reveals principles of adult stem cell organization and establishes a resource for the functional manipulation of neural stem cell subtypes.

Mitochondrial complex I activity in microglia sustains neuroinflammation

Sustained smouldering, or low-grade activation, of myeloid cells is a common hallmark of several chronic neurological diseases, including multiple sclerosis1. Distinct metabolic and mitochondrial features guide the activation and the diverse functional states of myeloid cells2. However, how these metabolic features act to perpetuate inflammation of the central nervous system is unclear. Here, using a multiomics approach, we identify a molecular signature that sustains the activation of microglia through mitochondrial complex I activity driving reverse electron transport and the production of reactive oxygen species. Mechanistically, blocking complex I in pro-inflammatory microglia protects the central nervous system against neurotoxic damage and improves functional outcomes in an animal disease model in vivo. Complex I activity in microglia is a potential therapeutic target to foster neuroprotection in chronic inflammatory disorders of the central nervous system3.

Differential expression and alternative splicing analyses of multiple tissues reveal albinism-associated genes in the Wels catfish (Silurus glanis)

Albinism is a widespread departure from a typical body colouration due to altered melanin production. The Wels catfish (Silurus glanis) is among the largest freshwater fish species in the world, and albino individuals occur both in the wild and in aquaculture. Here, we performed transcriptome-wide analysis of albino and normally pigmented S. glanis using four tissues (skin, dorsal fin, whole eye and liver) to identify genes associated with albinism by exploring patterns of differential expression (DE) and differential alternative splicing (DAS). Multi-tissue analyses revealed a large number of genes in skin (n = 1355) and fin (n = 614) tissue associated with the albino phenotype in S. glanis, while the number of DE genes in eye and liver tissues was lower (n = 188, n = 189, respectively). Several DE genes across multiple tissues were detected as the most promising candidates (e.g., hsp4, hsp90b1, raph1, uqcrfs1, adcy-family and wnt-family) potentially causally linked to the albino phenotype in Wels catfish. Moreover, our findings supported earlier observations of physiological differences between albino and normally pigmented individuals, particularly in energy metabolism and immune response. In contrast, there were only a few pigmentation-related genes observed among DAS genes (4 in skin, 2 in fin), the overlap between DAS and DE genes was low (n = 25) and did not include known pigmentation-related genes. This suggests that DAS and DE in Wels catfish are, to a large extent, independent processes, and the observed alternative splicing cases are probably not causally linked with albinism in S. glanis. This work provides the first transcriptome-wide multi-tissue insights into the albinism of Wels catfish and serves as a valuable resource for further understanding the genetic mechanisms of pigmentation in fish.

Silencing Apoe with divalent-siRNAs improves amyloid burden and activates immune response pathways in Alzheimer's disease

INTRODUCTION

The most significant genetic risk factor for late-onset Alzheimer's disease (AD) is APOE4, with evidence for gain- and loss-of-function mechanisms. A clinical need remains for therapeutically relevant tools that potently modulate APOE expression.

METHODS

We optimized small interfering RNAs (di-siRNA, GalNAc) to potently silence brain or liver Apoe and evaluated the impact of each pool of Apoe on pathology.

RESULTS

In adult 5xFAD mice, siRNAs targeting CNS Apoe efficiently silenced Apoe expression and reduced amyloid burden without affecting systemic cholesterol, confirming that potent silencing of brain Apoe is sufficient to slow disease progression. Mechanistically, silencing Apoe reduced APOE-rich amyloid cores and activated immune system responses.

DISCUSSION

These results establish siRNA-based modulation of Apoe as a viable therapeutic approach, highlight immune activation as a key pathway affected by Apoe modulation, and provide the technology to further evaluate the impact of APOE silencing on neurodegeneration.

Systemic proteome adaptions to 7-day complete caloric restriction in humans

Surviving long periods without food has shaped human evolution. In ancient and modern societies, prolonged fasting was/is practiced by billions of people globally for religious purposes, used to treat diseases such as epilepsy, and recently gained popularity as weight loss intervention, but we still have a very limited understanding of the systemic adaptions in humans to extreme caloric restriction of different durations. Here we show that a 7-day water-only fast leads to an average weight loss of 5.7 kg (±0.8 kg) among 12 volunteers (5 women, 7 men). We demonstrate nine distinct proteomic response profiles, with systemic changes evident only after 3 days of complete calorie restriction based on in-depth characterization of the temporal trajectories of ~3,000 plasma proteins measured before, daily during, and after fasting. The multi-organ response to complete caloric restriction shows distinct effects of fasting duration and weight loss and is remarkably conserved across volunteers with >1,000 significantly responding proteins. The fasting signature is strongly enriched for extracellular matrix proteins from various body sites, demonstrating profound non-metabolic adaptions, including extreme changes in the brain-specific extracellular matrix protein tenascin-R. Using proteogenomic approaches, we estimate the health consequences for 212 proteins that change during fasting across ~500 outcomes and identified putative beneficial (SWAP70 and rheumatoid arthritis or HYOU1 and heart disease), as well as adverse effects. Our results advance our understanding of prolonged fasting in humans beyond a merely energy-centric adaptions towards a systemic response that can inform targeted therapeutic modulation.

Targeting circulating labile heme as a defense strategy against malaria

Severe presentations of malaria emerge as Plasmodium (P.) spp. parasites invade and lyse red blood cells (RBC), producing extracellular hemoglobin (HB), from which labile heme is released. Here, we tested whether scavenging of extracellular HB and/or labile heme, by haptoglobin (HP) and/or hemopexin (HPX), respectively, counter the pathogenesis of severe presentations of malaria. We found that circulating labile heme is an independent risk factor for cerebral and non-cerebral presentations of severe P. falciparum malaria in children. Labile heme was negatively correlated with circulating HP and HPX, which were, however, not risk factors for severe P. falciparum malaria. Genetic Hp and/or Hpx deletion in mice led to labile heme accumulation in plasma and kidneys, upon Plasmodium infection This was associated with higher incidence of mortality and acute kidney injury (AKI) in ageing but not adult Plasmodium-infected mice, and was corroborated by an inverse correlation between heme and HPX with serological markers of AKI in P. falciparum malaria. In conclusion, HP and HPX act in an age-dependent manner to prevent the pathogenesis of severe presentation of malaria in mice and presumably in humans.

Single-cell multiplex chromatin and RNA interactions in ageing human brain

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an 'older' transcriptomic signature and Alzheimer's disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

The zebrafish heart harbors a thermogenic beige fat depot analog of human epicardial adipose tissue

Epicardial adipose tissue (eAT) is a metabolically active fat depot that has been associated with a wide array of cardiac homeostatic functions and cardiometabolic diseases. A full understanding of its diverse physiological and pathological roles is hindered by the dearth of animal models. Here, we show, in the heart of an ectothermic teleost, the zebrafish, the existence of a fat depot localized underneath the epicardium, originating from the epicardium and exhibiting the molecular signature of beige adipocytes. Moreover, a subset of adipocytes within this cardiac fat tissue exhibits primitive thermogenic potential. Transcriptomic profiling and cross-species analysis revealed elevated glycolytic and cardiac homeostatic gene expression with downregulated obesity and inflammatory hallmarks in the teleost eAT compared to that of lean aged humans. Our findings unveil epicardium-derived beige fat in the heart of an ectotherm considered to possess solely white adipocytes for energy storage and identify pathways that may underlie age-driven remodeling of human eAT.

A maternal-effect Padi6 variant causes nuclear and cytoplasmic abnormalities in oocytes, as well as failure of epigenetic reprogramming and zygotic genome activation in embryos

Maternal inactivation of genes encoding components of the subcortical maternal complex (SCMC) and its associated member, PADI6, generally results in early embryo lethality. In humans, SCMC gene variants were found in the healthy mothers of children affected by multilocus imprinting disturbances (MLID). However, how the SCMC controls the DNA methylation required to regulate imprinting remains poorly defined. We generated a mouse line carrying a Padi6 missense variant that was identified in a family with Beckwith-Wiedemann syndrome and MLID. If homozygous in female mice, this variant resulted in interruption of embryo development at the two-cell stage. Single-cell multiomic analyses demonstrated defective maturation of Padi6 mutant oocytes and incomplete DNA demethylation, down-regulation of zygotic genome activation (ZGA) genes, up-regulation of maternal decay genes, and developmental delay in two-cell embryos developing from Padi6 mutant oocytes but little effect on genomic imprinting. Western blotting and immunofluorescence analyses showed reduced levels of UHRF1 in oocytes and abnormal localization of DNMT1 and UHRF1 in both oocytes and zygotes. Treatment with 5-azacytidine reverted DNA hypermethylation but did not rescue the developmental arrest of mutant embryos. Taken together, this study demonstrates that PADI6 controls both nuclear and cytoplasmic oocyte processes that are necessary for preimplantation epigenetic reprogramming and ZGA.

Protocol for identifying differentially expressed genes using the RumBall RNA-seq analysis platform

Here, we present a protocol for the identification of differentially expressed genes through RNA sequencing analysis. Starting with FASTQ files from public datasets, this protocol leverages RumBall within a self-contained Docker system. We describe the steps for software setup, obtaining data, read mapping, sample normalization, statistical modeling, and gene ontology enrichment. We then detail procedures for interpreting results with plots and tables. RumBall internally utilizes popular tools, ensuring a comprehensive understanding of the analysis process.

Chronological transitions of hepatocyte growth factor treatment effects in spinal cord injury tissue

Inflammatory responses are known to suppress neural regeneration in patients receiving stem cell-based regenerative therapy for spinal cord injury (SCI). Consequently, pathways involved in neurogenesis and immunomodulation, such as the hepatocyte growth factor (HGF)/MET signaling cascade, have garnered significant attention. Notably, various studies, including our own, have highlighted the enhanced recovery of locomotor functions achieved in SCI animal models by combining HGF pretreatment and human induced stem cell-derived neural stem/progenitor cell (hiPSC-NS/PC) transplantation. However, these studies implicitly hypothesized that the functionality of HGF in SCI would be time consistent and did not elucidate its dynamics. In the present article, we investigated the time-course of the effect of HGF on SCI, aiming to uncover a more precise mechanism for HGF administration, which is indispensable for developing crystallizing protocols for combination therapy. To this end, we performed a detailed investigation of the temporal variation of HGF using the RNA-seq data we obtained in our most recent study. Leveraging the time-series design of the data, which we did not fully exploit previously, we identified three components in the effects of HGF that operate at different times: early effects, continuous effects, and delayed effects. Our findings suggested a concept where the three components together contribute to the acceleration of neurogenesis and immunomodulation, which reinforce the legitimacy of empirically fine-tuned protocols for HGF administration and advocate the novel possibility that the time-inconsistent effects of HGF progressively augment the efficacy of combined therapy.