Skip to main content

Transcriptome dynamics in developing testes of domestic cats and impact of age on tissue resilience to cryopreservation

Abstract

Background

Fundamental knowledge of cellular and molecular mechanisms in developing testicular tissues is critical to better understand gonadal biology and responses to non-physiological conditions. The objective of our study was to (1) analyze transcriptome dynamics in developing testis of the domestic cat and (2) characterize age effects on the initial response of the tissue to vitrification. Tissues from adult and juvenile cats were processed for histology, DNA integrity, and RNA sequencing analyses before and after vitrification.

Results

Transcriptomic findings enabled to further characterize juvenile period, distinguishing between early and late juvenile tissues. Changes in gene expression and functional pathways were extensive from early to late juvenile to adult development stages. Additionally, tissues from juvenile animals were more resilient to vitrification compared to adult counterparts, with early juvenile sample responding the least to vitrification and late juvenile sample response being closest to adult tissues.

Conclusions

This is the first study reporting comprehensive datasets on transcriptomic dynamic coupled with structural analysis of the cat testis according to the age and exposure to cryopreservation. It provides a comprehensive network of functional terms and pathways that are affected by age in the domestic cat and are either enriched in adult or juvenile testicular tissues.

Peer Review reports

Background

Mammalian puberty is a progressive process in which testis undergoes dramatic developmental and structural changes, and involves complex hormonal and molecular modulation to accomplish both somatic cell proliferation/maturation and the initiation of spermatogenesis [1]. Large differences exist between species in hormonal control of puberty and onset of spermatogenesis, as well as in regulation of spermatogenesis and steroidogenesis [2, 3]. Recent studies in rodents and humans using single-cell RNA sequencing (scRNA-seq) technology enriched our knowledge on the processes happening in testis between infancy and adulthood, and further emphasized high species-specificity and value of animal models [4, 5].

The domestic cat is an essential model for biomedical research as well as conservation of endangered felids [6, 7]. Cat spermatogenic function reaches maturity at the age of 8 to 10 months, with initial activation and first signs of spermatogenesis at 5 to 6 months [8, 9]. This is fundamentally different from the mouse model (spermatogonia begin to differentiate shortly after birth resulting in a synchronous first wave of spermatogenesis) and closer to humans, where spermatogonia are maintained in an undifferentiated state prior to the initiation of puberty [3].

In addition to extensive morphological studies of domestic cat testicular tissue [10, 11], the expression of specific genes has also been investigated in developing testis [12, 13]. However, we still lack deeper knowledge on molecular processes happening during testis maturation in cats, as well as the understanding of differences between immature (neonatal), maturing (pre-, peri-, pubertal), and fully mature (adult) testicular tissue. A deeper understanding of the molecular events happening in testis from infancy to adulthood would contribute greatly to development of male fertility preservation and fertility control in domestic cats and wild felids. Whole transcriptome of testicular tissue during its maturation has been sequenced in several mammalian species, including pig [14], Mongolian horse [15], domestic yak [16], human [17] and mouse [18]. Further studies using scRNA-seq were also recently performed for developing testis in human and mouse [4, 5]. For the domestic cat, there has been no study so far that looked on the whole transcriptome of testicular tissue at different stages of its maturation. The only existing RNA-seq study on cat testicular tissue used three adult males over 2 years old and focused on comparison with sterile hybrids [19].

Sperm cryopreservation remains the standard approach for preserving male fertility in many species, including humans. However, when sperm banking is not possible, preservation of small fragments of testicular tissue offers an alternate way for fertility preservation [20, 21]. Preservation of testicular tissue can benefit not only prepubertal individuals who do not yet produce mature sperm, but also post-pubertal and adult patients and animals who may be azoospermic at the moment of semen collection [22, 23]. A recent report on the cryopreservation of testicular tissue for patients across several centers showed that the age of individuals ranged from 5 months all the way to 35 years old [24]. It is therefore essential to develop testicular tissue preservation protocols optimized for different age groups.

Recently, our group observed a better recovery of immature testicular tissues compared to adult in domestic cat after preservation protocol using microwave-assisted drying [25]. Different separate studies also show that testicular tissue from immature cats tend to be more resilient to cryopreservation protocols compared to adults [26,27,28]. However, there is no study comparing age groups side by side. We also lack a deeper understanding of the cellular mechanisms occurring in testicular tissues in response to vitrification and warming. RNA-seq has been used to explore the effect of cryopreservation on mammalian sperm cells in ruminants [29], boars [30, 31] and giant panda [32], but no studies looked at transcriptome changes in whole testicular tissue. There is evidence that the intrinsic response of cells to cryopreservation is different depending upon whether the cells are part of a tissue or whether they are isolated in a cell suspension [21]. Scaling up of cryopreservation from a microscopic cellular level to a macroscopic tissue level will introduce heat and mass transfer phenomena [33]. Additionally, freeze-thaw procedures for tissue products must result in recovery of both cell viability and tissue structure [33]. Thus, getting information on the whole transcriptome changes in the tissue during preservation is essential for understanding the mechanisms of stress response and recovery of cells as a system. For instance, our group recently demonstrated how RNA-seq and resulting functional networks can help to understand the response of ovarian tissue to various stresses caused by preservation [34].

The objective of our study was to (1) analyze transcriptome dynamics in developing testis of the domestic cat and (2) characterize age effects on the initial response of the tissue to vitrification. RNA-seq was performed on adult and juvenile cats (juvenile being defined as the development period after the infant stage but prior to full adult sexual maturity [4]) using fresh and vitrified whole testicular tissues. Because this is the first study to look at the whole transcriptome of the domestic cat developing testis, our study focused mainly on characterization of global changes in transcriptome and functional pathways throughout testis maturation. We included conventional methods of histology and DNA integrity analysis to be able to relate our transcriptomic data to the available knowledge.

Results

Summary of the acquired dataset

Testicular tissues were obtained from five adult and five juvenile male cats and divided into the following four groups: adult, fresh tissue (AF), adult, vitrified/warmed tissue (AV), juvenile, fresh tissue (JF) and juvenile, vitrified/warmed tissue (JV; Fig. 1). Samples from all four groups were used for histomorphology, TUNEL and RNA-seq analyses. Whole transcriptome of 20 samples was sequenced with one library per each sample, five biological replicates per group, 30 million read depth per sample and 150 bp paired-end read length. The acquired sequence data in fastq format is deposited to NCBI Sequence Read Archive; BioProject accession number is PRJNA741252. Data retrieved after differential gene expression analysis is available in additional files. Data used to create networks is available in the interactive web session view for each network.

Fig. 1
figure 1

Experimental design

Histology and transcriptome dynamics in fresh testicular tissue from adult and juvenile domestic cats

