Pancreatic neuroendocrine tumors (PanNETs), the second most common type of primary pancreatic tumors, display notable heterogeneity in invasiveness. Current knowledge regarding genomic alterations, including DAXX/ATRX, MEN1 mutations, and copy number variations (CNVs), provides some insights into tumor invasiveness. However, the underlying reasons for the significant variation in invasiveness between insulinoma and other types of PanNETs remain unclear. To construct a comprehensive model for the stratification of prognosis, we employed analysis of both the well-established Rip1-Tag2 (RT2) mouse model of PanNETs and human PanNETs with various functional types. Firstly, by applying single-cell and bulk RNA sequencing in PanNETs from different ages and strains of RT2 mice and human PanNETs, we introduced a 2-dimensional (2D) classification system. Based on the 2-D classification system, human PanNETs were mainly classified as benign insulinomas or non-insulinomas subclusters. Non-insulinomas subtypes mainly included gastrinomas, glucagonomas, VIPomas, and NF-PanNETs, which all exhibited potential invasiveness. In addition, we discovered an enrichment of specific CNV patterns and mutations in corresponding human PanNET subclusters. Then we denoted somatic DAXX/ATRX as the \'second hit\' and confounding factors for invasiveness. Finally, by combining the 2D system, DAXX/ATRX mutation status, and tumor diameter, a group of indolent PanNETs with minimal recurrence risk was identified. In conclusion, our current work constructed a comprehensive model to elucidate the heterogeneity of invasiveness in PanNETs and improve prognostic stratification.

Depression, or major depressive disorder, poses a significant burden for both individuals and society, affecting approximately 10.8% of the general population. This psychiatric disorder leads to approximately 800,000 deaths per year. A combination of genetic and environmental factors such as early life stress (ELS) increase the risk for development of depression in humans, and a clear role for the hippocampus in the pathophysiology of depression has been shown. Nevertheless, the underlying mechanisms of depression remain poorly understood, resulting in a lack of effective treatments. To better understand the core mechanisms underlying the development of depression, we used a cross-species design to investigate shared hippocampal pathophysiological mechanisms in mouse ELS and human depression. Mice were subjected to ELS by a maternal separation paradigm, followed by RNA sequencing analysis of the adult hippocampal tissue. This identified persistent transcriptional changes linked to mitochondrial stress response pathways, with oxidative phosphorylation and protein folding emerging as the main mechanisms affected by maternal separation. Remarkably, there was a significant overlap between the pathways involved in mitochondrial stress response we observed and publicly available RNAseq data from hippocampal tissue of depressive patients. This cross-species conservation of changes in gene expression of mitochondria-related genes suggests that mitochondrial stress may play a pivotal role in the development of depression. Our findings highlight the potential significance of the hippocampal mitochondrial stress response as a core mechanism underlying the development of depression. Further experimental investigations are required to expand our understanding of these mechanisms.

