Single-cell atlas of domestic pig cerebral cortex and hypothalamus

Comparative transcriptome analysis of T lymphocyte subpopulations and identification of critical regulators defining porcine thymocyte identity

Introduction

The development and migration of T cells in the thymus and peripheral tissues are crucial for maintaining adaptive immunity in mammals. However, the regulatory mechanisms underlying T cell development and thymocyte identity formation in pigs remain largely underexplored.

Method

Here, by integrating bulk and single-cell RNA-sequencing data, we investigated regulatory signatures of porcine thymus and lymph node T cells.

Results

The comparison of T cell subpopulations derived from porcine thymus and lymph nodes revealed that their transcriptomic differences were influenced more by tissue origin than by T cell phenotypes, and that lymph node cells exhibited greater transcriptional diversity than thymocytes. Through weighted gene co-expression network analysis (WGCNA), we identified the key modules and candidate hub genes regulating the heterogeneity of T cell subpopulations. Further, we integrated the porcine thymocyte dataset with peripheral blood mononuclear cell (PBMC) dataset to systematically compare transcriptomic differences between T cell types from different tissues. Based on single-cell datasets, we further identified the key transcription factors (TFs) responsible for maintaining porcine thymocyte identity and unveiled that these TFs coordinately regulated the entire T cell development process. Finally, we performed GWAS of cell type-specific differentially expressed genes (DEGs) and 30 complex traits, and found that the DEGs in thymus-related and peripheral blood-related cell types, especially CD4_SP cluster and CD8-related cluster, were significantly associated with pig productive and reproductive traits.

Discussion

Our findings provide an insight into T cell development and lay a foundation for further exploring the porcine immune system and genetic mechanisms underlying complex traits in pigs.

Single-nucleus transcriptome sequencing reveals hepatic cell atlas in pigs

BACKGROUND

As the largest substantive organ of animals, the liver plays an essential role in the physiological processes of digestive metabolism and immune defense. However, the cellular composition of the pig liver remains poorly understood. This investigation used single-nucleus RNA sequencing technology to identify cell types from liver tissues of pigs, providing a theoretical basis for further investigating liver cell types in pigs.

RESULTS

The analysis revealed 13 cells clusters which were further identified 7 cell types including endothelial cells, T cells, hepatocytes, Kupffer cells, stellate cells, B cells, and cholangiocytes. The dominant cell types were endothelial cells, T cells and hepatocytes in the liver tissue of Dahe pigs and Dahe black pigs, which accounts for about 85.76% and 82.74%, respectively. The number of endothelial cells was higher in the liver tissue of Dahe pigs compared to Dahe black pigs, while the opposite tendency was observed for T cells. Moreover, functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic endothelial cells were significantly enriched in the protein processing in endoplasmic reticulum, MAPK signaling pathway, and FoxO signaling pathway. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic T cells were significantly enriched in the thyroid hormone signaling pathway, B cell receptor signaling pathway, and focal adhesion. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic hepatocytes were significantly enriched in the metabolic pathways.

CONCLUSIONS

In summary, this study provides a comprehensive cell atlas of porcine hepatic tissue. The number, gene expression level and functional characteristics of each cell type in pig liver tissue varied between breeds.

Single-cell dynamics of liver development in postnatal pigs

The postnatal development of the liver, an essential organ for metabolism and immunity, remains poorly characterized at the single-cell resolution. Here, we generated single-nucleus and single-cell transcriptomes of 84,824 pig liver cells at four postnatal time points: day 30, 42, 150, and 730. We uncovered 23 cell types, including three rare cell types: plasmacytoid dendritic cells, CAVIN3+IGF2+ endothelial cells, and EBF1+ fibroblasts. The latter two were verified by multiplex immunohistochemistry. Trajectory and gene regulatory analyses revealed 33 genes that encode transcription factors associated with hepatocyte development and function, including NFIL3 involved in regulating hepatic metabolism. We characterized the spatiotemporal heterogeneity of liver endothelial cells, identified and validated leucine zipper protein 2 (LUZP2) as a novel adult liver sinusoidal endothelial cell-specific transcription factor. Lymphoid cells (NK and T cells) governed the immune system of the pig liver since day 30. Furthermore, we identified a cluster of tissue-resident NK cells, which displayed virus defense functions, maintained proliferative features at day 730, and manifested a higher conservative transcription factor expression pattern in humans than in mouse liver. Our study presents the most comprehensive postnatal liver development single-cell atlas and demonstrates the metabolic and immune changes across the four age stages.

Comprehensive Genome and Transcriptome Analysis Identifies SLCO3A1 Associated with Aggressive Behavior in Pigs

