Single-cell transcriptomic and chromatin accessibility analyses of dairy cattle peripheral blood mononuclear cells and their responses to lipopolysaccharide

Single-cell RNA sequencing characterization of Holstein cattle blood and milk immune cells during a chronic Staphylococcus aureus mastitis infection

Mastitis remains the most prevalent and costly disease to dairy producers. Granulocytes are the primary host innate immune cell responders during infectious mastitis. Here we examine three mid-lactation Holsteins challenged with ~ 150 CFU of Staphylococcus aureus (Newbould) that developed chronic mastitis as assessed by bacteria and somatic cell counts in a single quarter. Single-cell RNA-sequencing (scRNA-seq) of blood and milk cells identified immune cell populations of interest from both tissues, and the proportion of cell types recovered via scRNA-seq were highly similar to those recovered via flow cytometry. Granulocytes were the predominating cell type in both blood and milk samples; however granulocytes identified via scRNA-seq revealed several clusters comprised primarily of milk-derived cells. Milk-enriched granulocyte clusters were further investigated to identify gene signatures indicative of the granulocyte-specific localized immune responses in the mammary gland during chronic mastitis infection. Biological process enrichment analysis of gene signatures further revealed relevant networks such as granulocyte migration, myeloid cell differentiation, and inflammatory responses. In total, the work describes the immune landscape occurring at both peripheral and local sites of cattle with mastitis and identified important granulocyte-specific features of the localized immune response occurring during chronic infection.

Comparative Genomics and Epigenomics of Transcriptional Regulation

Transcriptional regulation in response to diverse physiological cues involves complicated biological processes. Recent initiatives that leverage whole genome sequencing and annotation of regulatory elements significantly contribute to our understanding of transcriptional gene regulation. Advances in the data sets available for comparative genomics and epigenomics can identify evolutionarily constrained regulatory variants and shed light on noncoding elements that influence transcription in different tissues and developmental stages across species. Most epigenomic data, however, are generated from healthy subjects at specific developmental stages. To bridge the genotype-phenotype gap, future research should focus on generating multidimensional epigenomic data under diverse physiological conditions. Farm animal species offer advantages in terms of feasibility, cost, and experimental design for such integrative analyses in comparison to humans. Deep learning modeling and cutting-edge technologies in sequencing and functional screening and validation also provide great promise for better understanding transcriptional regulation in this dynamic field.

Comparative epigenetics of domestic animals: focusing on DNA accessibility and its impact on gene regulation and traits

IMPORTANCE

Chromatin accessibility is vital for gene regulation, determining the ability of DNA-binding proteins to access the genomic regions and drive transcriptional activity, reflecting environmental changes. Although human and murine studies have advanced the understanding of chromatin dynamics, domestic animals remain comparatively underexplored despite their importance in agriculture and veterinary medicine. Investigating the accessibility of chromatin in these species is crucial for improving traits such as productivity, disease resistance, and environmental adaptation. This review assessed chromatin accessibility research in domestic animals, highlighting its significance in understanding and improving livestock traits.

OBSERVATIONS

This review outlines chromatin accessibility research in domestic animals, focusing on critical developmental processes, tissue-specific regulation, and economically significant traits. Advances in techniques, such as Assay for Transposase-Accessible Chromatin using sequencing, have enabled detailed mapping of regulatory elements, shedding light on epigenetic regulation of traits, such as muscle development and productivity. Comparative studies have uncovered conserved and species-specific cis-regulatory elements across multiple species. These findings offer insights into regulatory mechanisms that can enhance breeding strategies and animal management. In addition, high-throughput techniques, such as single-cell analysis and deep-learning models, have advanced the study of chromatin accessibility in lesser-studied species.

CONCLUSIONS AND RELEVANCE

Chromatin accessibility is crucial in gene regulation in domestic animals, influencing development, immune response, and productivity. Despite the progress, more comprehensive epigenomic datasets and cross-species analytical tools are needed to harness chromatin accessibility in domestic animal research. Understanding these mechanisms has practical applications in improving livestock traits, advancing breeding programs, and developing disease-resistant animals, highlighting the importance of integrating epigenetic and genomic tools for enhancing animal health and productivity.

Alterations of the IKZF1-IKZF2 tandem in immune cells of schizophrenia patients regulate associated phenotypes