Cardiovascular research heavily relies on mouse (Mus musculus) models to study disease mechanisms and to test novel biomarkers and medications. Yet, applying these results to patients remains a major challenge and often results in noneffective drugs. Therefore, it is an open challenge of translational science to develop models with high similarities and predictive value. This requires a comparison of disease models in mice with diseased tissue derived from humans.
To compare the transcriptional signatures at single-cell resolution, we implemented an integration pipeline called OrthoIntegrate, which uniquely assigns orthologs and therewith merges single-cell RNA sequencing (scRNA-seq) RNA of different species. The pipeline has been designed to be as easy to use and is fully integrable in the standard Seurat workflow.We applied OrthoIntegrate on scRNA-seq from cardiac tissue of heart failure patients with reduced ejection fraction (HFrEF) and scRNA-seq from the mice after chronic infarction, which is a commonly used mouse model to mimic HFrEF. We discovered shared and distinct regulatory pathways between human HFrEF patients and the corresponding mouse model. Overall, 54% of genes were commonly regulated, including major changes in cardiomyocyte energy metabolism. However, several regulatory pathways (e.g., angiogenesis) were specifically regulated in humans.
The demonstration of unique pathways occurring in humans indicates limitations on the comparability between mice models and human HFrEF and shows that results from the mice model should be validated carefully. OrthoIntegrate is publicly accessible (https://github.com/MarianoRuzJurado/OrthoIntegrate) and can be used to integrate other large datasets to provide a general comparison of models with patient data.

The inner ear of mammals includes the cochlea and vestibule, which house specialized hair cells that are responsible for hearing and balance, respectively. While cochlear hair cells fail to regenerate following damage, those of the utricle, which is part of the vestibular apparatus, show partial regeneration. In birds, the macula lagena, a unique ear structure in this clade, has the ability to regenerate hair cells similarly to the utricle. Many studies have sought to explain regeneration in terms of evolution and species differences. However, it remains unclear what the cellular and molecular basis is behind the differences in inner ear structures and between avians and mammals. In the present study, we first investigated the anatomical structures of the inner ear of both chickens and rodents. We then performed RNA sequencing (RNA-Seq) and made cross-species analyses of the expression of homologous genes obtained from the inner ear tissue from both chickens and mice. Finally, we focused on the lagena, the basilar papilla, and the utricle in chickens and identified differentially expressed genes between tissues and determined the expression patterns of genes involved in inner ear structure formation by single-cell RNA sequencing and bulk RNA-Seq. We concluded that the cellular and molecular composition of the lagena is more similar to that of the utricle than the cochlea. Taken together, our study provides a valuable resource for the study of inner ear evolution and development.

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

Horticultural plants contribute immensely to the quality of human\'s life. The rapid development of omics studies on horticultural plants has resulted in large volumes of valuable growth- and development-related data. Genes that are essential for growth and development are highly conserved in evolution. Cross-species data mining reduces the impact of species heterogeneity and has been extensively used for conserved gene identification. Owing to the lack of a comprehensive database for cross-species data mining using multi-omics data from all horticultural plant species, the current resources in this field are far from satisfactory. Here, we introduce GERDH (https://dphdatabase.com), a database platform for cross-species data mining among horticultural plants, based on 12 961 uniformly processed publicly available omics libraries from more than 150 horticultural plant accessions, including fruits, vegetables and ornamental plants. Important and conserved genes that are essential for a specific biological process can be obtained by cross-species analysis module with interactive web-based data analysis and visualization. Moreover, GERDH is equipped with seven online analysis tools, including gene expression, in-species analysis, epigenetic regulation, gene co-expression, enrichment/pathway and phylogenetic analysis. By interactive cross-species analysis, we identified key genes contributing to postharvest storage. By gene expression analysis, we explored new functions of CmEIN3 in flower development, which was validated by transgenic chrysanthemum analysis. We believe that GERDH will be a useful resource for key gene identification and will allow for omics big data to be more available and accessible to horticultural plant community members.

Cetaceans (dolphins, whales, and porpoises) have large and anatomically sophisticated brains. To expand our understanding of the cellular makeup of cetacean brains and the similarities and divergence between the brains of cetaceans and terrestrial mammals, we report a short-finned pilot whale (Globicephala macrorhynchus) single-nucleus transcriptome atlas. To achieve this goal, we assembled a chromosome-scale reference genome spanning 2.25 Gb on 22 chromosomes and profiled the gene expression of five major anatomical cortical regions of the short-finned pilot whale by single-nucleus RNA-sequencing (snRNA-seq). We identified six major cell lineages in the cerebral cortex (excitatory neurons, inhibitory neurons, oligodendrocytes, oligodendrocyte precursor cells, astrocytes, and endothelial cells), eight molecularly distinct subclusters of excitatory neurons, and four subclusters of inhibitory neurons. Finally, a comparison of snRNA-seq data from the short-finned pilot whale, human, and rhesus macaque revealed a broadly conserved cellular makeup of brain cell types. Our study provides genomic resources and molecular insights into cetacean brain evolution.

Although the imbalance of circulating levels of Thyroid Hormones (THs) affects female fertility in vertebrates, its involvement in the promotion of Premature Ovarian Aging (POA) is debated. Therefore, altered synthesis of THs in both thyroid and ovary can be a trait of POA. We investigated the relationship between abnormal TH signaling, dysthyroidism, and POA in evolutionary distant vertebrates: from zebrafish to humans. Ovarian T3 signaling/metabolism was evaluated by measuring T3 levels, T3 responsive transcript, and protein levels along with transcripts governing T3 availability (deiodinases) and signaling (TH receptors) in distinct models of POA depending on genetic background and environmental exposures (e.g., diets, pesticides). Expression levels of well-known (Amh, Gdf9, and Inhibins) and novel (miR143/145 and Gas5) biomarkers of POA were assessed. Ovarian dysthyroidism was slightly influenced by genetics since very few differences were found between C57BL/6J and FVB/NJ females. However, diets exacerbated it in a strain-dependent manner. Similar findings were observed in zebrafish and mouse models of POA induced by developmental and long-life exposure to low-dose chlorpyrifos (CPF). Lastly, the T3 decrease in follicular fluids from women affected by diminished ovarian reserve, as well as of the transcripts modulating T3 signaling/availability in the cumulus cells, confirmed ovarian dysthyroidism as a common and evolutionary conserved trait of POA.

UNASSIGNED: Mouse models are the basis for primary Sjögren\'s syndrome (pSS) research. However, the depth of comparisons between mice and humans in salivary gland (SG) immune cells remains limited.
UNASSIGNED: The gene expression profiles of SGs from normal subjects and pSS patients were downloaded from the Gene Expression Comprehensive Database. The proportion of infiltrating immune cell subsets was then assessed by cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT). An experimental Sjögren\'s syndrome (ESS) mouse model was successfully constructed using SG protein. Based on mouse SG tissue RNA-Seq data, the seq-ImmuCC model was used to quantitatively analyze the compositional ratios of 10 immune cells in pSS patients and mouse model SG tissues.
UNASSIGNED: Computed and obtained 31 human data samples using the CIBERSORT deconvolution method. The immune cell infiltration results showed that, compared to normal human SG tissue, the content of gamma delta T cells was significantly different from naive CD4+ T cells and significantly increased, while the plasma cell content decreased. Principal component analysis indicated differences in immune cell infiltration between pSS patients and normal subjects. Meanwhile, for ESS model mouse data analysis, we found that the proportion of macrophages increased, while the proportion of CD4+ T cells, B cells, and monocytes decreased. Furthermore, we found that the proportion of monocytes was decreased, while the proportion of macrophages was increased in the SG tissues of pSS patients and model mice. The infiltration of CD4+ T, CD8+ T, and B cells also showed some differences.
UNASSIGNED: We comprehensively analyzed SG immune infiltration in pSS patients and model mice. We demonstrated conserved and nonconserved aspects of the immune system in mice and humans at the level of immune cells to help explain the primary regulation of immune mechanisms during the development of Sjögren\'s syndrome.