Copy number variation (CNV) represents a significant reservoir of genetic diversity within the genome and exhibits a strong association with economically valuable traits in livestock. The manifestation of aggressive behavior in pigs has detrimental effects on production efficiency, immune competency, and meat quality. Nevertheless, the impact of CNV on the aggressive behavior of pigs remains elusive. In this investigation, we employed an integrated analysis of genome and transcriptome data to investigate the interplay between CNV, gene expression changes, and indicators of aggressive behavior in weaned pigs. Specifically, a subset of pigs comprising the most aggressive pigs (MAP, n = 12) and the least aggressive pigs (LAP, n = 11) was purposefully selected from a herd of 500 weaned pigs following a mixing procedure based on their composite aggressive score (CAS). Subsequently, we thoroughly analyzed copy number variation regions (CNVRs) across the entire genome using next-generation sequencing techniques, ultimately revealing the presence of 6869 CNVRs. Using genome-wide association study (GWAS) analysis and evaluating variance-stabilizing transformation (VST) values, we successfully identified distinct CNVRs that distinguished the MAP and LAP counterparts. Among the prioritized CNVRs, CNVR-4962 (designated as the top-ranked p-value and VST value, No. 1) was located within the Solute Carrier Organic Anion Transporter Family Member 3A1 (SLCO3A1) gene. The results of our analyses indicated a significantly higher (p < 0.05) copy number of SLCO3A1 in the MAP compared to the LAP. Furthermore, this increased copy number exhibited a positive correlation with the CAS of the pigs (p < 0.05). Furthermore, we integrated genomic data with transcriptomic data from the temporal lobe to facilitate the examination of expression quantitative trait loci (eQTL). Importantly, these observations were consistent with the mRNA expression pattern of SLCO3A1 in the temporal lobe of both MAP and LAP (p < 0.05). Consequently, our findings strongly suggest that CNVs affecting SLCO3A1 may influence gene expression through a dosage effect. These results highlight the potential of SLCO3A1 as a candidate gene associated with aggressive traits in pig breeding programs.

A genome and single-nucleus cerebral cortex transcriptome atlas of the short-finned pilot whale Globicephala macrorhynchus

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.

Single-cell omics: A new direction for functional genetic research in human diseases and animal models

Over the past decade, with the development of high-throughput single-cell sequencing technology, single-cell omics has been emerged as a powerful tool to understand the molecular basis of cellular mechanisms and refine our knowledge of diverse cell states. They can reveal the heterogeneity at different genetic layers and elucidate their associations by multiple omics analysis, providing a more comprehensive genetic map of biological regulatory networks. In the post-GWAS era, the molecular biological mechanisms influencing human diseases will be further elucidated by single-cell omics. This review mainly summarizes the development and trend of single-cell omics. This involves single-cell omics technologies, single-cell multi-omics technologies, multiple omics data integration methods, applications in various human organs and diseases, classic laboratory cell lines, and animal disease models. The review will reveal some perspectives for elucidating human diseases and constructing animal models.

An optimized FACS-free single-nucleus RNA sequencing (snRNA-seq) method for plant science research

Recently, single-cell RNA sequencing (scRNA-seq) provides unprecedented power for accurately understanding gene expression regulatory mechanisms. However, scRNA-seq studies have limitations in plants, due to difficulty in protoplast isolation that requires enzymatic digestion of the cell walls from various plant tissues. Therefore, to overcome this problem, we developed a nuclei isolation approach that does not rely on Fluorescence Activated Cell Sorting (FACS). We validated the robustness of the FACS-free single-nucleus RNA sequencing (snRNA-seq) methodology in mature Arabidopsis plant tissue by comparing it to scRNA-seq results based on protoplasts extracted from the same batch of leaf materials. Sequencing results demonstrated the high quality of snRNA-seq data, as well as its utility in cell type classification and marker gene identification. This approach also showed several advantages, including the ability to use frozen samples, taking less suspension preparation time, and reducing biased cellular coverage and dissociation-induced transcriptional artifacts. Surprisingly, snRNA-seq detected two epidermal pavement cell clusters, while scRNA-seq only had one. Furthermore, we hypothesized that these two epidermal cells represent the top and lower epidermis based on differences in expression patterns of cluster-specific expressed genes. In summary, this study has advanced the application of snRNA-seq in Arabidopsis leaves and confirmed the advantages of snRNA-seq in plant research.

Integrative cross-species analysis of GABAergic neuron cell types and their functions in Alzheimer's disease

