Quality control in scRNA-Seq can discriminate pacemaker cells: the mtRNA bias.

Single-cell RNA-sequencing (scRNA-seq) provides high-resolution insights into complex tissues. Cardiac tissue, however, poses a major challenge due to the delicate isolation process and the large size of mature cardiomyocytes. Regardless of the experimental technique, captured cells are often impaired and some capture sites may contain multiple or no cells at all. All this refers to "low quality" potentially leading to data misinterpretation. Common standard quality control parameters involve the number of detected genes, transcripts per cell, and the fraction of transcripts from mitochondrial genes. While cutoffs for transcripts and genes per cell are usually user-defined for each experiment or individually calculated, a fixed threshold of 5% mitochondrial transcripts is standard and often set as default in scRNA-seq software. However, this parameter is highly dependent on the tissue type. In the heart, mitochondrial transcripts comprise almost 30% of total mRNA due to high energy demands. Here, we demonstrate that a 5%-threshold not only causes an unacceptable exclusion of cardiomyocytes but also introduces a bias that particularly discriminates pacemaker cells. This effect is apparent for our in vitro generated induced-sinoatrial-bodies (iSABs; highly enriched physiologically functional pacemaker cells), and also evident in a public data set of cells isolated from embryonal murine sinoatrial node tissue (Goodyer William et al. in Circ Res 125:379-397, 2019). Taken together, we recommend omitting this filtering parameter for scRNA-seq in cardiovascular applications whenever possible.

Insights into early ontogenesis: characterization of stress and development key genes of pikeperch (Sander lucioperca) in vivo and in vitro.

There are still numerous difficulties in the successful farming of pikeperch in the anthropogenic environment of various aquaculture systems, especially during early developmental steps in the hatchery. To investigate the physiological processes involved on the molecular level, we determined the basal expression patterns of 21 genes involved in stress and immune responses and early ontogenesis of pikeperch between 0 and 175 days post hatch (dph). Their transcription patterns most likely reflect the challenges of growth and feed conversion. The gene coding for apolipoprotein A (APOE) was strongly expressed at 0 dph, indicating its importance for yolk sac utilization. Genes encoding bone morphogenetic proteins 4 and 7 (BMP4, BMP7), creatine kinase M (CKM), and SRY-box transcription factor 9 (SOX9) were highly abundant during the peak phases of morphological changes and acclimatization processes at 4-18 dph. The high expression of genes coding for peroxisome proliferator-activated receptors alpha and delta (PPARA, PPARD) at 121 and 175 dph, respectively, suggests their importance during this strong growth phase of juvenile stages. As an alternative experimental model to replace further in vivo investigations of ontogenetically important processes, we initiated the first approach towards a long-lasting primary cell culture from whole pikeperch embryos. The present study provides a set of possible biomarkers to support the monitoring of pikeperch farming and provides a first basis for the establishment of a suitable cell model of this emerging aquaculture species.

Early response of salmonid head-kidney cells to stress hormones and toll-like receptor ligands.

The functional spectrum of the teleostean head kidney covers haematopoietic, immune and endocrine signalling pathways with physiological effects that are likely to conflict if activated at the same time. An in vivo experiment on the salmonid fish maraena whitefish (Coregonus maraena) revealed that the head kidney shows a remarkably strong response after injection of Aeromonas salmonicida within 48 h. In order to investigate the potential influence of endocrine signalling on the initiation of immune responses, we established a primary culture of head-kidney cells of maraena whitefish. For the characterisation of this model system, we used flow cytometry complemented with an extensive panel of immunological/haematological and stress-physiological/neuroendocrinological qPCR assays. More than one third of the cells expressed the characteristic signature of myeloid cells, while more than half of the cells expressed those genes typical for lymphocytes and monocytes. In parallel, we quantified the expression of genes encoding endocrine receptors and identified ADRA2D as by far the most highly expressed adrenergic-receptor gene in head-kidney cells. The stimulation of the head-kidney cells with toll-like receptor ligands induced the expression of typical immune genes (IL1B, CXCL8, TNF, SAA) after only 1 h. The incubation with the stress hormones cortisol, adrenaline and noradrenaline also had an immune-activating effect, though less pronounced. However, cortisol had the strongest suppressive effect on the stimulation-induced immune response, while adrenaline exerted a comparably weaker effect and noradrenaline was almost ineffective. Moreover, we found that cortisol reduced the expression of genes coding for adrenergic and some glucocorticoid receptors, while noradrenaline increased it. In conclusion, the primary head-kidney cells of maraena whitefish reflect the immunological and neuroendocrinological diversity of the entire organ. This in vitro system allowed thus identifying the correlative changes between the activities of hormones and immune factors in salmonid fish in order to contribute to a better understanding of the regulation circuit between stress and immune defence.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls.

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid ß-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

