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.

Single-cell resolution landscape of equine peripheral blood mononuclear cells reveals diverse cell types including T-bet+ B cells.

BACKGROUND: Traditional laboratory model organisms represent a small fraction of the diversity of multicellular life, and findings in any given experimental model often do not translate to other species. Immunology research in non-traditional model organisms can be advantageous or even necessary, such as when studying host-pathogen interactions. However, such research presents multiple challenges, many stemming from an incomplete understanding of potentially species-specific immune cell types, frequencies, and phenotypes. Identifying and characterizing immune cells in such organisms is frequently limited by the availability of species-reactive immunophenotyping reagents for flow cytometry, and insufficient prior knowledge of cell type-defining markers. RESULTS: Here, we demonstrate the utility of single-cell RNA sequencing (scRNA-Seq) to characterize immune cells for which traditional experimental tools are limited. Specifically, we used scRNA-Seq to comprehensively define the cellular diversity of equine peripheral blood mononuclear cells (PBMC) from healthy horses across different breeds, ages, and sexes. We identified 30 cell type clusters partitioned into five major populations: monocytes/dendritic cells, B cells, CD3+PRF1+ lymphocytes, CD3+PRF1- lymphocytes, and basophils. Comparative analyses revealed many cell populations analogous to human PBMC, including transcriptionally heterogeneous monocytes and distinct dendritic cell subsets (cDC1, cDC2, plasmacytoid DC). Remarkably, we found that a majority of the equine peripheral B cell compartment is comprised of T-bet+ B cells, an immune cell subpopulation typically associated with chronic infection and inflammation in human and mouse. CONCLUSIONS: Taken together, our results demonstrate the potential of scRNA-Seq for cellular analyses in non-traditional model organisms and form the basis for an immune cell atlas of horse peripheral blood.

Single-cell analysis of nonhuman primate preimplantation development in comparison to humans and mice.

BACKGROUND: Genetic programs underlying preimplantation development and early lineage segregation are highly conserved across mammals. It has been suggested that nonhuman primates would be better model organisms for human embryogenesis, but a limited number of studies have investigated the monkey preimplantation development. In this study, we collect single cells from cynomolgus monkey preimplantation embryos for transcriptome profiling and compare with single-cell RNA-seq data derived from human and mouse embryos. RESULTS: By weighted gene-coexpression network analysis, we found that cynomolgus gene networks have greater conservation with human embryos including a greater number of conserved hub genes than that of mouse embryos. Consistently, we found that early ICM/TE lineage-segregating genes in monkeys exhibit greater similarity with human when compared to mouse, so are the genes in signaling pathways such as LRP1 and TCF7 involving in WNT pathway. Last, we tested the role of one conserved pre-EGA hub gene, SIN3A, using a morpholino knockdown of maternal RNA transcripts in monkey embryos followed by single-cell RNA-seq. We found that SIN3A knockdown disrupts the gene-silencing program during the embryonic genome activation transition and results in developmental delay of cynomolgus embryos. CONCLUSION: Taken together, our study provided new insight into evolutionarily conserved and divergent transcriptome dynamics during mammalian preimplantation development.

Novel application of single-cell next-generation sequencing for determination of intratumoral heterogeneity of canine osteosarcoma cell lines.

Osteosarcoma (OSA) is a highly aggressive and metastatic neoplasm of both the canine and human patient and is the leading form of osseous neoplasia in both species worldwide. To gain deeper insight into the heterogeneous and genetically chaotic nature of OSA, we applied single-cell transcriptome (scRNA-seq) analysis to 4 canine OSA cell lines. This novel application of scRNA-seq technology to the canine genome required uploading the CanFam3.1 reference genome into an analysis pipeline (10X Genomics Cell Ranger); this methodology has not been reported previously in the canine species, to our knowledge. The scRNA-seq outputs were validated by comparing them to cDNA expression from reverse-transcription PCR (RT-PCR) and Sanger sequencing bulk analysis of 4 canine OSA cell lines (COS31, DOUG, POS, and HMPOS) for 11 genes implicated in the pathogenesis of canine OSA. The scRNA-seq outputs revealed the significant heterogeneity of gene transcription expression patterns within the cell lines investigated (COS31 and DOUG). The scRNA-seq data showed 10 distinct clusters of similarly shared transcriptomic expression patterns in COS31; 12 clusters were identified in DOUG. In addition, cRNA-seq analysis provided data for integration into the Qiagen Ingenuity Pathway Analysis software for canonical pathway analysis. Of the 81 distinct pathways identified within the clusters, 33 had been implicated in the pathogenesis of OSA, of which 18 had not been reported previously in canine OSA.