Figure 2 A presents histological overview of samples used for RNA-sEq. All samples from adult cats had Johnsen score above 7, while scores of juvenile samples ranged from 3 to 5 (Table S1, [35]). Unsupervised hierarchical clustering of all genes separated adult and juvenile samples, showing clear transcriptomic division of these age groups (Fig. 2B). Volcano plot in Fig. 2 C represents the results of differential expression (DE) analysis in adult vs. juvenile samples; a total of 8,732 differentially expressed genes (DEGs; Additional file 1). Principle component analysis revealed further division of juvenile samples into seemingly two age groups which we labeled as late (JF1, JF2) and early (JF3, JF4. JF5) juvenile ages (Fig. 2D). DE analysis then was performed using the new three groups with the following comparison pairs: adult vs. late juvenile (Additional file 2), late vs. early juvenile (Additional file 3) and adult vs. early juvenile. Out of all DEGs, 824 were shared between all three comparison pairs (Fig. 2E). All except 2 shared genes had the same expression direction being either continuously upregulated (653 protein coding genes, 140 lncRNAs) or downregulated (11 protein coding genes) from early to late juvenile to adult (Additional file 4). The 19 genes shared only between comparison pairs of late vs. early juvenile and adult vs. late juvenile all had different expression directions (Fig. 2E, Additional file 4).

Fig. 2
figure 2

Histology and transcriptome dynamics in fresh testicular tissue of adult and juvenile domestic cats. A Histology of testicular tissue samples used for transcriptome analysis; AF6 represents the histology of adult samples. Here and further, sample IDs represent age (A, adult; J, juvenile), tissue condition (F, fresh) and animal (1 to 10). Sg, spermatogonia; Sc, spermatocytes; Sd, round spermatids, Sz, spermatozoa. Scale bar: 50 μm. B Heatmap of one-way hierarchical clustering analysis (Euclidean method, complete linkage) using Z-score for RLE normalized values of all genes expressed in testicular tissue from adult and juvenile cats. C Volcano plot showing differentially expressed genes in adult vs. juvenile samples (Wald test, adjusted p-value < 0.05, absolute fold change > 2). D Principal component analysis plot representing variation in samples from adult and juvenile cats. E Venn diagram of differentially expressed genes from 3 comparison pairs: Adult, Late Juvenile (JF1 and JF2) and Early Juvenile (JF3, JF4 and JF5); Wald test, adjusted p-value < 0.05, absolute fold change > 2. F Heatmap of one-way hierarchical clustering analysis (Euclidean method, complete linkage) using Z-score for RLE normalized values of selected differentially expressed cell markers. Cell markers have been organized to mark more differentiated germ cell populations as they move from the top to the bottom rows. * - differentially expressed in adult vs. late juvenile; ^ - differentially expressed in late vs. early juvenile

Heatmap in Fig. 2 F presents clustering of samples based on germ and somatic cell markers (identified from comprehensive data generated in human testicular tissue [4, 36, 37]) that were differentially expressed in adult vs. juvenile. Genes known to be expressed in more differentiated germ cells (SPO11 to PRM3) had higher expression in adult samples when compared to overall juvenile group, while HOOK1, SPAG6 and CCNA1 were also higher expressed in late compared to early juvenile samples showing continuous increase of these markers (Fig. 2 F, Additional file 4). The expression of spermatogonia markers (UTF1 to KIT) was higher in juvenile compared to adult samples with no differences within the overall juvenile group, while the expression of GFRA1 was also higher in late juvenile compared to adult samples (Fig. 2 F, Additional file 2). Somatic cell markers had higher expression in juvenile compared to adult group, with the expression of Sertoli cell marker AMH continuously decreasing from early to late juvenile to adult (Fig. 2 F, Additional file 4). Results from Fig. 2 F in addition with hierarchical clustering based on all genes (Fig. 2B) allowed us to further separate early juvenile group into early I (JF5) and early II (JF3, JF4).

In sum, early I juvenile age group (JF5) presented testis lacking apparent lamina or lumen with only spermatogonia and Sertoli cells present in tubules (Fig. 2 A, Table S1), high expression of undifferentiated spermatogonia marker UTF1, and low expression of differentiating spermatogonia marker KIT (Fig. 2 F). A tubular structure became progressively apparent in early II juvenile (JF3, JF4) with presence of spermatocytes and only few spermatids in some tubules, while a clear defined lamina and lumen were observed across tubules in late juvenile samples (JF1, JF2) with higher number of spermatids present in more tubules (Fig. 2 A, Table S1), as well as significantly higher expression of spermatocyte markers CCNA1 and SPAG6 compared to early juvenile (Fig. 2 F, Additional file 3). Testicular tissue from adult group showed full establishment of spermatogenesis with many spermatozoa present in tubules (Fig. 2 A, Table S1), increased expression of great number of lncRNA compared to juvenile (Fig. 2 C), as well as higher expression of spermatocyte and spermatid markers (SPO11 to PRM3) and lower expression of spermatogonia (UTF1 to KIT) and somatic cell (VIM to IGF1) markers compared to juvenile age group (Fig. 2 F, Additional file 1).

Changes in functional terms and pathways throughout testicular tissue development in adult and juvenile cats

Figure 3 A visualizes DAVID results of functional analysis in an Enrichment map for DEGs from comparison pair adult vs. juvenile. Out of 4,834 protein coding DEGs, 4,550 were annotated in DAVID database for the domestic cat; Gene Ontology and KEGG Pathways databases were used for gene set enrichment. Figure 3B visualizes DAVID results for DEGs from adult vs. late juvenile and late vs. early juvenile comparison pairs. Additional file 5 contains the web session of both networks with interactive view and data table. Functional terms enriched in testicular tissue from adult animals (upregulated in adult compared to juvenile) were mainly related to sperm motility, processes involved in spermatogenesis and fertilization, and cell cycle and division (Fig. 3 A). Majority of these terms were also enriched in adult samples when comparing to late juvenile and in late compared to early juvenile samples, indicating a continuous increase in function of these terms and pathways with testicular development (Fig. 3B). Functional terms enriched in testicular tissues from juvenile cats (downregulated in adult compared to juvenile) were mainly related to cell signaling, cell adhesion, cell migration, as well as terms associated with extracellular region, cell surface, membrane and response to stimulus (Fig. 3 A). Terms related to immune response were specifically enriched in early juvenile samples when comparing with late juvenile (Fig. 3B). Terms related to cell adhesion, membrane and PI3K-Akt signaling were specifically enriched in late juvenile samples when comparing to adult (Fig. 3B).

Fig. 3
figure 3

Functional terms and pathways enriched in fresh testicular tissues from different development stages. Enrichment map displays enriched gene-sets based on differentially expressed genes in testicular tissues of (A) adult vs. juvenile cats and (B) adult vs. late juvenile and late vs. early juvenile cats. Additional file 5 contains the web session of both networks for interactive view

Effect of vitrification on histology, DNA integrity, and transcriptome in testicular tissues from adult and juvenile cats

Main types of morphological damage caused by vitrification were cellular disorganization, shrinkage from the basal membrane, and nuclear condensation (Fig. 4 A). The percentage of damaged seminiferous tubules was higher in tissues from adult compared to juvenile cats (Fig. 4B). After 24 h of in vitro culture, warmed tissues from juvenile cats retained a higher percentage of morphologically normal seminiferous tubules compared to adult (Fig. 4B). DNA damage could not be attributed to specific cell types and was observed in different parts of seminiferous tubules (Fig. 4 C). Percentage of cells with damaged DNA increased after vitrification of tissues from both adult and juvenile cats (Fig. 4D). After 24 h in vitro culture, warmed tissues from adult cats presented higher percentage of DNA fragmentation compared to juvenile (Fig. 4D).