Establishment and Characterization of an SV40 Large T Antigen-Transduced Porcine Colonic Epithelial Cell Line.

Continuous cell lines have become indispensable tools that have enabled investigations into cellular mechanisms by increasing experimental reproducibility and sample availability, and decreasing the use of experimental animals. To facilitate studies of epithelial barrier function of the porcine colon, we aimed to establish an epithelial cell line with an extended replicative capacity. Cells were isolated from the proximal colon of a 3-week-old piglet and transduced using a recombinant retroviral vector construct containing the simian virus 40 large T antigen (SV40 TAg). We established a clonal epithelial cell line, referred to as PoCo83-3, that stably expressed the SV40 TAg, verified at mRNA and protein levels. PoCo83-3 showed epithelial cell-specific features, such as cobblestone-like morphology, dome structure formation, the presence of apical microvilli, and the expression of keratin 18, E-cadherin and the tight junction-associated proteins zonula occludens-1, occludin, and claudin-1. To validate PoCo83-3 as an in vitro model in epithelial barrier research, proinflammatory cytokine-inducible alterations in barrier integrity were demonstrated by incubating the cells with TNF-α and IFN-γ for 48 h. These cytokine treatments promoted a decreased transepithelial electrical resistance. In summary, PoCo83-3 exhibited an extended life span and a differentiated phenotype while maintaining epithelial characteristics. Based on these results, we present this cell line as a valuable in vitro model for investigations of epithelial barrier function in the porcine colon.

Immediate and long-term transcriptional response of hind muscle tissue to transient variation of incubation temperature in broilers.

BACKGROUND: In oviparous species accidental variation of incubation temperatures may occur under natural conditions and mechanisms may have evolved by natural selection that facilitate coping with these stressors. However, under controlled artificial incubation modification of egg incubation temperature has been shown to have a wide-ranging impact on post-hatch development in several poultry species. Because developmental changes initiated in-ovo can affect poultry production, understanding the molecular routes and epigenetic alterations induced by incubation temperature differences may allow targeted modification of phenotypes. RESULTS: In order to identify molecular pathways responsive to variable incubation temperature, broiler eggs were incubated at a lower or higher temperature (36.8 °C, 38.8 °C) relative to control (37.8 °C) over two developmental intervals, embryonic days (E) 7-10 and 10-13. Global gene expression of M. gastrocnemius was assayed at E10, E13, and slaughter age [post-hatch day (D) 35] (6 groups; 3 time points; 8 animals each) by microarray analysis and treated samples were compared to controls within each time point. Transcript abundance differed for between 113 and 738 genes, depending on treatment group, compared to the respective control. In particular, higher incubation temperature during E7-10 immediately affected pathways involved in energy and lipid metabolism, cell signaling, and muscle development more so than did other conditions. But lower incubation temperature during E10-13 affected pathways related to cellular function and growth, and development of organ, tissue, and muscle as well as nutrient metabolism pathways at D35. CONCLUSION: Shifts in incubation temperature provoke specific immediate and long-term transcriptional responses. Further, the transcriptional response to lower incubation temperature, which did not affect the phenotypes, mediates compensatory effects reflecting adaptability. In contrast, higher incubation temperature triggers gene expression and has long-term effects on the phenotype, reflecting considerable phenotypic plasticity.