Analysis of canine myeloid-derived suppressor cells (MDSCs) utilizing fluorescence-activated cell sorting, RNA protection mediums to yield quality RNA for single-cell RNA sequencing.

Fluorescence-activated cell sorting (FACS) is a branch of flow cytometry that allows for the isolation of specific cell populations that can then be further analyzed by single-cell RNA sequencing (scRNA-seq). When utilizing FACS for population isolation prior to sequencing, it is essential to consider the protection of RNA from RNase activity, environmental conditions, and the sorting efficiency to ensure optimum sample quality. This study aimed to optimize a previously published MDSC flow cytometry strategy to FACS sort canine Myeloid-Derived Suppressor Cells (MDSC) with various permutations of RNAlater ™ and RiboLock™ before and after FACS sorting. Concentrations of RNAlater™ greater than 2 % applied before flow analysis affected cell survival and fluorescence, whereas concentrations ≤ 2 % and time ≤ 4 h had little to no effect on cells. To shorten the procedural time and to enhance the sorting of rare populations, we used a primary PE-conjugated CD11b antibody and magnetic column. The combination of RiboLock™ pre- and post-sorting for FACS provided the best quality RNA as determined by the RNA integrity number (RIN ≥ 7) for scRNA-seq in a normal and dog and a dog with untreated oral melanoma dog. As proof of principle, we sequenced two samples, one from a normal dog another from a dog with untreated oral melanoma. Applying scRNA-Seq analysis using the 10X Genomic platform, we identified 6 clusters in the Seurat paired analysis of MDSC sorted samples. Two clusters, with the majority of the cells coming from the melanoma sample, had genes that were upregulated (> log2); these included MMP9, MMP1, HPGD, CPA3, and GATA3 and CYBB, CSTB, COX2, ATP6, and COX 17 for cluster 5 and 6 respectively. All genes have known associations with MDSCs. Further characterization using pathway analysis tools was not attempted due to the lower number of cells sequenced in the normal sample. The benefit deriving from the results of the study helped to gain data consistency when working with cells prone to RNase activity, and the scRNA-seq provided data showing transcriptional heterogeneity in MDSC populations and potentially identifying previously unreported or rare cell populations.

Oocyte IVM or vitrification significantly impairs DNA methylation patterns in blastocysts as analysed by single-cell whole-genome methylation sequencing.

To explore the mechanisms leading to the poor quality of IVF blastocysts, the single-cell whole-genome methylation sequencing technique was used in this study to analyse the methylation patterns of bovine blastocysts derived from invivo, fresh (IVF) or vitrified (V_IVF) oocytes. Genome methylation levels of blastocysts in the IVF and V_IVF groups were significantly lower than those of the invivo group (P<0.05). In all, 1149 differentially methylated regions (DMRs) were identified between the IVF and invivo groups, 1578 DMRs were identified between the V_IVF and invivo groups and 151 DMRs were identified between the V_IVF and IVF groups. For imprinted genes, methylation levels of insulin-like growth factor 2 receptor (IGF2R) and protein phosphatase 1 regulatory subunit 9A (PPP1R9A) were lower in the IVF and V_IVF groups than in the invivo group, and the methylation level of paternally expressed 3 (PEG3) was lower in the V_IVF group than in the IVF and invivo groups. Genes with DMRs between the IVF and invivo and the V_IVF and IVF groups were primarily enriched in oocyte maturation pathways, whereas DMRs between the V_IVF and invivo groups were enriched in fertilisation and vitrification-vulnerable pathways. The results of this study indicate that differences in the methylation of critical DMRs may contribute to the differences in quality between invitro- and invivo-derived embryos.

Identification of Pro-Fibrotic Macrophage Populations by Single-Cell Transcriptomic Analysis in West Highland White Terriers Affected With Canine Idiopathic Pulmonary Fibrosis.