Fig. 4
figure 4

Effect of vitrification on histology, DNA integrity and transcriptomics in testicular tissues from adult and juvenile domestic cats. A Representation of morphological damage in vitrified testicular tissue sections showing mature tubule with lumen (upper picture) and immature tubule without lumen (lower picture). Black arrowhead: separation of the basal membrane; white arrowhead: nuclear condensation; black arrows: shrinkage from the basal membrane; asterisk: tubular cell disorganization. Scale bar: 20 μm. B Changes in percentage of morphologically normal seminiferous tubules after vitrification and/or 24 h culture. C Representation of DNA fragmentation in testicular tissue sections. White arrowhead: intact DNA (Hoechst); white arrow: fragmented DNA (TUNEL). Scale bar: 20 μm. D Changes in percentage of DNA fragmentation after vitrification and/or 24 h culture. E Heatmap of one-way hierarchical clustering analysis (Euclidean method, complete linkage) using Z-score for RLE normalized values of genes differentially expressed in vitrified vs. fresh samples, independent of age effect (likelihood ratio test, condition effect, adjusted p-value < 0.05). Sample IDs represent age (A, adult; J, juvenile), tissue condition (F, fresh; V, vitrified) and animal (1 to 10). F Enrichment map displaying enriched gene-sets based on genes responding to vitrification (p-value < 0.05, absolute fold change > 1.5) in vitrified vs. fresh samples from adult and juvenile cats. G STRING network of predicted protein interactions for genes responding to vitrification (p-value < 0.05, absolute fold change > 1.5) in vitrified vs. fresh samples from adult and juvenile cats. Protein interaction clusters were analyzed for functional enrichment within STRING app; clusters with no enrichment were removed from visualization. Additional file 7 contains the web session of both Enrichment map and STRING networks for interactive view

Vitrification led to differential expression of 49 genes independent of age (adjusted p-value < 0.05), with 46 downregulated (7 lncRNA, 40 protein coding, 2 pseudogenes) and 3 upregulated (all protein coding) in vitrified vs. fresh tissue (Additional file 6). Hierarchical clustering of samples based on these DEGs separated adult vitrified samples from all the rest representing the biggest vitrification effect in adult tissues (Fig. 4E). One vitrified late juvenile sample (JV1) clustered with vitrified adult, indicating a bigger vitrification effect in that sample compared to the rest juvenile (Fig. 4E). Vitrified and fresh tissues from early juvenile sample(JV5, JF5) clustered together, indicating the smallest vitrification effect in that sample compared to the rest juvenile (Fig. 4E). Clustering of JF5 together with other vitrified tissues may indicate that the transcriptome of this sample is more similar to vitrified rather than fresh tissues (Fig. 4E).

Analysis of vitrification effect separately in adult and juvenile samples identified 21 DEGs in adult and no DEGs in juvenile vitrified vs. fresh tissue (adjusted p-value < 0.05; Additional file 6). To compare the effect of vitrification between adult and juvenile samples, genes satisfying condition of p-value < 0.05 and absolute fold change > 1.5 were selected for functional analysis as genes responsive to vitrification. Out of these selected genes, 18 were shared between the age groups and included 17 downregulated and 1 upregulated gene (Table 1). Out of 116 selected protein coding genes, 110 were annotated in DAVID database for domestic cat; Gene Ontology, KEGG Pathways, InterPro Domains, UniProt Keywords and SMART domains databases were used for gene set enrichment. Figure 4 F visualizes DAVID results in an Enrichment Map network. Figure 4G visualizes networks of predicted protein-protein interactions built using 109 out of 116 selected protein coding genes annotated in STRING database for domestic cat; functional enrichment was performed for each interaction cluster using STRING app, which allowed us to consider pairwise relationships among interacting genes when checking for biological significance. Additional file 7 contains the web session of both networks with interactive view and data table. Functional terms related to transmembrane were enriched in fresh testicular tissue (downregulated in vitrified vs. fresh) for both adult and juvenile ages (Fig. 4 F) and some of the DEGs from that enrichment gene set formed a protein-protein interaction cluster (Fig. 4G). Terms enriched in vitrified tissues (upregulated in vitrified vs. fresh) were related to transcription activity and TNF signaling for adult samples and to nucleus and chromosome in juvenile samples, with shared term of DNA binding enriched for both ages (Fig. 4 F). The biggest protein-protein interaction cluster included both upregulated and downregulated DEGs in vitrified vs. fresh tissue from adult and juvenile samples (Fig. 4G). Genes with the most predicted protein-protein interactions were FOS, JUNB and CREB1 (8 interactions each), JUN and EGR1 (6 interactions each), and FOSB and NFAT5 (5 interactions each).

Table 1 Genes responsive to vitrification in both adult and juvenile testicular tissues

Discussion

This is the first overview of transcriptome dynamics coupled with structural analysis in testicular tissue from adult and juvenile domestic cats. Transcriptomic findings enabled to further characterize juvenile period, distinguishing between early and late juvenile tissues. Changes in gene expression and functional pathways were extensive from early to late juvenile to adult development stages. Additionally, tissues from juvenile animals were more resilient to vitrification compared to adult counterparts, with early juvenile sample responding the least to vitrification and late juvenile sample response being closest to adult tissues.

Gradual testis maturation during puberty in cat

In the domestic cat, the spermatogenic function is suggested to become mature at 8 to 10 months of age with initial activation and first signs of spermatogenesis at 5 to 6 months [8, 9]. In our study, we defined juvenile period as the development period after the infant stage but prior to full adult sexual maturity to make it comparable with the recent scRNA-seq study on testis development during human puberty [4]. In juvenile samples, various degrees of maturation and germ cell differentiation were observed, from tubules with only spermatogonia in early I juvenile stage to the presence of some spermatocytes in early II juvenile to many spermatocytes and round spermatids with occasional spermatozoa in late juvenile. The expression of cellular markers associated with more differentiated germ cells also increased from early to late juvenile period, while the expression of markers associated with undifferentiated spermatogonia was highest in early juvenile period. Functional terms and pathways related to spermatogenesis and sperm motility were enriched already in late juvenile and stayed enriched in adult samples. All of this supports the gradual progressive establishment of spermatogenesis in the domestic cat beginning in late juvenile period. This is similar to humans, where low levels of incomplete spermatogenesis are observed in portions of the testis in juveniles prior to puberty, and the appearance of first spermatozoa during juvenile period does not mark completion of spermatogenic development, but the beginning of the final stages of puberty, where gradual improvement in spermatogenic efficiency is the mechanism that leads to complete maturation [38]. It is suggested that rather than sudden activation of the testis at puberty there is a slow and progressive increase of activity from mid-childhood or even earlier [38]. This also is aligned with recent scRNA-seq studies, showing that spermatogenic cell phenotypes exist on a continuum rather than in distinct subgroups separated by large transcriptome changes [39].

Sertoli cell maturation