Altered incubation temperatures between embryonic Days 7 and 13 influence the weights and the mitochondrial respiratory and enzyme activities in breast and leg muscles of broiler embryos.

Altering incubation temperature during embryogenesis has an impact on chicken embryo growth, but the underlying molecular mechanisms are not understood; the present study was performed to address these changes. Broiler eggs were incubated at low (36.8°C), control (37.8°C), and high (38.8°C) temperatures between Embryonic Day (ED) 7 and 10 or ED 10 and 13, which cover critical periods of embryonic myogenesis. The embryos were then dissected immediately after treatment on ED 10 or 13 to assess body, liver, and heart weights as well as to analyze breast and leg muscle fibers for their mitochondrial respiratory activity (MRA). Breast muscle samples were additionally used to evaluate the activity of enzymes involved in energy metabolism and cell-cycle progression. ED-10 embryos incubated at 38.8°C showed elevated weights (body, liver, and heart), MRA, and activities of lactate dehydrogenase and cytochrome oxidase compared to the ED-10 embryos incubated at 36.8°C. Similarly, the ED-13 embryos incubated at 38.8°C showed elevated body weight, MRA, and activities of glycogen phosphorylase, phosphofructokinase, and cytochrome oxidase compared to their 36.8°C counterparts. Embryos incubated at the normal temperature (37.8°C), however, showed variable differences from those incubated at 38.8°C versus 36.8°C. Cell-cycle enzyme activities were not impacted by the different temperature treatments. Thus, an increase or decrease in the incubation temperature during embryonic broiler myogenesis results in altered embryo activity, muscle energy metabolism, and activity-dependent muscle growth.

Cloning and characterization of the proximal promoter region of rainbow trout (Oncorhynchus mykiss) interleukin-6 gene.

Interleukin-6 (IL6) is a pleiotropic cytokine with important immunoregulatory functions. Its expression is inducible in immune cells and tissues of several fish species. We also found that IL6 mRNA abundance was significantly increased in spleen, liver, and gill of rainbow trout after experimental infection with Aeromonas salmonicida. Genomic DNA sequences of IL6 orthologs from three salmonid species revealed a conserved exon/intron structure and a high overall nucleotide identity of >88%. To uncover key mechanisms regulating IL6 expression in salmonid fish, we amplified a fragment of the proximal IL6 promoter from rainbow trout and identified in-silico conserved binding sites for NF-κB and CEBP. The activity of this IL6 promoter fragment was analyzed in the established human embryonic kidney line HEK-293. Luciferase- and GFP-based reporter systems revealed that the proximal IL6 promoter is activated by Escherichia coli. Essentially, both reporter systems proved that NF-κB p50, but not NF-κB p65 or CEBP, activates the IL6 promoter fragment. Truncation of this fragment caused a significant decrease in IL6 promoter activation. This characterization of the proximal promoter of the IL6-encoding gene provides basic knowledge about the IL6 gene expression in rainbow trout.

Haplotype-resolved and near-T2T genome assembly of the African catfish (Clarias gariepinus).

Airbreathing catfish are stenohaline freshwater fish capable of withstanding various environmental conditions and farming practices, including breathing atmospheric oxygen. This unique ability has enabled them to thrive in semi-terrestrial habitats. However, the genomic mechanisms underlying their adaptation to adverse ecological environments remain largely unexplored, primarily due to the limited availability of high-quality genomic resources. Here, we present a haplotype-resolved and near telomere-to-telomere (T2T) genome assembly of the African catfish (Clarias gariepinus), utilizing Oxford Nanopore, PacBio HiFi, Illumina and Hi-C sequencing technologies. The primary assembly spans 969.62 Mb with only 47 contigs, achieving a contig N50 of 33.71 Mb. Terminal telomeric signals were detected in 22 of 47 contigs, suggesting T2T assembled chromosomes. BUSCO analysis confirmed gene space completeness of 99% against the Actinopterygii dataset, highlighting the high quality of the assembly. Genome annotation identified 25,655 protein-coding genes and estimated 43.94% genome-wide repetitive elements. This data provides valuable genomic resources to advance aquaculture practices and to explore the genomic underpinnings of the ecological resilience of airbreathing catfish and related teleosts.