Schizophrenia is a complex multifactorial disorder and increasing evidence suggests the involvement of immune dysregulations in its pathogenesis. We observed that IKZF1 and IKZF2, classic immune-related transcription factors (TFs), were both downregulated in patients' peripheral blood mononuclear cells (PBMCs) but not in their brain. We generated a new mutant mouse model with a reduction in Ikzf1 and Ikzf2 to study the impact of those changes. Such mice developed deficits in the three dimensions (positive-negative-cognitive) of schizophrenia-like phenotypes associated with alterations in structural synaptic plasticity. We then studied the secretomes of cultured PBMCs obtained from patients and identified potentially secreted molecules, which depended on IKZF1 and IKZF2 mRNA levels, and that in turn have an impact on neural synchrony, structural synaptic plasticity and schizophrenia-like symptoms in in vivo and in vitro models. Our results point out that IKZF1-IKZF2-dependent immune signals negatively impact on essential neural circuits involved in schizophrenia.

Cross-species single-cell analysis reveals divergence and conservation of peripheral blood mononuclear cells

BACKGROUND

Single-cell transcriptome sequencing (scRNA-seq) has revolutionized the study of immune cells by overcoming the limitations of traditional antibody-based identification and isolation methods. This advancement allows us to obtain comprehensive gene expression profiles from a diverse array of vertebrate species, facilitating the identification of various cell types. Comparative immunology across vertebrates presents a promising approach to understanding the evolution of immune cell types. In this study, we conducted a comparative transcriptome analysis of peripheral blood mononuclear cells (PBMCs) at the single-cell level across 12 species.

RESULTS

Our findings shed light on the cellular compositional features of PBMCs, spanning from fish to mammals. Notably, we identified genes that exhibit vertebrate universality in characterizing immune cells. Moreover, our investigation revealed that monocytes have maintained a conserved transcriptional regulatory program throughout evolution, emphasizing their pivotal role in orchestrating immune cells to execute immune programs.

CONCLUSIONS

This comprehensive analysis provides valuable insights into the evolution of immune cells across vertebrates.

Single cell RNA-sequencing of feline peripheral immune cells with V(D)J repertoire and cross species analysis of T lymphocytes

Introduction

The domestic cat (Felis catus) is a valued companion animal and a model for virally induced cancers and immunodeficiencies. However, species-specific limitations such as a scarcity of immune cell markers constrain our ability to resolve immune cell subsets at sufficient detail. The goal of this study was to characterize circulating feline T cells and other leukocytes based on their transcriptomic landscape and T-cell receptor repertoire using single cell RNA-sequencing.

Methods

Peripheral blood from 4 healthy cats was enriched for T cells by flow cytometry cell sorting using a mouse anti-feline CD5 monoclonal antibody. Libraries for whole transcriptome, αβ T cell receptor transcripts and γδ T cell receptor transcripts were constructed using the 10x Genomics Chromium Next GEM Single Cell 5' reagent kit and the Chromium Single Cell V(D)J Enrichment Kit with custom reverse primers for the feline orthologs.

Results

Unsupervised clustering of whole transcriptome data revealed 7 major cell populations - T cells, neutrophils, monocytic cells, B cells, plasmacytoid dendritic cells, mast cells and platelets. Sub cluster analysis of T cells resolved naive (CD4+ and CD8+), CD4+ effector T cells, CD8+ cytotoxic T cells and γδ T cells. Cross species analysis revealed a high conservation of T cell subsets along an effector gradient with equitable representation of veterinary species (horse, dog, pig) and humans with the cat. Our V(D)J repertoire analysis identified a subset of CD8+ cytotoxic T cells with skewed TRA and TRB gene usage, conserved TRA and TRB junctional motifs, restricted TRA/TRB pairing and reduced diversity in TRG junctional length. We also identified a public γδ T cell subset with invariant TRD and TRG chains and a CD4+ TEM-like phenotype. Among monocytic cells, we resolved three clusters of classical monocytes with polarization into pro- and anti-inflammatory phenotypes in addition to a cluster of conventional dendritic cells. Lastly, our neutrophil sub clustering revealed a larger mature neutrophil cluster and a smaller exhausted/activated cluster.

Discussion

Our study is the first to characterize subsets of circulating T cells utilizing an integrative approach of single cell RNA-sequencing, V(D)J repertoire analysis and cross species analysis. In addition, we characterize the transcriptome of several myeloid cell subsets and demonstrate immune cell relatedness across different species.

Advances in single-cell transcriptomics in animal research