Understanding the phenotypic and functional diversity of cerebral cortical GABAergic neurons requires a comprehensive analysis of key transcriptional signatures and neuronal subtype identity. However, the diversity and conservation of GABAergic neurons across multiple mammals remain unclear. Here, we collected the single-nucleus RNA sequencing (snRNA-seq) datasets of cerebral cortex from human, macaque, mouse, and pig to identify the conserved neuronal cell types across species. After systematic analysis of the heterogeneity of GABAergic neurons, we defined four major conserved GABAergic neuron subclasses (Inc SST, Inc LAMP5, Inc PVALB, and Inc VIP) across species. We characterized the species-enriched subclasses of GABAergic neurons from four mammals, such as Inc Meis2 in mouse. Then, we depicted the genetic regulatory network (GRNs) of GABAergic neuron subclasses, which showed the conserved and species-specific GRNs for GABAergic neuron cell types. Finally, we investigated the GABAergic neuron subclass-specific expression modules of Alzheimer's disease (AD)-related genes in GABAergic neuron cell types. Overall, our study reveals the conserved and divergent GABAergic neuron subclasses and GRNs across multiple species and unravels the gene expression modules of AD-risk genes in GABAergic neuron subclasses, facilitating the GABAergic neurons research and clinical treatment.

SPEED: Single-cell Pan-species atlas in the light of Ecology and Evolution for Development and Diseases