KRT8, FAF1 and PTH1R gene polymorphisms are associated with leg weakness traits in pigs.

The liability to lesions of dysfunctions of bone and joints in pigs, summarized as leg weakness and mostly expressed as osteochondrosis, is an animal welfare and economic issue in pig production. The objective of this study was to identify polymorphisms in the functional and positional candidate genes keratin 8 (KRT8), Fas-associated factor 1 (FAF1) and parathyroid hormone type I receptor (PTH1R) and to evaluate their association with leg weakness traits. Therefore, osteochondrosis lesions were scored in animals of a Duroc × Pietrain F2 population (DuPi; n = 310) and commercial herds of the breed Large White (n = 298). In addition, bone mineralization traits were observed in DuPi population. SNPs were identified in genes KRT8 (g.8,039G > A), FAF1 (g.380,914T > C) and PTH1R (c.1,672C > T). KRT8 showed significant association with bone mineral density and content (P ≤ 0.05). FAF1 was association with OC lesions score of all joints inspected (P ≤ 0.05). PTH1R showed significant dominance effects on OC lesion scores of the distal femur articular cartilage (P = 0.01) and epiphysis of the distal ulna (P = 0.05) as well as sums of scores of all joints (OCsum, P = 0.04) and assignment to groups of either severely or gently affected animals (OCcat, P = 0.01). This study reveals clear genetic-statistical evidence for a link of KRT8, FAF1 and PTH1R with some of leg weakness related traits in pigs.

Experimental Handling Challenges Result in Minor Changes in the Phagocytic Capacity and Transcriptome of Head-Kidney Cells of the Salmonid Fish Coregonus maraena.

Aquaculture management involves regular handling procedures, but these can evoke stress responses in farmed fish. We compiled an extensive list of published parameters that indicate the most likely handling-induced physiological deviations from the norm. However, since these parameters are based almost exclusively on studies of rainbow trout and Atlantic salmon, we conducted a handling-challenge experiment with maraena whitefish (Coregonus maraena). This salmonid fish was sampled at either 3 or 24 h after a single 1-min handling or after 10 days of daily repeated 1-min handling. The cortisol levels were strongly elevated in some individuals at 3 h after the single handling challenge, but these elevations were not significantly different between the challenged and control cohorts. The phagocytic capacity of myeloid head-kidney cells stimulated with fluorophore-labeled, inactivated Aeromonas salmonicida was significantly decreased in maraena whitefish at 3 h after the handling challenge compared to control fish. Microarray analysis of head-kidney samples from the challenged and control fish revealed 12 differentially expressed genes at 3 h and 70 at 24 h after the single handling episode, but only 5 differentially expressed genes after 10 days of repeated daily handling. The identified genes were assigned to numerous stress- and immune-relevant functional pathways, including "glucocorticoid receptor signaling" (3 h post-challenge), "HIF1A signaling" (24 h post-challenge), or "complement system" (10 days of repeated challenge). Our data reveal the tight interconnection of immune and stress pathways in the head kidney of maraena whitefish and corroborate several parameters previously found regulated in other tissues of handling-stressed rainbow trout. These findings indicate that handling may compromise the health and welfare of maraena whitefish in aquaculture.

Evaluation of blood cell viability rate, gene expression, and O-GlcNAcylation profiles as indicative signatures for fungal stimulation of salmonid cell models.