The present study reports downregulation of functional terms and pathways related to immune response in late compared to early juvenile period, which also coincided with the emergence of more spermatocytes and spermatids, as well as the decrease in anti-Mullerian hormone (AMH) expression. During the Sertoli cell differentiation in puberty, neighboring Sertoli cells form tight junctions, which contribute to the blood-testis barrier that permits the establishment of a special microenvironment needed for spermatogenesis [40, 41]. This barrier, together with the expression by Sertoli cells of immunoregulatory factors that actively suppress innate, humoral and cell-mediated immune responses, makes the whole testis immune privileged [42]. Around the same time, the expression of AMH in Sertoli cells decreases, which is complemented by meiotic entry of germ cells [1]. Based on this, we hypothesize that maturation of Sertoli cells begins in late juvenile period in the domestic cat. The recent scRNA-seq study in human developing testis showed that the maturation of Sertoli cells in different tubules is asynchronous and proceeds gradually, with first mature Sertoli cells starting to emerge in samples from 11 years old onward [4].

Long non-coding RNAs in developing testis

Testis is a tissue with the most tissue-specific genes by far [17], higher fraction of genes and more diverse mRNA (less dominance of a few highly expressed genes) [43]. It is also a tissue with the highest amount of long non-coding RNAs (lncRNAs) even when comparing to brain or liver [44, 45]. lncRNAs are a large class of noncoding RNAs more than 200 nucleotides in length. Unlike mRNAs, lncRNAs exhibit unique cellular localization patterns highly correlated with the functions they perform in the cell [46, 47]. Recent studies revealed critical roles of lncRNAs during spermatogenesis [48] and accumulation of lncRNAs during meiotic and post-meiotic stages of spermatogenesis [49, 50]. The present study reports a progressive increase in expression of lncRNAs from early to late juvenile period to adult in cat testicular tissue, which may indicate on the role of lncRNAs in domestic cat spermatogenesis as well. We could not conduct functional analysis on the set of differentially expressed lncRNAs because there currently is no comprehensive functional database for domestic cat lncRNAs and lncRNAs have low sequence conservation. However, this list will inform future studies for further analysis of lncRNA functions and their interactions with mRNA and other non-coding RNAs.

Juvenile testes are more resilient to vitrification protocol

Our results show that vitrification led to more detrimental changes in mature adult tissues compared to immature juvenile which only progressed after 24 h culture and were expressed in increased percentage of damaged seminiferous tubules and DNA fragmentation. Mature testicular tissue from adult cats also was more responsive transcriptionally to vitrification compared to juvenile. Similarly, immature testes from domestic cat were reported to be more resilient to microwave-assisted dehydration [25] and more successful in surviving and establishing spermatogenesis in xenografts [51, 52], compared to mature testes from adult cats. Specifically, we observed the downregulation of functional terms related to focal and cell-cell adhesion in adult compared to late juvenile, as well as to the whole juvenile group. Thus, the detrimental effect of vitrification and 24 h culture on testicular tissue in adult might be related to its increased sensitivity to mechanical damage due to decreased connections between cells and cells and extracellular matrix. At the same time, one sample from late juvenile period that exhibited the highest number of meiotic and post-meiotic germ cells out of all juvenile group had a similar transcriptomic response to adult samples, which might indicate the effect of the presence of these cells on tissue sensitivity to vitrification stress.

Initial response of testicular tissue to vitrification and warming

We identified a cluster of genes responsive to vitrification in testicular tissue of adult and juvenile cats and forming predicted protein-protein interactions and enriched in transcription and TNF signaling pathway. These genes were upregulated in response to vitrification and included early response transcription factors FOS, FOSB, JUN, JUNB, and EGR1. Interestingly, the same list of genes has been reported in studies analyzing the response of different tissues to warm ischemia [53,54,55,56]. In those studies, normal and/or tumor tissues that underwent warm ischemia due to delayed processing after surgical incision demonstrated increased expression of FOS, FOSB, JUN, EGR1 as well as upregulation in immune system pathways compared to tissues processes immediately [53,54,55,56]. We hypothesize that the damage occurring in testicular tissue during vitrification is caused by the warming step and may be similar to the response of tissues to warm ischemia. Interestingly, in one study that looked at short intervals of warm ischemia, the pattern of gene and protein expression in the tissue changed within minutes following surgical excision [57]. In our protocol the initial 5-sec warming step is followed by cryoprotectant removal at room temperature for 15 min, which might be the period when testicular tissue undergoes stress similar to warm ischemia resulting in damage.

Study limitations and next steps

In our current study, we performed bulk RNA-seq on a whole testicular tissue to have the first look on transcriptome dynamics in the developing testis of domestic cat. The addition of cell-level transcriptome information would have been invaluable to our study, however, there has been no scRNA-seq performed so far on any of the cat tissues and, therefore, no protocols established yet. We recognize that the bulk analysis of mature RNA limits our interpretations. For example, downregulation of somatic cell markers in adult compared to juvenile samples is most likely related to the decrease in somatic to germ cell ratio (in mammals, mature testes germ cells make up to 90 % of the tubular mass, compared to less than 5–10 % in immature) [2]. Our study, however, is the first necessary step in unraveling transcriptome dynamics of the developing testis in the domestic cat and the effect of age on tissue response to preservation protocols. Future studies focusing on cell-specific expression and post-transcriptional mechanisms, as well as linking transcriptomic and proteomics [58] would be essential to understand the complex process of puberty in cats, as well as tissue resilience.

Due to unavailability of domestic cat samples from the infant stage, we could only analyze testis development starting from early juvenile period. In future, adding neonatal tissues to the analysis, as well as including more samples from juvenile period, would contribute greatly to our understanding of the full timeline of testis development in domestic cat.

Conclusions

Our study generated a high-quality transcriptomic data for testicular tissue of adult and juvenile domestic cats, which provides an important resource for future studies on testis development, spermatogenesis and fertility in cats, as well as new insights into tissue resilience. Transcriptomic data can also contribute to identification of druggable protein targets in male reproductive tracts and development of male contraception [59]. Our study provides a network of functional terms and pathways that are affected by age in the domestic cat and are either enriched in adult or juvenile testicular tissues. The interactive view of our network allows navigation through many enriched terms, their interconnection, and the associated set of genes from our study. This is the first study providing data on transcriptomic dynamic coupled with structural analysis in the cat testis according to the age and exposure to cryopreservation. Collective findings will also enable the optimization of testicular tissue preservation.

Materials and methods

Collection of testicular tissues

Testes from juvenile (3 to 6 months old) and adult (over 1 year old) male domestic cats were collected on different days after routine orchiectomy at local veterinary clinics and transported in phosphate buffered saline (PBS) at 4 °C to the laboratory within 6 h of excision. Testes were washed once with PBS, dissected from surrounding tissues and cut in pieces of approximately 2–3 mm3 in handling medium composed of Hepes-Ham’s F10 medium (Irvine Scientific, Santa Ana, CA) supplemented with 1mM pyruvate, 2 mM L-glutamine, 100 IU/mL penicillin, 100 µg/mL streptomycin, 2.5 % fetal bovine serum (FBS). For each animal, tissue pieces were either (1) fixed overnight in Bouin’s solution (histomorphology) or 4 % paraformaldehyde in PBS (TUNEL analysis), embedded in paraffin and sectioned at a thickness of 5 μm, (2) incubated overnight in RNAlatertm solution (Invitrogen, Carlsbad, CA) for RNA isolation, or (3) processed for vitrification and/or in vitro culture as described below (Fig. 1). All chemicals and reagents were purchased from Sigma-Aldrich (St. Louis, MO), unless otherwise indicated.