Understanding biological mechanisms is fundamental for improving animal production and health to meet the growing demand for high-quality protein. As an emerging biotechnology, single-cell transcriptomics has been gradually applied in diverse aspects of animal research, offering an effective method to study the gene expression of high-throughput single cells of different tissues/organs in animals. In an unprecedented manner, researchers have identified cell types/subtypes and their marker genes, inferred cellular fate trajectories, and revealed cell‒cell interactions in animals using single-cell transcriptomics. In this paper, we introduce the development of single-cell technology and review the processes, advancements, and applications of single-cell transcriptomics in animal research. We summarize recent efforts using single-cell transcriptomics to obtain a more profound understanding of animal nutrition and health, reproductive performance, genetics, and disease models in different livestock species. Moreover, the practical experience accumulated based on a large number of cases is highlighted to provide a reference for determining key factors (e.g., sample size, cell clustering, and cell type annotation) in single-cell transcriptomics analysis. We also discuss the limitations and outlook of single-cell transcriptomics in the current stage. This paper describes the comprehensive progress of single-cell transcriptomics in animal research, offering novel insights and sustainable advancements in agricultural productivity and animal health.

Pathogen stimulations and immune cells synergistically affect the gene expression profile characteristics of porcine peripheral blood mononuclear cells

BACKGROUND

Pigs serve as a crucial source of protein in the human diet and play a fundamental role in ensuring food security. However, infectious diseases caused by bacteria or viruses are a major threat to effective global pig farming, jeopardizing human health. Peripheral blood mononuclear cells (PBMCs) are a mixture of immune cells that play crucial roles in immunity and disease resistance in pigs. Previous studies on the gene expression regulation patterns of PBMCs have concentrated on a single immune stimulus or immune cell subpopulation, which has limited our comprehensive understanding of the mechanisms of the pig immune response.

RESULTS

Here, we integrated and re-analyzed RNA-seq data published online for porcine PBMC stimulated by lipopolysaccharide (LPS), polyinosinic acid (PolyI:C), and various unknown microorganisms (EM). The results revealed that gene expression and its functional characterization are highly specific to the pathogen, identifying 603, 254, and 882 pathogen-specific genes and 38 shared genes, respectively. Notably, LPS and PolyI:C stimulation directly triggered inflammatory and immune-response pathways, while exposure to mixed microbes (EM) enhanced metabolic processes. These pathogen-specific genes were enriched in immune trait-associated quantitative trait loci (QTL) and eGenes in porcine immune tissues and were implicated in specific cell types. Furthermore, we discussed the roles of eQTLs rs3473322705 and rs1109431654 in regulating pathogen- and cell-specific genes CD300A and CD93, using cellular experiments. Additionally, by integrating genome-wide association studies datasets from 33 complex traits and diseases in humans, we found that pathogen-specific genes were significantly enriched for immune traits and metabolic diseases.

CONCLUSIONS

We systematically analyzed the gene expression profiles of the three stimulations and demonstrated pathogen-specific and cell-specific gene regulation across different stimulations in porcine PBMCs. These findings enhance our understanding of shared and distinct regulatory mechanisms of genetic variants in pig immune traits.

Single-Cell Sequencing Technology in Ruminant Livestock: Challenges and Opportunities

Advancements in single-cell sequencing have transformed the genomics field by allowing researchers to delve into the intricate cellular heterogeneity within tissues at greater resolution. While single-cell omics are more widely applied in model organisms and humans, their use in livestock species is just beginning. Studies in cattle, sheep, and goats have already leveraged single-cell and single-nuclei RNA-seq as well as single-cell and single-nuclei ATAC-seq to delineate cellular diversity in tissues, track changes in cell populations and gene expression over developmental stages, and characterize immune cell populations important for disease resistance and resilience. Although challenges exist for the use of this technology in ruminant livestock, such as the precise annotation of unique cell populations and spatial resolution of cells within a tissue, there is vast potential to enhance our understanding of the cellular and molecular mechanisms underpinning traits essential for healthy and productive livestock. This review intends to highlight the insights gained from published single-cell omics studies in cattle, sheep, and goats, particularly those with publicly accessible data. Further, this manuscript will discuss the challenges and opportunities of this technology in ruminant livestock and how it may contribute to enhanced profitability and sustainability of animal agriculture in the future.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

Single cell RNA-sequencing of feline peripheral immune cells with V(D)J repertoire and cross species analysis of T lymphocytes