Fungal diseases of fish are a significant economic problem in aquaculture. Using high-throughput expression analysis, we identified potential transcript markers in primary head kidney and secondary embryonic cells from salmonid fish after stimulation with the inactivated fungi Mucor hiemalis and Fusarium aveneacium and with purified fungal molecular patterns. The transcript levels of most of the 45 selected genes were altered in head-kidney cells after 24 h of stimulation with fungal antigens. Stimulation with the inactivated fungus M. hiemalis induced the most pronounced transcriptional changes, including the pathogen receptor-encoding genes CLEC18A and TLR22, the cytokine-encoding genes IL6 and TNF, and the gene encoding the antimicrobial peptide LEAP2. In parallel, we analyzed the total GlcNAcylation status of embryonic salmonid cells with or without stimulation with inactivated fungi. O-GlcNAcylation modulates gene expression, intracellular protein, and signal activity, but we detected no significant differences after a 3-h stimulation. A pathway analysis tool identified the "apoptosis of leukocytes" based on the expression profile 24 h after fungal stimulation. Fluorescence microscopy combined with flow cytometry revealed apoptosis in 50 % of head-kidney leukocytes after 3 h stimulation with M. hiemalis, but this level decreased by > 5% after 24 h of stimulation. The number of apoptotic cells significantly increased in all blood cells after a 3-h stimulation with fungal molecular patterns compared to unstimulated controls. This in vitro approach identified transcript-based parameters that were strongly modulated by fungal infections of salmonid fish.

Annotation of novel transcripts putatively relevant for bovine fat metabolism.

Two bovine transcripts encoded by the interleukin-1 receptor-associated kinase 1 (IRAK1) gene and the locus LOC618944 predicted as similar to human chromosome 6 open reading frame 52 (C6orf52) gene had indicated divergent expression in bovine skeletal muscle containing different amount of intramuscular fat in a pilot screening experiment. However, for both loci any role in the regulation of energy or fat metabolism is not yet described. In this study, we validated and refined gene structure, screened for mRNA splice variants and analyzed the tissue-specific gene expression patterns of both loci as a prerequisite to elucidate their potential physiological function. Based on comparative sequence analysis, a new full-length gene model for the bovine IRAK1 gene was developed and confirmed experimentally. Expression of IRAK1 mRNA was found in a variety of tissues, and a splice variant was identified in skeletal muscle caused by an in-frame deleted segment of 210 bp affecting regions of intrinsic disorder in the respective protein. For the locus LOC618944, our data contributed to a revised gene model and its assignment to BTA23 (bovine chromosome 23) on the current bovine genome assembly supported by comparative similarity analysis between the bovine and human genomes and experimental data. Furthermore, we identified several splice variants in mammary gland, fat and skeletal muscle tissue and detected a highly similar processed pseudogene on BTA26. All transcript variants of LOC618944 detected in the analyzed tissues represent noncoding RNAs. For both loci, our results suggest yet undetected physiological functions in tissues relevant for fat or energy metabolism in cattle.

Comprehensive Characterization of Multitissue Expression Landscape, Co-Expression Networks and Positive Selection in Pikeperch.

Promising efforts are ongoing to extend genomics resources for pikeperch (Sander lucioperca), a species of high interest for the sustainable European aquaculture sector. Although previous work, including reference genome assembly, transcriptome sequence, and single-nucleotide polymorphism genotyping, added a great wealth of genomic tools, a comprehensive characterization of gene expression across major tissues in pikeperch still remains an unmet research need. Here, we used deep RNA-Sequencing of ten vital tissues collected in eight animals to build a high-confident and annotated trancriptome atlas, to detect the tissue-specificity of gene expression and co-expression network modules, and to investigate genome-wide selective signatures in the Percidae fish family. Pathway enrichment and protein-protein interaction network analyses were performed to characterize the unique biological functions of tissue-specific genes and co-expression modules. We detected strong functional correlations and similarities of tissues with respect to their expression patterns-but also significant differences in the complexity and composition of their transcriptomes. Moreover, functional analyses revealed that tissue-specific genes essentially play key roles in the specific physiological functions of the respective tissues. Identified network modules were also functionally coherent with tissues' main physiological functions. Although tissue specificity was not associated with positive selection, several genes under selection were found to be involved in hypoxia, immunity, and gene regulation processes, that are crucial for fish adaption and welfare. Overall, these new resources and insights will not only enhance the understanding of mechanisms of organ biology in pikeperch, but also complement the amount of genomic resources for this commercial species.