Canine idiopathic pulmonary fibrosis (CIPF) affects old dogs from the West Highland white terrier (WHWT) breed and mimics idiopathic pulmonary fibrosis (IPF) in human. The disease results from deposition of fibrotic tissue in the lung parenchyma causing respiratory failure. Recent studies in IPF using single-cell RNA sequencing (scRNA-seq) revealed the presence of profibrotic macrophage populations in the lung, which could be targeted for therapeutic purpose. In dogs, scRNA-seq was recently validated for the detection of cell populations in bronchoalveolar lavage fluid (BALF) from healthy dogs. Here we used the scRNA-seq to characterize disease-related heterogeneity within cell populations of macrophages/monocytes (Ma/Mo) in the BALF from five WHWTs affected with CIPF in comparison with three healthy WHWTs. Gene set enrichment analysis was also used to assess pro-fibrotic capacities of Ma/Mo populations. Five clusters of Ma/Mo were identified. Gene set enrichment analyses revealed the presence of pro-fibrotic monocytes in higher proportion in CIPF WHWTs than in healthy WHWTs. In addition, monocyte-derived macrophages enriched in pro-fibrotic genes in CIPF compared with healthy WHWTs were also identified. These results suggest the implication of Ma/Mo clusters in CIPF processes, although, further research is needed to understand their role in disease pathogenesis. Overexpressed molecules associated with pulmonary fibrosis processes were also identified that could be used as biomarkers and/or therapeutic targets in the future.

Make it spin: individual trapping of sperm for analysis and recovery using micro-contact printing.

In this article, we describe the development of a high throughput platform to spatially manipulate viable sperm for motility measurements and recovery of the best single sperm for fertilization purposes. Micro-contact printing was used to pattern islands of adhesive proteins (fibronectin) separated by sperm repellent species (Pluronic acid F-127) on commercially available polystyrene substrates. Following washing, arrays of viable single sperm were captured onto the islands demonstrating for the first time that sperm can be trapped by micro-contact printing with patterning efficiency of 90% while retaining 100% viability. These were then subjected to motility analysis whilst remaining spatially confined to the islands. Single sperm motility was assessed (n = 37) by software analysis measuring the number of rotations per second (degrees s⁻¹). The assignment of array coordinates allows the more active single sperm to be easily identified and recovered by a simple micromanipulator pipette aspiration step with automated possibility for assisted reproductive technologies or further quality correlation analysis. Taken together, we show for the first time a technique to simultaneously screen thousands of viable single sperm for motility assessment while retaining the ability for single species recovery for enhanced fertilization purposes.

Zona pellucida birefringence correlates with developmental capacity of bovine oocytes classified by maturational environment, COC morphology and G6PDH activity.

In the present study we aimed to analyse structural changes during in vitro maturation of the bovine zona pellucida (ZP) by scanning electron microscopy (SEM) ands zona pellucida birefringence (ZPB). Here we show that alterations during in vitro maturation invasively analysed by SEM are reflected in ZPB. In vivo-matured oocytes displayed significantly lower birefringence parameters and significantly higher blastocyst rates compared with in vitro-derived oocytes (39.1% vs 21.6%). The same was observed for in vitro-matured oocytes with cumulus-oocyte complex (COC) Quality 1 (Q1) compared with Q3-COCs with respect to zona birefringence and developmental capacity. Immature oocytes with Q1-COCs displayed higher ZPB values and a higher developmental capacity to the blastocyst stage (27.7% vs 16.9%) compared with immature Q3-COCs. Considering in vitro-matured oocytes, only those with Q1-COC showed a trend for ZPB similar to in vivo-matured oocytes. Therefore, a decreasing trend for ZPB during in vitro maturation seems to be typical for high-quality oocytes and successful cytoplasmic maturation. In accordance, fully-grown immature oocytes reached significantly higher blastocyst rates (32.0% vs 11.5%) and lower ZPB values compared with still-growing ones. In conclusion, we successfully evaluated the applicability of zona imaging to bovine oocytes: alterations during in vitro maturation invasively analysed by scanning electron microscopy were reflected in the birefringence of the zona pellucida of bovine oocytes affecting developmental capacity at the same value. Therefore ZPB measurement by live zona imaging has potential to become a new tool to assess correctness of in vitro maturation and to predict developmental competence.