It is a challenge to efficiently integrate and present the tremendous amounts of single-cell data generated from multiple tissues of various species. Here, we create a new database named SPEED for single-cell pan-species atlas in the light of ecology and evolution for development and diseases (freely accessible at http://8.142.154.29 or http://speedatlas.net). SPEED is an online platform with 4 data modules, 7 function modules and 2 display modules. The 'Pan' module is applied for the interactive analysis of single cell sequencing datasets from 127 species, and the 'Evo', 'Devo', and 'Diz' modules provide comprehensive analysis of single-cell atlases on 18 evolution datasets, 28 development datasets, and 85 disease datasets. The 'C2C', 'G2G' and 'S2S' modules explore intercellular communications, genetic regulatory networks, and cross-species molecular evolution. The 'sSearch', 'sMarker', 'sUp', and 'sDown' modules allow users to retrieve specific data information, obtain common marker genes for cell types, freely upload, and download single-cell datasets, respectively. Two display modules ('HOME' and 'HELP') offer easier access to the SPEED database with informative statistics and detailed guidelines. All in all, SPEED is an integrated platform for single-cell RNA sequencing (scRNA-seq) and single-cell whole-genome sequencing (scWGS) datasets to assist the deep-mining and understanding of heterogeneity among cells, tissues, and species at multi-levels, angles, and orientations, as well as provide new insights into molecular mechanisms of biological development and pathogenesis.

Intestinal single-cell atlas reveals novel lymphocytes in pigs with similarities to human cells

Single-cell RNA sequencing of porcine ileal lymphocytes reveals similarities to human cells and discovery of porcine intestinal innate lymphoid cells.

Lymphocytes can heavily influence intestinal health, but resolving intestinal lymphocyte function is challenging as the intestine contains a vastly heterogeneous mixture of cells. Pigs are an advantageous biomedical model, but deeper understanding of intestinal lymphocytes is warranted to improve model utility. Twenty-six cell types were identified in the porcine ileum by single-cell RNA sequencing and further compared with cells in human and murine ileum. Though general consensus of cell subsets across species was revealed, some porcine-specific lymphocyte subsets were identified. Differential tissue dissection and in situ analyses conferred spatial context, revealing similar locations of lymphocyte subsets in Peyer’s patches and epithelium in pig-to-human comparisons. Like humans, activated and effector lymphocytes were abundant in the ileum but not periphery of pigs, suggesting tissue-specific and/or activation-associated gene expression. Gene signatures for peripheral and ileal innate lymphoid cells newly discovered in pigs were defined and highlighted similarities to human innate lymphoid cells. Overall, we reveal novel lymphocyte subsets in pigs and highlight utility of pigs for intestinal research applications.

Screening of cell-virus, cell-cell, gene-gene crosstalk among animal kingdom at single cell resolution

BACKGROUND

The exact animal origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains obscure and understanding its host range is vital for preventing interspecies transmission.

METHODS

Herein, we applied single-cell sequencing to multiple tissues of 20 species (30 data sets) and integrated them with public resources (45 data sets covering 26 species) to expand the virus receptor distribution investigation. While the binding affinity between virus and receptor is essential for viral infectivity, understanding the receptor distribution could predict the permissive organs and tissues when infection occurs.

RESULTS

Based on the transcriptomic data, the expression profiles of receptor or associated entry factors for viruses capable of causing respiratory, blood, and brain diseases were described in detail. Conserved cellular connectomes and regulomes were also identified, revealing fundamental cell-cell and gene-gene cross-talks from reptiles to humans.

CONCLUSIONS

Overall, our study provides a resource of the single-cell atlas of the animal kingdom which could help to identify the potential host range and tissue tropism of viruses and reveal the host-virus co-evolution.

Endothelial cell heterogeneity and microglia regulons revealed by a pig cell landscape at single-cell level

Pigs are valuable large animal models for biomedical and genetic research, but insights into the tissue- and cell-type-specific transcriptome and heterogeneity remain limited. By leveraging single-cell RNA sequencing, we generate a multiple-organ single-cell transcriptomic map containing over 200,000 pig cells from 20 tissues/organs. We comprehensively characterize the heterogeneity of cells in tissues and identify 234 cell clusters, representing 58 major cell types. In-depth integrative analysis of endothelial cells reveals a high degree of heterogeneity. We identify several functionally distinct endothelial cell phenotypes, including an endothelial to mesenchymal transition subtype in adipose tissues. Intercellular communication analysis predicts tissue- and cell type-specific crosstalk between endothelial cells and other cell types through the VEGF, PDGF, TGF-β, and BMP pathways. Regulon analysis of single-cell transcriptome of microglia in pig and 12 other species further identifies MEF2C as an evolutionally conserved regulon in the microglia. Our work describes the landscape of single-cell transcriptomes within diverse pig organs and identifies the heterogeneity of endothelial cells and evolutionally conserved regulon in microglia.

Pigs are important large animal models for biomedical research. Here, the authors construct a single-cell landscape of pig tissues, unravelling the phenotypic heterogeneity of blood endothelial cells in adipose tissues and the evolutionally conserved regulons of microglia in brains.

VThunter: a database for single-cell screening of virus target cells in the animal kingdom

Viral infectious diseases are a devastating and continuing threat to human and animal health. Receptor binding is the key step for viral entry into host cells. Therefore, recognizing viral receptors is fundamental for understanding the potential tissue tropism or host range of these pathogens. The rapid advancement of single-cell RNA sequencing (scRNA-seq) technology has paved the way for studying the expression of viral receptors in different tissues of animal species at single-cell resolution, resulting in huge scRNA-seq datasets. However, effectively integrating or sharing these datasets among the research community is challenging, especially for laboratory scientists. In this study, we manually curated up-to-date datasets generated in animal scRNA-seq studies, analyzed them using a unified processing pipeline, and comprehensively annotated 107 viral receptors in 142 viruses and obtained accurate expression signatures in 2 100 962 cells from 47 animal species. Thus, the VThunter database provides a user-friendly interface for the research community to explore the expression signatures of viral receptors. VThunter offers an informative and convenient resource for scientists to better understand the interactions between viral receptors and animal viruses and to assess viral pathogenesis and transmission in species. Database URL: https://db.cngb.org/VThunter/.

Single cell atlas for 11 non-model mammals, reptiles and birds

The availability of viral entry factors is a prerequisite for the cross-species transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Large-scale single-cell screening of animal cells could reveal the expression patterns of viral entry genes in different hosts. However, such exploration for SARS-CoV-2 remains limited. Here, we perform single-nucleus RNA sequencing for 11 non-model species, including pets (cat, dog, hamster, and lizard), livestock (goat and rabbit), poultry (duck and pigeon), and wildlife (pangolin, tiger, and deer), and investigated the co-expression of ACE2 and TMPRSS2. Furthermore, cross-species analysis of the lung cell atlas of the studied mammals, reptiles, and birds reveals core developmental programs, critical connectomes, and conserved regulatory circuits among these evolutionarily distant species. Overall, our work provides a compendium of gene expression profiles for non-model animals, which could be employed to identify potential SARS-CoV-2 target cells and putative zoonotic reservoirs.

A high-resolution cell atlas of the domestic pig lung and an online platform for exploring lung single-cell data

The genetically engineered pig is regarded as an optimal source of organ transplantation for humans and an excellent model for human disease research, given its comparable physiology to human beings. A myriad of single-cell RNA sequencing (scRNA-seq) data on humans has been reported, but such data on pigs are scarce. Here, we apply scRNA-seq technology to study the cellular heterogeneity of 3-month-old pig lungs, generating the single-cell atlas of 13,580 cells covering 16 major cell types. Based on these data, we systematically characterize the similarities and differences in the cellular cross-talk and expression patterns of respiratory virus receptors in each cell type of pig lungs compared with human lungs. Furthermore, we analyze pig lung xenotransplantation barriers and reported the cell-type expression patterns of 10 genes associated with pig-to-human immunobiological incompatibility and coagulation dysregulation. We also investigate the conserved transcription factors (TFs) and their candidate target genes and constructed five conserved TF regulatory networks in the main cell types shared by pig and human lungs. Finally, we present a comprehensive and openly accessible online platform, ScdbLung. Our scRNA-seq atlas of the domestic pig lung and ScdbLung database can guide pig lung research and clinical applicability.