Effects of Chronic Hypoxia on the Immune Status of Pikeperch (Sander lucioperca Linnaeus, 1758).

Inadequate oxygen saturation can induce stress responses in fish and further affect their immunity. Pikeperch, recently introduced in intensive aquaculture, is suggested to be reared at nearly 100% DO (dissolved oxygen), yet this recommendation can be compromised by several factors including the water temperature, stocking densities or low circulation. Herein, we aimed to investigate the effect of low oxygen saturation of 40% DO (±3.2 mg/L) over 28 days on pikeperch farmed in recirculating aquaculture systems. The obtained data suggest that-although the standard blood and health parameters did not reveal any significant differences at any timepoint-the flow cytometric analysis identified a slightly decreased proportion of lymphocytes in the HK (head kidney) of fish exposed to hypoxia. This has been complemented by marginally downregulated expression of investigated immune and stress genes in HK and liver (including FTH1, HIF1A and NR3C1). Additionally, in the model of acute peritoneal inflammation induced with inactivated Aeromonas hydrophila, we observed a striking dichotomy in the sensitivity to the low DO between innate and adaptive immunity. Thus, while the mobilization of myeloid cells from HK to blood, spleen and peritoneal cavity, underlined by changes in the expression of key proinflammatory cytokines (including MPO, IL1B and TNF) was not influenced by the low DO, hypoxia impaired the influx of lymphocytes to the peritoneal niche in the later phases of the immune reaction. Taken together, our data suggest high robustness of pikeperch towards the low oxygen saturation and further encourage its introduction to the intensive aquaculture systems.

Cytogenetic anchoring of radiation hybrid and virtual maps of sheep chromosome X and comparison of X chromosomes in sheep, cattle, and human.

A comprehensive physical map was generated for Ovis aries chromosome X (OARX) based on a cytogenomics approach. DNA probes were prepared from bacterial artificial chromosome (BAC) clones from the CHORI-243 sheep library and were assigned to G-banded metaphase spreads via fluorescence in-situ hybridization (FISH). A total of 22 BACs gave a single hybridization signal to the X chromosome and were assigned out of 32 tested. The positioned BACs contained 16 genes and a microsatellite marker which represent new cytogenetically mapped loci in the sheep genome. The gene and microsatellite loci serve to anchor between the existing radiation hybrid (RH) and virtual sheep genome (VSG) maps to the cytogenetic OARX map, whilst the BACs themselves also serve as anchors between the VSG and the cytogenetic maps. An additional 17 links between the RH and cytogenetic maps are provided by BAC end sequence (BES) derived markers that have also been positioned on the RH map. Comparison of the map orders for the cytogenetic, RH, and virtual maps reveals that the orders for the cytogenetic and RH maps are most similar, with only one locus, represented by BAC CH243-330E18, mapping to relatively different positions. Several discrepancies, including an inverted segment are found when comparing both the cytogenetic and RH maps with the virtual map. These discrepancies highlight the value of using physical mapping methods to inform the process of future in silico map construction. A detailed comparative analysis of sheep, human, and cattle mapping data allowed the construction of a comparative map that confirms and expands the knowledge about evolutionary conservation and break points between the X chromosomes of the three mammalian species.

Indication of Premelanosome Protein (PMEL) Expression Outside of Pigmented Bovine Skin Suggests Functions Beyond Eumelanogenesis.