Introduction

The domestic cat (Felis catus) is a valued companion animal and a model for virally induced cancers and immunodeficiencies. However, species-specific limitations such as a scarcity of immune cell markers constrain our ability to resolve immune cell subsets at sufficient detail. The goal of this study was to characterize circulating feline T cells and other leukocytes based on their transcriptomic landscape and T-cell receptor repertoire using single cell RNA-sequencing.

Methods

Peripheral blood from 4 healthy cats was enriched for T cells by flow cytometry cell sorting using a mouse anti-feline CD5 monoclonal antibody. Libraries for whole transcriptome, alpha/beta T cell receptor transcripts and gamma/delta T cell receptor transcripts were constructed using the 10x Genomics Chromium Next GEM Single Cell 5' reagent kit and the Chromium Single Cell V(D)J Enrichment Kit with custom reverse primers for the feline orthologs.

Results

Unsupervised clustering of whole transcriptome data revealed 7 major cell populations - T cells, neutrophils, monocytic cells, B cells, plasmacytoid dendritic cells, mast cells and platelets. Sub cluster analysis of T cells resolved naive (CD4+ and CD8+), CD4+ effector T cells, CD8+ cytotoxic T cells and gamma/delta T cells. Cross species analysis revealed a high conservation of T cell subsets along an effector gradient with equitable representation of veterinary species (horse, dog, pig) and humans with the cat. Our V(D)J repertoire analysis demonstrated a skewed T-cell receptor alpha gene usage and a restricted T-cell receptor gamma junctional length in CD8+ cytotoxic T cells compared to other alpha/beta T cell subsets. Among myeloid cells, we resolved three clusters of classical monocytes with polarization into pro- and anti-inflammatory phenotypes in addition to a cluster of conventional dendritic cells. Lastly, our neutrophil sub clustering revealed a larger mature neutrophil cluster and a smaller exhausted/activated cluster.

Discussion

Our study is the first to characterize subsets of circulating T cells utilizing an integrative approach of single cell RNA-sequencing, V(D)J repertoire analysis and cross species analysis. In addition, we characterize the transcriptome of several myeloid cell subsets and demonstrate immune cell relatedness across different species.

The single-cell transcriptome and chromatin accessibility datasets of peripheral blood mononuclear cells in Chinese holstein cattle

OBJECTIVES

This study was performed in the frame of a more extensive study dedicated to the integrated analysis of the single-cell transcriptome and chromatin accessibility datasets of peripheral blood mononuclear cells (PBMCs) with a large-scale GWAS of 45 complex traits in Chinese Holstein cattle. Lipopolysaccharide (LPS) is a crucial mediator of chronic inflammation to modulate immune responses. PBMCs include primary T and B cells, natural killer (NK) cells, monocytes (Mono), and dendritic cells (DC). How LPS stimulates PBMCs at the single-cell level in dairy cattle remains largely unknown.

DATA DESCRIPTION

We sequenced 30,756 estimated single cells and mapped 26,141 of them (96.05%) with approximately 60,075 mapped reads per cell after quality control for four whole-blood treatments (no, 2 h, 4 h, and 8 h LPS) by single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq). Finally, 7,107 (no), 9,174 (2 h), 6,741 (4 h), and 3,119 (8 h) cells were generated with ~ 15,000 total genes in the whole population. Therefore, the single-cell transcriptome and chromatin accessibility datasets in this study enable a further understanding of the cell types and functions of PBMCs and their responses to LPS stimulation in vitro.

Analysis of CircRNA Expression in Peripheral Blood of Holstein Cows in Response to Heat Stress