Vitrification and warming

Vitrification and warming was performed using modified protocol reported previously for domestic cat testicular tissue [26]. Tissue biopsies were exposed to equilibration solution (1.4 M dimethyl sulfoxide (DMSO) + 1.4 M glycerol + 0.25 M sucrose in Ham’s F10) for 10 min at room temperature (~ 22 °C) followed by a vitrification solution (2.8 M dimethyl sulfoxide (DMSO) + 2.8 M glycerol + 0.5 M sucrose + 10 % FBS in Ham’s F10) for 5 min at room temperature (~ 22 °C), then placed in cryotubes and plunged directly into liquid nitrogen and stored for at least one week in liquid nitrogen. Warming was performed by immersing cryotubes in a water bath at 50 °C for 5 s. Tissue fragments were then transferred to a sucrose gradient (0.50 M; 0.25 M; 0.00 M in Hepes-Ham’s F10 and 20 % of FBS) for 5 min at each step at room temperature (~ 22 °C) in order to remove the cryoprotectants. Warmed tissues were fixed for histology and RNA isoltation as described above or processed for in vitro culture.

In vitro culture

Tissue culture was performed using the same protocol reported previously for domestic cat [26]. Tissue fragments were placed into 1cm3 pieces of 1.5 % agarose gel that were pre-conditioned by immersion in culture medium composed of Hepes-Ham`s F10 (supplemented with 2mM L-glutamine, 1mM pyruvate, 100 IU/ml penicillin, 100 µg/ml streptomycin and 5 % FBS). Two tissue biopsies on each gel were incubated for 24 h in a 4-well culture plate with 400 µl of culture medium at 38.5 °C in a humidified atmosphere of 5 % CO2 in air.

Assessment of tissue histomorphology

Testicular tissue morphology was assessed via hematoxylin-eosin staining. Samples used for RNA-seq were scored using Johnsen method for registration of spermatogenesis [35]. In short, each seminiferous tubule in testicular tissue section was given a score from 10 to 1 according to the presence or absence of the main cell types arranged in the order of maturity: presence of spermatozoa scores 10, 9 or 8; spermatids (and no further) 7 or 6; spermatocytes (and no further) 5 or 4; only spermatogonia 3, only Sertoli cells 2 and no cells 1. The mean score for each sample was calculated by multiplying the number of tubuli recorded at each score with the score and then dividing the sum of all 10 multiplications by the total number of tubuli recorded [35].

Integrity of seminiferous tubules and cells was evaluated according to criteria previously established [26]. Intact tubules with no detachment of cells from the basement membrane, no rupture of stroma, no swelling of the lamina propria and normal junctions between cells were considered as a normal structure (score 1). Score 0 was attributed to tubules with changes in any of the previous criteria. A total of 30 randomly selected seminiferous tubules for each animal in each experimental group were classified as normal structure (score 1) or damage structure (score 0) totaling 180 tubules per group. Percentage of normal seminiferous tubules was calculated relative to the total number of observed tubules.

Assessment of DNA integrity

DNA integrity was assessed using the protocol reported previously for domestic cat testicular tissue using the In-Situ Cell Death Detection kit (Roche, Basel, Switzerland) [25]. Deparaffinized rehydrated sections were rinsed twice with 0.05 % Triton X-100 in PBS for 5 min each, permeabilized with 0.5 % Triton X-100 in PBS for 30 min, rinsed once with 0.05 % Triton X-100 in PBS for 5 min and incubated in TUNEL reaction mixture (enzyme solution with terminal deoxynucleotidyl transferase (TdT) and label solution with nucleotide polymers) for 1 h at 37 °C within a humidified darkened container. Negative control (omitted TdT) was included in each run. For positive control, tissue was incubated with recombinant DNase I for 10 min before labeling. The nucleus of all cells was stained with Hoechst 33,342 (1:100) in a humidified chamber for 10 min at room temperature and then, the slides were mounted with Vectashield mounting medium (Vector laboratories, Burlingame, CA). We evaluated 60 images per experimental group, which were captured using an Olympus BX41 epifluorescence microscope (Olympus Corporation, Tokyo Japan) with SPOT advanced software 5.0 (Diagnostic Instruments, Inc., Sterling Heights, MI). Percentage of DNA damage (TUNEL positive cells) was calculated relative to the total number of observed cells.

RNA preparation

Twenty samples were selected for transcriptomic analysis and assigned to the following 4 groups: adult, fresh tissue (AF, n = 5), adult, vitrified/warmed tissue (AV, n = 5), juvenile, fresh tissue (JF, n = 5) and juvenile, vitrified/warmed tissue (JV, n = 5; Fig. 1). Total RNA was isolated from up to 10 mg of tissue using PureLink RNA Mini Kit with on-column DNase Set (Invitrogen); tissue was homogenized in RNA lysis buffer using TissueLyser (Qiagen, Hilden, Germany; 2 × 2 min at 30 Hz; 5 mm stainless steel beads). Concentration and purity of isolated RNA was measured with NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA); RNA integrity was assessed using 2100 Bioanalyzer instrument (Agilent Technologies, Santa Clara, CA). Purified RNA was stored in nuclease-free water at -80 °C until library preparation.

Library preparation and transcriptome sequencing

Only samples with RIN ≥ 7 were used for library preparation. Sequencing libraries were generated using TruSeq Stranded mRNA LT Sample Prep Kit from Illumina (San Diego, CA) according to the manufacturer’s recommendations. In short, the workflow included randomly fragmenting total RNA for short read sequencing, reverse transcribing fragmented RNA into cDNA, ligating adaptors onto both ends of the cDNA fragments, amplifying cDNA and selecting fragments with insert sizes between 200 and 400 bp. The libraries were sequenced 150 bp paired-end using an Illumina NovaSeq 6000 System at the Psomagen Inc. (formerly Macrogen Corp., Rockville, MD) attempting 30 million reads depth per sample and 5 biological replicates per group.

Quality control

The quality of produced data was determined by the phred quality score at each cycle using FastQC (v. 0.11.7). Trimmomatic (v. 0.38) [60] program was used to remove adapter sequences and bases with base quality lower than three from the ends. Using sliding window method, bases of reads that did not qualify for window size 4 and mean quality 15 were trimmed. Afterwards, reads with length shorter than 36 bp were dropped to produce trimmed data. Quality information for each sample after trimming is provided in Table S2.

Reads mapping and gene expression levels quantification

Trimmed reads were mapped to reference genome GCF_000181335.3_Felis_catus_9.0 with HISAT2 (v. 2.1.0) [61], which is known to handle spliced read mapping through Bowtie2 (v. 2.3.4.1) aligner, splice-aware aligner. Table S3 shows the statistic obtained from HISAT2. After the read mapping, known genes and transcripts were assembled with StringTie (v. 1.3.4d) [62, 63] based on reference genome model. After assembly, the abundance of gene/transcript was calculated in the read count and normalized value as FPKM (Fragment per Kilobase of transcript per Million mapped reads) for a sample.

Differential expression analysis

Differentially expression analysis was performed on data obtained from 20 samples using the DESeq2 R package [64]. For visualization, size factors were estimated from the count data and the Relative Log Expression (RLE) normalization was used to obtain regularized log transformed values. These normalized values were then used for principal component analysis (plotPCA function in DESeq2 R package) and creation of clustered heatmaps (pheatmap R package). Wald test (age effect) and likelihood ratio test (vitrification effect) were used on genes that passed an independent filtering step and resulting P values were adjusted for multiple testing using the Benjamini-Hochberg procedure.

Functional enrichment analysis using DAVID and Enrichment Map

Selected genes from differential expression analysis were used for gene-set functional enrichment analysis with DAVID tool [65], setting species to domestic cat. For each comparison pair, total number of genes and separately up- and downregulated genes were analyzed. EASE score (modified Fisher Exact p-value of enrichment) was set to 0.05. Functional enrichment network was built based on DAVID output charts of gene-set enrichment for each comparison pair using Enrichment Map app (v. 3.3.0) [66] in Cytoscape software (v. 3.8.0) [67, 68] with Overlap set to 0.5. Autoannotate App (v. 1.3.3) with MLC algorithm based on similarity coefficient was used to create annotated groups.

In silico protein-protein interaction analysis using STRING

In silico protein-protein interaction analysis of selected genes was performed on the basis of the STRING database for the domestic cat [69, 70]. Interaction network was built based on the list of selected genes from each comparison pair using stringApp (v. 1.5.1, [71]) in Cytoscape with confidence cutoff score set to 0.4. Functional enrichment of formed clusters was performed using domestic can genome as a background, enriched terms were analyzed with varying redundancy cutoff settings.

Statistical analysis

Statistical analysis of RNA-seq data is provided above. For histomorphology and TUNEL analyses, data were expressed as mean and standard error and analyzed using the statistical software graphpad prism version 5.01 (GraphPad Software Inc., San Diego, CA). Data was tested for normality (Shapiro–Wilk test) and homoscedasticity (Levene’s test). Analysis of variance (ANOVA) followed by Tukey test was used to compare the effect of vitrification and 24 h in vitro culture in different experiment groups.

Availability of data and materials

All sequence data was generated in the current study and deposited to the NCBI SRA repository, BioProject (PRJNA662384). The differential gene expression datasets supporting the conclusions of this article were generated in the current study and are included within the article as additional files.

References

  1. Koskenniemi JJ, Virtanen HE, Toppari J. Testicular growth and development in puberty. Curr Opin Endocrinol Diabetes Obes. 2017;24:215–24.

    Article  CAS  PubMed  Google Scholar 

  2. Schlatt S, Ehmcke J. Regulation of spermatogenesis: An evolutionary biologist’s perspective. Seminars in Cell & Developmental Biology. 2014;29:2–16.

    Article  CAS  Google Scholar 

  3. Tharmalingam MD, Jorgensen A, Mitchell RT. Experimental models of testicular development and function using human tissue and cells. Molecular and Cellular Endocrinology. 2018;468:95–110.

    Article  CAS  PubMed  Google Scholar 

  4. Guo J, Nie X, Giebler M, Mlcochova H, Wang Y, Grow EJ, et al. The Dynamic Transcriptional Cell Atlas of Testis Development during Human Puberty. Cell Stem Cell. 2020;26:262–276.e4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Grive KJ, Hu Y, Shu E, Grimson A, Elemento O, Grenier JK, et al. Dynamic transcriptome profiles within spermatogonial and spermatocyte populations during postnatal testis maturation revealed by single-cell sequencing. PLOS Genetics. 2019;15:e1007810.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Comizzoli P, Paulson EE, McGinnis LK. The mutual benefits of research in wild animal species and human-assisted reproduction. J Assist Reprod Genet. 2018;35:551–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Comizzoli P. Biobanking efforts and new advances in male fertility preservation for rare and endangered species. Asian J Androl. 2015;17:640–5.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tsutsui T, Kuwabara S, Kuwabara K, Kugota Y, Kinjo T, Hori T. Development of spermatogenic function in the sex maturation process in male cats. J Vet Med Sci. 2004;66:1125–7.

    Article  CAS  PubMed  Google Scholar 

  9. Sánchez B, Pizarro M, García P, Flores J. Postnatal development of seminiferous tubules in the cat. J Reprod Fertil Suppl. 1993;47:343–8.

    PubMed  Google Scholar 

  10. D Francisco F, Lapuente C, Merlo ML, Barbeito C, Gobello C. Feline testicular ultrasonogram differentiates pre vs. postpubertal and normal vs. disrupted spermatogenesis. Theriogenology. 2020;157:503–7.

    Article  CAS  PubMed  Google Scholar 

  11. França LR, Godinho CL. Testis morphometry, seminiferous epithelium cycle length, and daily sperm production in domestic cats (Felis catus). Biol Reprod. 2003;68:1554–61.

    Article  PubMed  Google Scholar 

  12. Braun BC, Okuyama MW, Müller K, Dehnhard M, Jewgenow K. Steroidogenic enzymes, their products and sex steroid receptors during testis development and spermatogenesis in the domestic cat (Felis catus). The Journal of Steroid Biochemistry and Molecular Biology. 2018;178:135–49.

    Article  CAS  PubMed  Google Scholar 

  13. Braun BC, Müller K, Jewgenow K. Expression profiles of relaxin family peptides and their receptors indicate their influence on spermatogenesis in the domestic cat (Felis catus). Domest Anim Endocrinol. 2015;52:25–34.

    Article  CAS  PubMed  Google Scholar 

  14. Li Y, Li J, Fang C, Shi L, Tan J, Xiong Y, et al. Genome-wide differential expression of genes and small RNAs in testis of two different porcine breeds and at two different ages. Sci Rep. 2016;6:26852.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Li B, He X, Zhao Y, Bai D, Du M, Song L, et al. Transcriptome profiling of developing testes and spermatogenesis in the Mongolian horse. BMC Genetics. 2020;21:46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang G, Li Y, Yang Q, Xu S, Ma S, Yan R, et al. Gene expression dynamics during the gonocyte to spermatogonia transition and spermatogenesis in the domestic yak. Journal of Animal Science and Biotechnology. 2019;10:64.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Djureinovic D, Fagerberg L, Hallström B, Danielsson A, Lindskog C, Uhlén M, et al. The human testis-specific proteome defined by transcriptomics and antibody-based profiling. Mol Hum Reprod. 2014;20:476–88.

    Article  CAS  PubMed  Google Scholar 

  18. Laiho A, Kotaja N, Gyenesei A, Sironen A. Transcriptome Profiling of the Murine Testis during the First Wave of Spermatogenesis. PLOS ONE. 2013;8:e61558.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Davis BW, Seabury CM, Brashear WA, Li G, Roelke-Parker M, Murphy WJ. Mechanisms Underlying Mammalian Hybrid Sterility in Two Feline Interspecies Models. Molecular Biology and Evolution. 2015;32:2534–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Silva AM da, Pereira AF, Comizzoli P, Silva AR. Cryopreservation and Culture of Testicular Tissues: An Essential Tool for Biodiversity Preservation. Biopreservation and Biobanking. 2020;18:235–43.

  21. Onofre J, Baert Y, Faes K, Goossens E. Cryopreservation of testicular tissue or testicular cell suspensions: a pivotal step in fertility preservation. Hum Reprod Update. 2016;22:744–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Picton HM, Wyns C, Anderson RA, Goossens E, Jahnukainen K, Kliesch S, et al. A European perspective on testicular tissue cryopreservation for fertility preservation in prepubertal and adolescent boys. Hum Reprod. 2015;30:2463–75.

    Article  PubMed  Google Scholar 

  23. Dohle GR. Male infertility in cancer patients: Review of the literature. International Journal of Urology. 2010;17:327–31.

    Article  PubMed  Google Scholar 

  24. Valli-Pulaski H, Peters KA, Gassei K, Steimer SR, Sukhwani M, Hermann BP, et al. Testicular tissue cryopreservation: 8 years of experience from a coordinated network of academic centers. Hum Reprod. 2019;34:966–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Silva HVR, da Silva AM, Lee P-C, Brito BF, Silva AR, da Silva LDM, et al. Influence of Microwave-Assisted Drying on Structural Integrity and Viability of Testicular Tissues from Adult and Prepubertal Domestic Cats. Biopreservation and Biobanking. 2020;18:415–24.

    Article  CAS  PubMed  Google Scholar 

  26. Lima DBC, Silva LDM da, Comizzoli P. Influence of warming and reanimation conditions on seminiferous tubule morphology, mitochondrial activity, and cell composition of vitrified testicular tissues in the domestic cat model. PLOS ONE. 2018;13:e0207317.

  27. Lima DBC, Silva TFP da, Aquino-Cortez A, Leiva-Revilla J, Silva LDM da. Vitrification of testicular tissue from prepubertal cats in cryotubes using different cryoprotectant associations. Theriogenology. 2018;110:110–5.

  28. Macente BI, Apparicio M, Mansano CFM, Tavares MR, Fonseca-Alves CE, Sousa BP, et al. Effect of cryopreservation on sperm DNA fragmentation and apoptosis rates in the testicular tissue of domestic cats. Anim Reprod Sci. 2019;211:106224.

    Article  CAS  PubMed  Google Scholar 

  29. Peris-Frau P, Soler AJ, Iniesta-Cuerda M, Martín-Maestro A, Sánchez-Ajofrín I, Medina-Chávez DA, et al. Sperm Cryodamage in Ruminants: Understanding the Molecular Changes Induced by the Cryopreservation Process to Optimize Sperm Quality. International Journal of Molecular Sciences. 2020;21:2781.

    Article  CAS  PubMed Central  Google Scholar 

  30. Fraser L, Brym P, Pareek CS, Mogielnicka-Brzozowska M, Paukszto Ł, Jastrzębski JP, et al. Transcriptome analysis of boar spermatozoa with different freezability using RNA-SEq. Theriogenology. 2020;142:400–13.

    Article  CAS  PubMed  Google Scholar 

  31. Dai D-H, Qazi IH, Ran M-X, Liang K, Zhang Y, Zhang M, et al. Exploration of miRNA and mRNA Profiles in Fresh and Frozen-Thawed Boar Sperm by Transcriptome and Small RNA Sequencing. International Journal of Molecular Sciences. 2019;20:802.

    Article  CAS  PubMed Central  Google Scholar 

  32. Ran M-X, Li Y, Zhang Y, Liang K, Ren Y-N, Zhang M, et al. Transcriptome Sequencing Reveals the Differentially Expressed lncRNAs and mRNAs Involved in Cryoinjuries in Frozen-Thawed Giant Panda (Ailuropoda melanoleuca) Sperm. International Journal of Molecular Sciences. 2018;19:3066.

    Article  PubMed Central  Google Scholar 

  33. Karlsson JOM, Toner M. Long-term storage of tissues by cryopreservation: critical issues. Biomaterials. 1996;17:243–56.

    Article  CAS  PubMed  Google Scholar 

  34. Amelkina O, Comizzoli P. Initial response of ovarian tissue transcriptome to vitrification or microwave-assisted dehydration in the domestic cat model. BMC Genomics. 2020;21:828.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Johnsen SG. Testicular biopsy score count–a method for registration of spermatogenesis in human testes: normal values and results in 335 hypogonadal males. Hormones. 1970;1:2–25.

    CAS  PubMed  Google Scholar 

  36. Guo J, Grow EJ, Mlcochova H, Maher GJ, Lindskog C, Nie X, et al. The adult human testis transcriptional cell atlas. Cell Research. 2018;28:1141–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Pontén F, Schwenk JM, Asplund A, Edqvist P-HD. The Human Protein Atlas as a proteomic resource for biomarker discovery. Journal of Internal Medicine. 2011;270:428–46.

    Article  PubMed  Google Scholar 

  38. Chemes HE. Infancy is not a quiescent period of testicular development. International Journal of Andrology. 2001;24:2–7.

    Article  CAS  PubMed  Google Scholar 

  39. Suzuki S, Diaz VD, Hermann BP. What has single-cell RNA-seq taught us about mammalian spermatogenesis? Biology of Reproduction. 2019;101:617–34.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Griswold MD. 50 years of spermatogenesis: Sertoli cells and their interactions with germ cells. Biology of Reproduction. 2018;99:87–100.

    Article  PubMed  PubMed Central  Google Scholar 

  41. França LR, Hess RA, Dufour JM, Hofmann MC, Griswold MD. The Sertoli cell: one hundred fifty years of beauty and plasticity. Andrology. 2016;4:189–212.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Mital P, Hinton BT, Dufour JM. The Blood-Testis and Blood-Epididymis Barriers Are More than Just Their Tight Junctions1. Biology of Reproduction. 2011;84:851–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ramsköld D, Wang ET, Burge CB, Sandberg R. An Abundance of Ubiquitously Expressed Genes Revealed by Tissue Transcriptome Sequence Data. PLOS Computational Biology. 2009;5:e1000598.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Soumillon M, Necsulea A, Weier M, Brawand D, Zhang X, Gu H, et al. Cellular source and mechanisms of high transcriptome complexity in the mammalian testis. Cell Rep. 2013;3:2179–90.

    Article  CAS  PubMed  Google Scholar 

  45. Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, et al. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011;25:1915–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Chen L-L. Linking Long Noncoding RNA Localization and Function. Trends Biochem Sci. 2016;41:761–72.

    Article  CAS  PubMed  Google Scholar 

  47. Ulitsky I. Evolution to the rescue: using comparative genomics to understand long non-coding RNAs. Nature Reviews Genetics. 2016;17:601–14.

    Article  CAS  PubMed  Google Scholar 

  48. Wichman L, Somasundaram S, Breindel C, Valerio DM, McCarrey JR, Hodges CA, et al. Dynamic expression of long noncoding RNAs reveals their potential roles in spermatogenesis and fertility†. Biology of Reproduction. 2017;97:313–23.

    Article  PubMed  Google Scholar 

  49. Chalmel F, Lardenois A, Evrard B, Rolland AD, Sallou O, Dumargne M-C, et al. High-Resolution Profiling of Novel Transcribed Regions During Rat Spermatogenesis1. Biology of Reproduction. 2014;91. doi:https://doi.org/10.1095/biolreprod.114.118166.

  50. Rolland AD, Evrard B, Darde TA, Le Béguec C, Le Bras Y, Bensalah K, et al. RNA profiling of human testicular cells identifies syntenic lncRNAs associated with spermatogenesis. Human Reproduction. 2019;34:1278–90.

    Article  CAS  PubMed  Google Scholar 

  51. Mota PC, Ehmcke J, Westernströer B, Gassei K, Ramalho-Santos J, Schlatt S. Effects of different storage protocols on cat testis tissue potential for xenografting and recovery of spermatogenesis. Theriogenology. 2012;77:299–310.

    Article  PubMed  Google Scholar 

  52. Kim Y, Selvaraj V, Pukazhenthi B, Travis AJ, Kim Y, Selvaraj V, et al. Effect of donor age on success of spermatogenesis in feline testis xenografts. Reprod Fertil Dev. 2007;19:869–76.

    Article  PubMed  Google Scholar 

  53. Pedersen IS, Thomassen M, Tan Q, Kruse T, Thorlacius-Ussing O, Garne JP, et al. Differential effect of surgical manipulation on gene expression in normal breast tissue and breast tumor tissue. Molecular Medicine. 2018;24:57.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Pearce DA, Arthur LM, Turnbull AK, Renshaw L, Sabine VS, Thomas JS, et al. Tumour sampling method can significantly influence gene expression profiles derived from neoadjuvant window studies. Scientific Reports. 2016;6:29434.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Almeida A, Paul Thiery J, Magdelénat H, Radvanyi F. Gene expression analysis by real-time reverse transcription polymerase chain reaction: influence of tissue handling. Analytical Biochemistry. 2004;328:101–8.

    Article  CAS  PubMed  Google Scholar 

  56. Huang J, Qi R, Quackenbush J, Dauway E, Lazaridis E, Yeatman T. Effects of Ischemia on Gene Expression. Journal of Surgical Research. 2001;99:222–7.

    Article  CAS  Google Scholar 

  57. Spruessel A, Steimann G, Jung M, Lee SA, Carr T, Fentz A-K, et al. Tissue ischemia time affects gene and protein expression patterns within minutes following surgical tumor excision. BioTechniques. 2004;36:1030–7.

    Article  CAS  PubMed  Google Scholar 

  58. Chalmel F, Rolland AD. Linking transcriptomics and proteomics in spermatogenesis. Reproduction. 2015;150:R149–57.

    Article  CAS  PubMed  Google Scholar 

  59. Robertson MJ, Kent K, Tharp N, Nozawa K, Dean L, Mathew M, et al. Large-scale discovery of male reproductive tract-specific genes through analysis of RNA-seq datasets. BMC Biol. 2020;18:103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nature Methods. 2015;12:357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Pertea M, Pertea GM, Antonescu CM, Chang T-C, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nature Biotechnology. 2015;33:290–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016;11:1650–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Protocols. 2009;4:44–57.

    Article  CAS  Google Scholar 

  66. Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0013984. Accessed 1 May 2021.

  67. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003;13:2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, et al. Integration of biological networks and gene expression data using Cytoscape. Nature Protocols. 2007;2:2366–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Research. 2015;43:D447–52.

    Article  CAS  PubMed  Google Scholar 

  70. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Research. 2017;45:D362–8.

    Article  CAS  PubMed  Google Scholar 

  71. Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J Proteome Res. 2019;18:623–32.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr Brent Whitaker (Animal Rescue Inc.) and Dr Keiko Antoku, and their staff for providing domestic cat testes. We thank Dr. Robert Fleischer and Nancy McInerney from Center for Conservation Genomics (Smithsonian Institution) for providing laboratory space and equipment for RNA-seq experiment. Smithsonian Institution High Performance Cluster (SI/HPC) was used for sequence data analysis.

Funding

This research was supported by Smithsonian Scholarly Studies Awards Program, CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazilian Federal Agency for the Support and Evaluation of Graduate Education), and PDSE (Programa de Doutorado Sanduíche no Exterior— internship program PhD abroad).

Author information

Authors and Affiliations

Authors

Contributions

OA and AMS designed the study, collected, analyzed and interpreted the data, and wrote the whole manuscript; ARS participated in the data analysis and writing of the manuscript; PC designed the study, interpreted the data, and wrote the whole manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Pierre Comizzoli.

Ethics declarations

Ethics approval and consent to participate

The study did not require the approval of the Animal Care and Use Committee of the Smithsonian Conservation Biology Institute because cat ovaries were collected at local veterinary clinics as byproducts from owner-requested routine ovario-hysterectomies.

Consent for publication

Not applicable.

Competing interests

The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific research.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

Differentially expressed genes in fresh testicular tissue of adult compared to juvenile cat; Wald test, adjusted p-value < 0.05, absolute fold change > 2.

Additional file 2

Differentially expressed genes in fresh testicular tissue of adult compared to late juvenile cat; Wald test, adjusted p-value < 0.05, absolute fold change > 2.

Additional file 3

Differentially expressed genes in fresh testicular tissue of late compared to early juvenile cat; Wald test, adjusted p-value < 0.05, absolute fold change > 2.

Additional file 4

Differentially expressed genes in fresh tissue shared between all comparison pairs adult vs. late juvenile, adult vs. early juvenile and late vs. early juvenile (Sheet-1) or only between comparison pairs adult vs. late juvenile and late vs. early juvenile (Sheet-2); Wald test, adjusted p-value < 0.05, absolute fold change > 2.

Additional file 5

Web session interactive view with data table of Enrichment Map networks for comparison pairs adult vs. juvenile, adult vs. late juvenile and late vs. early juvenile. For description and legends refer to Fig. 3 A (adult vs. juvenile) and Fig. 3B (adult vs. late juvenile, late vs. early juvenile).

Additional file 6

Genes affected by vitrification independent of age; likelihood ratio test, condition effect, adjusted p-value < 0.05 (Sheet-1). Genes responsive to vitrification protocol in adult (Sheet-2) and juvenile (Sheet-3) samples; likelihood ratio test, p-value < 0.05, absolute fold change > 1.5.

Additional file 7

Web session interactive view with data table for Enrichment map and STRING network for genes responsive to vitrification in adult and juvenile samples. For description and legend refer to Fig. 4 F (Enrichment map) and Fig. 4G (STRING network).

Additional file 8

Table S1. Johnsen score counts for testicular samples used in RNA-seq analysis.

Additional file 9

Table S2. Trimmed data statistics.

Additional file 10

Table S3. Mapped data statistics.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Amelkina, O., Silva, A.M.d., Silva, A.R. et al. Transcriptome dynamics in developing testes of domestic cats and impact of age on tissue resilience to cryopreservation. BMC Genomics 22, 847 (2021). https://doi.org/10.1186/s12864-021-08099-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08099-8

Keywords