The premelanosome protein (PMEL) is important for fibril formation within melanosomes during vertebrate melanogenesis. Fibrils form a matrix for pigment deposition within pigmented tissues such as skin and hair. PMEL mutations are known to modulate eumelanic pigmentation in vertebrates. However, in bovines, PMEL mutations were also found to alter pheomelanic pigmentation resulting in coat color dilution. Furthermore, epistatic effects of a mutated PMEL allele were detected in the phenotypic expression of the bovine hair defect "rat-tail syndrome" (RTS) characterized by charcoal coat color and hair deformation. Reports about PMEL gene expression in non-pigmented tissues raised the hypothesis that there may be unknown functions of the PMEL protein beyond eumelanin deposition to PMEL fibrils. In our study, we analysed the PMEL protein expression in pigmented skin and non-pigmented bovine tissues (non-pigmented skin, thyroid gland, rumen, liver, kidney, and adrenal gland cortex). We found that a processed form of the bovine PMEL protein is expressed in pigmented as well as in non-pigmented tissues, which is in line with gene expression data from targeted RT-PCR and whole transcriptome RNAseq analysis. The PMEL protein is located in membranes and within the cytosol of epithelial cells. Based on our data from bovine tissues, we concluded that at least in cattle PMEL potentially has additional, yet unexplored functions, which might contribute to effects of PMEL mutations on pheomelanin coat color dilution and charcoal coat color in RTS animals. However, indication of PMEL protein in unpigmented cells and tissues will require further confirmation in the future, because there have been no confirmed reports before, which had detected bovine PMEL protein with specific antibodies either in pigmented or unpigmented tissue.

Integrative Cluster Analysis of Whole Hearts Reveals Proliferative Cardiomyocytes in Adult Mice.

The recent development and broad application of sequencing techniques at the single-cell level is generating an unprecedented amount of data. The different techniques have their individual limits, but the datasets also offer unexpected possibilities when utilized collectively. Here, we applied snRNA-seq in whole adult murine hearts from an inbred (C57BL/6NRj) and an outbred (Fzt:DU) mouse strain to directly compare the data with the publicly available scRNA-seq data of the tabula muris project. Explicitly choosing a single-nucleus approach allowed us to pin down the typical heart tissue-specific technical bias, coming up with novel insights on the mammalian heart cell composition. For our integrated dataset, cardiomyocytes, fibroblasts, and endothelial cells constituted the three main cell populations accounting for about 75% of all cells. However, their numbers severely differed between the individual datasets, with cardiomyocyte proportions ranging from about 9% in the tabula muris data to around 23% for our BL6 data, representing the prime example for cell capture technique related bias when using a conventional single-cell approach for these large cells. Most strikingly in our comparison was the discovery of a minor population of cardiomyocytes characterized by proliferation markers that could not be identified by analyzing the datasets individually. It is now widely accepted that the heart has an, albeit very restricted, regenerative potential. However there is still an ongoing debate where new cardiomyocytes arise from. Our findings support the idea that the renewal of the cardiomyocyte pool is driven by cytokinesis of resident cardiomyocytes rather than differentiation of progenitor cells. We thus provide data that can contribute to an understanding of heart cell regeneration, which is a prerequisite for future applications to enhance the process of heart repair.

Single-Nucleus Sequencing of an Entire Mammalian Heart: Cell Type Composition and Velocity.

: Analyses on the cellular level are indispensable to expand our understanding of complex tissues like the mammalian heart. Single-nucleus sequencing (snRNA-seq) allows for the exploration of cellular composition and cell features without major hurdles of single-cell sequencing. We used snRNA-seq to investigate for the first time an entire adult mammalian heart. Single-nucleus quantification and clustering led to an accurate representation of cell types, revealing 24 distinct clusters with endothelial cells (28.8%), fibroblasts (25.3%), and cardiomyocytes (22.8%) constituting the major cell populations. An additional RNA velocity analysis allowed us to study transcription kinetics and was utilized to visualize the transitions between mature and nascent cellular states of the cell types. We identified subgroups of cardiomyocytes with distinct marker profiles. For example, the expression of Hand2os1 distinguished immature cardiomyocytes from differentiated cardiomyocyte populations. Moreover, we found a cell population that comprises endothelial markers as well as markers clearly related to cardiomyocyte function. Our velocity data support the idea that this population is in a trans-differentiation process from an endothelial cell-like phenotype towards a cardiomyocyte-like phenotype. In summary, we present the first report of sequencing an entire adult mammalian heart, providing realistic cell-type distributions combined with RNA velocity kinetics hinting at interrelations.