The present study aimed to identify key circRNAs and pathways associated with heat stress in blood samples of Holstein cows, which will provide new insights into the molecular mechanisms driving heat stress in cows. Hence, we evaluated changes in milk yield, rectal temperature, and respiratory rate of experimental cows between heat stress (summer) and non-heat stress (spring) conditions with two comparisons, including Sum1 vs. Spr1 (same lactation stage, different individuals, 15 cows per group) and Sum1 vs. Spr2 (same individual, different lactation stages, 15 cows per group). Compared to both Spr1 and Spr2, cows in the Sum1 group had a significantly lower milk yield, while rectal temperature and respiratory rate were significantly higher (p < 0.05), indicating that cows in the Sum1 group were experiencing heat stress. In each group, five animals were chosen randomly to undergo RNA-seq. The results reveal that 140 and 205 differentially expressed (DE) circRNAs were screened in the first and second comparisons, respectively. According to the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, these DE circRNAs were mainly enriched in five signaling pathways, including choline metabolism, the PI3K/AKT signaling pathway, the HIF-1 signaling pathway, the longevity-regulating pathway, and autophagy. Then, we obtained the top 10 hub source genes of circRNAs according to the protein-protein interaction networks. Among them, ciRNA1282 (HIF1A), circRNA4205 (NR3C1), and circRNA12923 (ROCK1) were enriched in multiple pathways and identified as binding multiple miRNAs. These key circRNAs may play an important role in the heat stress responses of dairy cows. These results provide valuable information on the involvement of key circRNAs and their expression pattern in the heat stress response of cows.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Gut microbiome is linked to functions of peripheral immune cells in transition cows during excessive lipolysis

BACKGROUND

Postpartum dairy cows experiencing excessive lipolysis are prone to severe immunosuppression. Despite the extensive understanding of the gut microbial regulation of host immunity and metabolism, its role during excessive lipolysis in cows is largely unknown. Herein, we investigated the potential links between the gut microbiome and postpartum immunosuppression in periparturient dairy cows with excessive lipolysis using single immune cell transcriptome, 16S amplicon sequencing, metagenomics, and targeted metabolomics.

RESULTS

The use of single-cell RNA sequencing identified 26 clusters that were annotated to 10 different immune cell types. Enrichment of functions of these clusters revealed a downregulation of functions in immune cells isolated from a cow with excessive lipolysis compared to a cow with low/normal lipolysis. The results of metagenomic sequencing and targeted metabolome analysis together revealed that secondary bile acid (SBA) biosynthesis was significantly activated in the cows with excessive lipolysis. Moreover, the relative abundance of gut Bacteroides sp. OF04 - 15BH, Paraprevotella clara, Paraprevotella xylaniphila, and Treponema sp. JC4 was mainly associated with SBA synthesis. The use of an integrated analysis showed that the reduction of plasma glycolithocholic acid and taurolithocholic acid could contribute to the immunosuppression of monocytes (CD14⁺MON) during excessive lipolysis by decreasing the expression of GPBAR1.

CONCLUSIONS

Our results suggest that alterations in the gut microbiota and their functions related to SBA synthesis suppressed the functions of monocytes during excessive lipolysis in transition dairy cows. Therefore, we concluded that altered microbial SBA synthesis during excessive lipolysis could lead to postpartum immunosuppression in transition cows. Video Abstract.

Comparative transcriptome analyses of immune responses to LPS in peripheral blood mononuclear cells from the giant panda, human, mouse, and monkey

Gram-negative bacteria are major pathogens that can cause illnesses in giant pandas. Lipopolysaccharides (LPS), components of Gram-negative bacteria, can activate immune responses in mammals (i.e., humans and mice) through recognition by toll-like receptors (TLRs). However, the giant pandas' immune response to LPS stimulation and the differences between the giant panda and other mammals are not fully known. In this study, we administrated peripheral blood mononuclear cells (PBMCs) from giant pandas, humans, C57BL/6 mice, and rhesus monkeys by LPS treatment at 6 h followed by RNA sequencing (RNA-seq), respectively, with control of non-stimulation. KEGG analyses of differentially expressed genes (DEGs) pathways indicated that LPS could activate the classic signaling pathway of NF-κB in PBMCs from those four tested species. Thus, similar to the other three species, NF-κB is an LPS-responsive regulator of innate immune responses in giant pandas. Furthermore, the expression patterns of adapter genes, inflammatory cytokine genes, chemokines, interferon genes, cytokine genes related to cell growth and development, costimulatory molecules, Th1/Th2 cytokine genes, Th17 cytokine genes, Th9, and Th22 cytokine genes were compared among giant pandas and three other species. Our data indicated that in addition to the similar expression patterns of certain genes among giant pandas and other species, the unique expression pattern response to LPS in giant pandas was also discovered. Furthermore, Th9, Th17, and Th22 cells might be involved in the response to LPS in giant pandas at this tested time point. This study reveals that LPS-induced immune responses have different sensitivities and response timelines in giant pandas compared with other mammals. This study facilitates further understanding of the role of the TLR signaling pathway and the immune system in giant pandas, which might be helpful for disease prevention and protection.