Despite the global prevalence of nonalcoholic fatty liver disease (NAFLD), its pathophysiology remains unclear. In this study, we established highly confident nonalcoholic steatohepatitis (NASH) gene signatures and evaluated the pathological mechanisms underlying NASH through a systematic meta-analysis of transcriptome and proteome datasets obtained from NASH patients and mouse models.
We analyzed NASH transcriptome datasets from 539 patients and 99 mice. A whole-liver tissue proteome dataset was used to confirm the protein level dysregulation of NASH signatures significant in both humans and mice.
In total, 254 human and 1,917 mouse NASH gene signatures were established. Up-regulated genes of 254 human signatures were associated with inflammation, steatosis, apoptosis, and extracellular matrix organization, whereas down-regulated genes were associated with response to metal ions and lipid and amino acid metabolism. When different mouse models were compared against humans, models with high fat and high fructose diet most closely resembled the genetic features of human NAFLD. Cross-species analysis revealed 66 genes that were concordantly dysregulated between human and mouse NASH. Among these, 14 genes were further validated to be dysregulated at the protein level. The resulting 14 genes included some of the well-established NASH associated genes and a promising NASH drug target. Functional enrichment analysis revealed that dysregulation of amino acid metabolism was the most significant hepatic perturbation in both human and mouse NASH.
We established the most comprehensive hepatic gene signatures for NASH in humans and mice to date. To the best of our knowledge, this is the first study to collectively analyze the common signatures between human and mouse NASH on a transcriptome-proteome scale.