[go: up one dir, main page]
Next Article in Journal
Isolation of Aerobic Bacterial Species Associated with Palpable Udder Defects in Non-Dairy Ewes
Previous Article in Journal
Reconstruction of the Quadriceps Extensor Mechanism with a Calcaneal Tendon–Bone Allograft in a Dog with a Resorbed Tibial Tuberosity Fracture
Previous Article in Special Issue
Specificity of Key Sex Determination Genes in a Mammal with Ovotestes: The European Mole Talpa europaea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 14
Issue 16
10.3390/ani14162316
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria

by
Yukako Katsura
1,*,
Shuji Shigenobu
2 and
Yoko Satta
3
1
Center for the Evolutionary Origins of Human Behavior, Kyoto University, Inuyama 484-8506, Japan
2
NIBB Core Research Facilities, National Institute for Basic Biology, Okazaki 444-0867, Japan
3
Research Center for Integrative Evolutionary Science, SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193, Japan
*
Author to whom correspondence should be addressed.
Animals 2024, 14(16), 2316; https://doi.org/10.3390/ani14162316
Submission received: 2 July 2024 / Revised: 29 July 2024 / Accepted: 30 July 2024 / Published: 9 August 2024
(This article belongs to the Special Issue Recent Advances in Genetic Mechanism and Evolution of Sex Determination and Chromosomes)
Browse Figures
Figure 1
<p>Acquisition/loss of testis-expressed genes in 10 species. The genes at each clade are listed in <a href="#animals-14-02316-t001" class="html-table">Table 1</a>, <a href="#animals-14-02316-t002" class="html-table">Table 2</a>, <a href="#app1-animals-14-02316" class="html-app">Tables S2 and S3</a>, and the left and right of the ‘/’ symbols mean acquisition and loss of testis-expressed genes, respectively.</p> "> Figure 2
<p>The phylogenetic trees, branch lengths, and dn/ds ratios for three therian-specific testis genes in 10 species. (<b>A</b>) <span class="html-italic">ARHGAP28</span> genes (1833 bp). (<b>B</b>) <span class="html-italic">SYNM</span> genes (2235 bp). (<b>C</b>) <span class="html-italic">PDZRN3</span> genes (2256 bp). The trees were constructed by the NJ method using the number of nucleotide differences. The average and standard deviation of branch lengths in six eutherians are shown above the black bold branch, and those in two marsupials are shown above the gray bold branch. The dn/ds ratio in six eutherian or two marsupial pairs is shown under the branch. ** and * mean <span class="html-italic">p</span> &lt; 0.01 and <span class="html-italic">p</span> &lt; 0.05, respectively, and are supported by Fisher’s exact test.</p> "> Figure 2 Cont.
<p>The phylogenetic trees, branch lengths, and dn/ds ratios for three therian-specific testis genes in 10 species. (<b>A</b>) <span class="html-italic">ARHGAP28</span> genes (1833 bp). (<b>B</b>) <span class="html-italic">SYNM</span> genes (2235 bp). (<b>C</b>) <span class="html-italic">PDZRN3</span> genes (2256 bp). The trees were constructed by the NJ method using the number of nucleotide differences. The average and standard deviation of branch lengths in six eutherians are shown above the black bold branch, and those in two marsupials are shown above the gray bold branch. The dn/ds ratio in six eutherian or two marsupial pairs is shown under the branch. ** and * mean <span class="html-italic">p</span> &lt; 0.01 and <span class="html-italic">p</span> &lt; 0.05, respectively, and are supported by Fisher’s exact test.</p> ">
Review Reports Versions Notes

Abstract

:

Simple Summary

Simple Summary: A transcriptome landscape of therian mammals is known, but it remains unclear how transcriptomic patterns have evolved in a marsupial- or eutherian-specific way. It is important to understand marsupial- or eutherian-specific transcriptomic patterns since their fitness and sex differentiation are different. This study examines therian testis transcriptomes to elucidate marsupial and eutherian uniqueness in male differentiation. Using the massive transcriptomic data, we show the evolutionary tempo and mode of testis-expressed genes in Theria and identify candidate genes involved in the specificity of marsupial or eutherian testes.

Abstract

Gene expression patterns differ in different tissues, and the expression pattern of genes in the mammalian testis is known to be extremely variable in different species. To clarify how the testis transcriptomic pattern has evolved in particular species, we examined the evolution of the adult testis transcriptome in Theria using 10 species: two marsupials (opossum and Tasmanian devil), six eutherian (placental) mammals (human, chimpanzee, bonobo, gorilla, rhesus macaque, and mouse), and two outgroup species (platypus and chicken). We show that 22 testis-expressed genes are marsupial-specific, suggesting their acquisition in the stem lineage of marsupials after the divergence from eutherians. Despite the time length of the eutherian stem lineage being similar to that of the marsupial lineage, acquisition of testis-expressed genes was not found in the stem lineage of eutherians; rather, their expression patterns differed by species, suggesting rapid gene evolution in the eutherian ancestors. Fifteen testis-expressed genes are therian-specific, and for three of these genes, the evolutionary tempo is markedly faster in eutherians than in marsupials. Our phylogenetic analysis of Rho GTPase-activating protein 28 (ARHGAP28) suggests the adaptive evolution of this gene in the eutherians, probably together with the expression pattern differentiation.

1. Introduction

Theria is a mammalian subclass composed of Eutheria (eutherian or placental mammals) and Metatheria (marsupials or pouched mammals), and the molecular clock suggests that these two infraclasses diverged 148 million years ago [1]. Eutheria have spread around the world, but marsupials have lived only in limited places such as the Australian and American continents, suggesting that the two groups have adapted to particular ecological environments. Differences in fitness to these environments should result from the divergence of genome and transcriptome between eutherians and marsupials. The transcriptome of reproductive tissues is directly related to the efficiency of producing offspring, and it is thus important to understand transcriptomic differences in reproductive tissues as the evolution of the directly influential system in fitness.
The therian Y chromosome is a testis determiner and possesses a sex determination gene (SRY). In marsupials, the molecular system of sex determination and differentiation is not as well understood as that in eutherians. We have shown that the marsupial SRY gene is functionally similar to eutherian SRY and upregulates SRY-box9 (SOX9) synergistically with steroidogenic factor 1 (SF1) [2]. This means that the system of sex determination is common in marsupials and eutherians, but the differentiation of sexual characteristics is distinct in the two groups. For example, in marsupials, gonadal differentiation occurs after birth, but in eutherians, the gonads differentiate in the embryo before birth [3]. Moreover, marsupial mammary glands and scrotum, which are sexually dimorphic structures, have usually started to differentiate before ovary and testis differentiation, independent of endogenous estrogen or testosterone [3,4,5,6,7]. Graves, Renfree, and their colleagues suggest that the marsupial Y chromosome does not control all aspects of sex differentiation, as it does in eutherians through the initiation of androgen pathways that control most sexual dimorphisms [4,6,8].
RNA sequencing (RNA-seq) using next-generation sequencers is useful for understanding the entire profile of gene expression in different tissues. Marsupial transcriptomes have been investigated by several groups [9,10,11]. Brawand et al. (2011) investigated transcriptomes of six tissues in 10 species including opossum (Monodelphis domestica) and showed rapid evolution of the transcriptome in testis compared to other tissues [9]. Murat et al. (2023) published transcriptomes of testes using single-cell RNA-seq and showed temporal expression changes of genes in spermatogenesis across 11 species including opossum [12]. These previous studies revealed a transcriptome landscape of mammals, but it remains unclear how transcriptomic patterns have evolved in a marsupial-specific way. In particular, we do not know which transcriptomic characteristics were acquired or lost in the stem lineage of marsupials.
In this study, we examine therian testis transcriptomes to elucidate marsupial and eutherian uniqueness in adult male differentiation. Using the massive transcriptomic data, it is now possible to understand what evolutionary events happened during therian evolution. Thus, we show the evolutionary tempo and mode of testis-expressed genes in Theria and identify candidate genes involved in the specificity of marsupial or eutherian testes.

2. Materials and Methods

2.1. Sample and Library Preparation

We used two RNA samples of Tasmanian devil testis from a single individual; these samples were kind gifts from Dr. Jenny Graves. RNA was extracted with a GenElute™ Mammalian Total RNA Miniprep Kit (Sigma, Castle Hill, NSW, Australia). cDNA libraries were generated from purified RNA (1 μg of each sample) using a TruSeq RNA Sample Preparation Kit v2 (Illumina, San Diego, CA, USA).

2.2. RNA Sequencing and Assembling

The two multiplexed libraries were sequenced in a single lane using Hiseq2000 (Illumina) with 101 bp paired-end reads. The short reads were modified by removing adapter sequences and were trimmed to remove low-quality sequences shorter than 64 bp and below QV (quality value) 20. The cleaned high-quality reads were mapped on the Tasmanian devil nuclear genome (Devil_ref v7.0: GCA_000189315.1) and mitochondrion genome (FN666604) using bowtie2 software [13], with three mismatches being allowed in the software. The quality of mapping sequences was evaluated by SAMtools [14], and the mapped sequences were visualized using Integrative Genomics Viewer [15]. The short-read sequences were also de novo-assembled using Trinity [16].
Transcripts were assembled and annotated using Cufflinks [17] and quantified with RSEM software [18]. Fragments per kb of exon per million mapped fragments (FPKM) and reads per kb of exon per million mapped reads (RPKM) were calculated, and genes expressed differentially between samples were found using Cuffdiff [17].
Table S4 shows the number of RNA-seq reads in the Tasmanian devil testis and brain. We identified 42,964 genes in total, of which 16,438 transcripts were previously reported by Murchison et al. (2010) [10]. Of the 42,964 genes, 20,414 loci are registered in Ensembl.

2.3. Comparative Evolutionary Analyses

We identified 1:1 orthologous genes among 10 species and created a testis expression gene list of 5627 orthologous genes expressed in the testis from at least one other of the 10 species. The orthologous list was created using Ensembl biomart. The 10 species are Tasmanian devil in this study and nine species from a previous study [9]: human (Homo sapiens), chimpanzee (Pan troglodytes), bonobo (Pan paniscus), gorilla (Gorilla gorilla), rhesus macaque (Macaca mulatta), mouse (Mus musculus), opossum (Monodelphis domestica), platypus (Ornithorhynchus anatinus), and chicken (Gallus gallus). If the RPKM value of a gene is <1, the gene is considered not to be expressed [19]. Thus, we assume that RPKM < 1 indicates loss of the gene’s expression, and RPKM ≥ 1 indicates acquisition of the gene’s expression. All RPKM values can be requestable. Of the 5627 genes, 2989 are expressed in the testes of all 10 species (RPKM ≥ 1).
Our comparison showed that the 22 genes acquired testis expression in the stem lineage of marsupials based on bulk RNA-seq data. In human testes, however, expression of one of the 22 marsupial testis-expressed genes, PRDM1 (1.2 nTPM: normalized transcripts per million), is shown, and medium (25–75% fraction of stained cells) protein expression is reported as well (the Human Protein Atlas: https://www.proteinatlas.org/, accessed on 29 November 2023). Since we categorized the expressed genes into a binary state (present or non-present) based on an expression level (RPKM ≥ 1 or RPKM < 1), genes situated in the cut-off limit could be identified as non-expressed. In PRDM1, for example, its eutherian RPKM is 0.53 ± 0.34 (human; RPKM = 0.80 ± 0.50, TPM = 1.75 ± 0.37), suggesting it is not likely expressed in this category, although the expression level is slightly below the cut-off. Using this cut-off, we found that the expression level of eutherian PRDM1 is much lower than the marsupial expression (RPKM = 192.02 ± 130.21), while IGFBP1 is detected in Sertori cells by only human single-cell RNA sequencing since the single-cell RNA-seq is known to be highly sensitive, more than the bulk RNA-seq (the Human Protein Atlas).

2.4. Molecular Evolutionary Analyses

The translated nucleotide sequences were aligned using multiple sequence comparisons by log expectation (MUSCLE) [20], and the entire alignment is available in the Supplementary Materials (data sheet 1–4). Phylogenetic trees were constructed using neighbor-joining (NJ) and maximum likelihood (ML) methods implemented in the MEGA7 program [21]. The reliability of the trees was assessed by bootstrap re-sampling with 1000 replications.

3. Results and Discussion

3.1. The Evolutionary Tempo and Mode of Testis Genes in Theria

First, we identified 2989 testis-expressed genes that are evolutionary conserved in 10 species using RNA-seq data (see Section 2). Of 2989 testis-expressed genes in mammals as well as chickens, 53 are expressed in testes of the mammals (eutherians, marsupials, and monotremes) only (Section 2, Figure 1, and Table S1). Of 53, 22 are expressed in adult testes of the marsupials, and this means that the testis expression of these genes was acquired in the stem lineage of marsupials (Figure 1 and Table 1). The 22 genes are located on the autosome and include a serine protease (TRAF3IP3), protease inhibitors (IGFBP1, SERPINA10), a glycosyltransferase (FUT9), a zinc finger transcription factor (BCL11B), an oxidase (TYRP1), an actin-binding protein (LMOD3), a deaminase (AMPD1), and a signaling molecule (DAB1) (Table 1).
There are no genes expressed commonly in the testes of only the eutherians (Figure 1). One reason for this may be a rapid evolution of gene expression in eutherians. If gene expression has changed rapidly in one group, common gene expression in all groups cannot be found from the present transcriptomic information. To test whether the rapid evolution of genes has occurred in eutherians, we investigated phylogenetic trees of genes shared by Theria. Fifteen genes were expressed in the testes of only the therian species, and we inferred that this testis expression pattern had been obtained in the common ancestor of Theria (Table 2). We constructed phylogenetic trees of the 15 genes in 10 species used and calculated branch lengths in the eutherian and marsupial stem clades. In three genes (ARHGAP28, SYNM, and PDZRN3), the average of the eutherian branch lengths is significantly greater than that of the marsupial branch lengths (Figure 2A–C; Fisher’s exact test, p < 0.05). The longer branch lengths reflect faster gene evolution in Eutheria. In Rho GTPase-activating protein 28 (ARHGAP28) the dn/ds ratio (0.39) at the eutherian stem lineage is significantly higher than 0.10 at the marsupial stem lineage (Figure 2A; Fisher’s exact test, p < 0.01). In the other two genes, the ratios are not different statistically (Figure 2B,C).
These observations are consistent with previous reports: in rodents and primates, the evolution of testis-expressed genes is accelerated compared to that of genes expressed in other tissues [22]. In our study, the rapid evolution of genes is not observed in marsupials. The evolutionary tempo and mode of testis-expressed genes may be different in marsupials and eutherians.

3.2. Therian-Specific Testis-Expressed Genes

Our results support the hypothesis that at least 15 genes obtained testis expression in the common ancestor of Theria, and then three (ARHGAP28, SYNM, and PDZRN3) evolved rapidly in only the eutherians (Table 2, Figure 1 and Figure 2). One reason is that the mutation rate at the genome level is lower in marsupials than in eutherians. Feigin et al. (2018) estimated that the mutation rate is 1.17 × 10−9 mutations per base per generation in the Tasmanian devil genome, compared with 6.98 × 10−8 mutations per base per generation in the human genome [23]. Our phylogenetic analyses suggest that the evolutionary rate of the three genes did not slow down in the stem lineage of marsupials compared to the ancestral state of Theria; rather, the evolutionary rate increased in the stem lineage of eutherians. However, in ARHGAP28 the increasing rate (longer branch in eutherians) is more likely the acceleration of the nonsynonymous substitution rate due to adaptive evolution in the stem lineage of eutherians. The phylogenetic tree using the amino acid sequences shows that the branch at the eutherian stem is six times longer than that at the marsupial stem (Figure S1; Fisher’s exact test, p < 0.01), and this is consistent with the higher dn/ds ratio at the branch of the eutherian stem (Figure 2A). ARHGAP28 is one of the GTPase-activating proteins (GAPs) that activate small monomeric Rho GTPases, and it is highly expressed in developing eutherian testes [24,25,26]. In the eutherian ancestor, ARHGAP28 may have acquired adaptive changes that enabled it to control molecules expressed in testis but the precise function of ARHGAP28 in testis has not yet been studied.
Synemin (SYNM) is a member of the intermediate filament family and encodes a 230-kDa polypeptide having a cytoskeletal role mainly in muscle [27]. PDZ domain containing ring finger 3 (PDZRN3) is a ubiquitin ligase required for vascular morphogenesis in endothelium [28]. These ubiquitous genes are likely expressed not only in testis (mRNA/protein data available from the Human Protein Atlas database, accessed on 27 October 2022), although they probably have a fundamental function in testes. Cytochrome P450 family 2 subfamily R member 1 (CYP2R1) is also a ubiquitous gene listed in the therian testis-expressed genes (Table 2). CYP2R1, also known as vitamin D 25-hydroxylase, acts in the vitamin D synthesis pathway and drug metabolism, and vitamin D is positively associated with sperm motility, affecting male fertility in humans [29].

3.3. 22 Genes: Functional Differentiation in a Member of a Large Gene Family

Of the 22 marsupial and 15 therian testis-expressed genes, four are type II cadherin genes (CDH7, 18, 19, and 20). CDH9 is also one of the platypus testis-expressed genes (Table S2). Cadherin is a transmembrane or membrane-associated glycoprotein that mediates cell adhesion in a Ca2+-dependent manner and is essential for spermatogenesis [30]. Type II cadherin proteins do not contain a conserved cell adhesion recognition sequence, His-Ala-Val, in their ectodomain module [31]. Multiple members of the cadherin family (CDH6, 8, 10, and 11) are expressed in the rat testis [32] and are involved in germ cell adhesion to Sertoli cells [33,34]. However, the function of type II cadherin in testes is not known.
As noted above, protease inhibitors (IGFBP1, SERPINA10) and a serine protease (TRAF3IP3) are included in the marsupial testis-expressed genes and may be required for spermatogenesis [35]. SERPINA comprises a large gene family, one of whose members SERPINA5 (also known as protein C inhibitor) is known to be upregulated by testosterone [35]. Testis expression of the gene is shared only by therian mammals, as the chicken SERPINA5 is expressed in the liver and kidney. While the marsupial SERPINA10 is expressed in testes, the human and chicken genes are highly expressed in the liver (AceView) [36]. This implies that in the therian ancestor, SERPINA5 acquired a new function as a testis-expressed gene, and SERPINA10 became a testis-specific protease in the marsupial ancestor. A 5,6-dihydroxyindole-2-carboxylic acid oxidase (TYRP1) and an actin-binding protein (LMOD3) are the two other marsupial testis-expressed genes. It has been suggested that male infertility is caused by oxidative stress and downregulation or mutations of actin-related genes (e.g., ACTL7B) in sperms [37,38,39].
These 22 genes newly acquired testis expression in the stem lineage of marsupials, but their function in testis is not yet known. Each of the genes is a member of a gene family, one of whose other members reportedly have functions in testis differentiation and spermatogenesis of eutherians. The particular subsets of gene family members that function in testis may differ slightly between marsupials and eutherians, and these functional differences between genes may contribute to the different patterns of sex differentiation observed in marsupials and eutherians.

3.4. Perspectives on the Function of PRDM1

Using previous data [40], we investigated whether the 22 marsupial-adult testis-expressed genes are expressed in opossum newborn testes, and one gene, PR/SET Domain 1 (PRDM1), was also found to be expressed in the neonatal period (postnatal 6 days; RPKM = 18.59 ± 10.41). PRDM is a gene family with an ancient origin, and most of its members emerged early in animal evolution [41,42]. The number of PRDM genes is 13–15 in marsupials and 12–17 in eutherians [41]. The paralogues to PRDM1 are crucial for the development of testis. PRDM14 is expressed in mouse PGCs, and its knockout mice are infertile [43]. PRDM8 regulates mouse testis steroidogenesis by repressing p450c17 (CYP17A1) and luteinizing hormone receptor (LHR) [44]. The accelerated evolution of PRDM9 is supported by site-specific positive selection [45], but we could not identify any statistically significant positive selection sites in PRDM1.
PRDM1 is a transcription factor with five zinc finger motifs and represses β-interferon (β-INF) gene expression by binding to the positive regulatory domain 1 (PRD1) element of the β-INF promoter [46]. In humans and mice, PRDM1 is highly expressed in immune system cells such as B/T lymphocytes and natural killer (NK) cells. In mouse studies, β-INF overexpression in the testis disrupts spermatogenesis [47]. PRDM1 is also known as B-lymphocyte-induced maturation protein 1 (Bimp1) and is crucial for primordial germ cell (PGC) differentiation [48]. The function of PRDM1 in testes is not well understood, and their functional analysis is necessary in the future.

4. Conclusions

We show the evolutionary tempo of three genes (ARHGAP28, SYNM, and PDZRN3) is fast in eutherians compared to marsupials, and especially the amino acid substitutions of ARHGAP28 are accumulated in the common ancestor of eutherians. Those three genes are probably involved in the specificity of eutherian testes. We also identified 22 genes that are involved in the specificity of marsupial adult testes with high expression based on the comparative bulk RNA-seq data. We observed the functional differentiation and sub-functionalization, that is, marsupial-specific high testis expression, of a member of a large gene family such as PRDM, CDH, and SERPINA.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ani14162316/s1: Figure S1: Phylogenetic tree of ARHGAP28 amino acid sequences; Table S1: The number of genes with RPKM > 1 in 10 species; Table S2: Acquisition of testis expression in each clade; Table S3: Loss of testis expression in each clade; Table S4: RNA-seq data for the Tasmanian devil.

Author Contributions

Conceptualization, Y.K.; methodology, Y.K. and S.S.; validation, Y.K., S.S. and Y.S.; formal analysis, Y.K. and S.S.; investigation, Y.K.; resources, Y.K. and S.S.; data curation, Y.K. and S.S.; writing—original draft preparation, Y.K.; writing—review and editing, Y.K., S.S. and Y.S.; visualization, Y.K.; supervision, Y.K.; project administration, Y.K.; funding acquisition, Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by JSPS KAKENHI Grant Number 22H05070 and 22H05069.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used were deposited in a repository (DRA016090).

Acknowledgments

We thank Jenny Graves for sample collection, Katsushi Yamaguchi for supporting RNA-seq sequencing, Erich Smith for data organization, and Henrik Kaessmann and Margarida Cardoso Moreira for sharing expression data. This work was supported by Trans-Omics Facility, NIBB Trans-Scale Biology Center. We thank Trans-Omics Facility, NIBB Trans-Scale Biology Center for technical support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bininda-Emonds, O.R.P.; Cardillo, M.; Jones, K.E.; MacPhee, R.D.E.; Beck, R.M.D.; Grenyer, R.; Price, S.A.; Vos, R.A.; Gittleman, J.L.; Purvis, A. The delayed rise of present-day mammals. Nature 2007, 446, 507–512. [Google Scholar] [CrossRef]
  2. Katsura, Y.; Kondo, H.X.; Ryan, J.; Harley, V.; Satta, Y. The Evolutionary Process of Mammalian Sex Determination Genes focusing on Marsupial SRYs. BMC Evol. Biol. 2018, 18, 3. [Google Scholar]
  3. Renfree, M.B.; Chew, K.Y.; Shaw, G. Inducing sex reversal of the urogenital system of marsupials. Differentiation 2014, 87, 23–31. [Google Scholar] [CrossRef] [PubMed]
  4. Graves, J.A.M.; Renfree, M.B. Marsupials in the Age of Genomics. Annu. Rev. Genom. Hum. Genet. 2013, 14, 393–420. [Google Scholar] [CrossRef] [PubMed]
  5. O, W.S.; Short, R.V.; Renfree, M.B.; Shaw, G. Primary genetic control of somatic sexual differentiation in a mammal. Nature 1988, 331, 716–717. [Google Scholar] [CrossRef] [PubMed]
  6. Renfree, M.B.; Shaw, G. Diapause. Annu. Rev. Physiol. 2000, 62, 353–375. [Google Scholar] [CrossRef]
  7. Shaw, G.; Renfree, M.B.; Short, R.V. OWS Experimental manipulation of sexual differentiation in wallaby pouch young treated with exogenous steroids. Development 1988, 104, 689–701. [Google Scholar] [CrossRef] [PubMed]
  8. Cooper, D.; Johnston, P.; Watson, J.; Graves, J.A.M. X-inactivation in marsupials and monotremes. Semin. Dev. Biol. 1993, 4, 117–128. [Google Scholar] [CrossRef]
  9. Brawand, D.; Soumillon, M.; Necsulea, A.; Julien, P.; Csardi, G.; Harrigan, P.; Weier, M.; Liechti, A.; Aximu-Petri, A.; Kircher, M.; et al. The evolution of gene expression levels in mammalian organs. Nature 2011, 478, 343–348. [Google Scholar] [CrossRef]
  10. Murchison, E.P.; Tovar, C.; Hsu, A.; Bender, H.S.; Kheradpour, P.; Rebbeck, C.A.; Obendorf, D.; Conlan, C.; Bahlo, M.; Blizzard, C.A.; et al. The Tasmanian devil transcriptome reveals Schwann cell origins of a clonally transmissible cancer. Science 2010, 327, 84–87. [Google Scholar] [CrossRef]
  11. Renfree, M.B.; Papenfuss, A.T.; Deakin, J.E.; Lindsay, J.; Heider, T.; Belov, K.; Rens, W.; Waters, P.D.; Pharo, E.A.; Shaw, G.; et al. Genome sequence of an Australian kangaroo, Macropus eugenii, provides insight into the evolution of mammalian reproduction and development. Genome. Biol. 2011, 12, R81. [Google Scholar] [CrossRef]
  12. Murat, F.; Mbengue, N.; Winge, S.B.; Trefzer, T.; Leushkin, E.; Sepp, M.; Cardoso-Moreira, M.; Schmidt, J.; Schneider, C.; Mößinger, K.; et al. The molecular evolution of spermatogenesis across mammals. Nature 2023, 613, 308–316. [Google Scholar] [CrossRef] [PubMed]
  13. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef]
  14. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. 1000 Genome Project Data Processing Subgroup. The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  15. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative Genomics Viewer. Nat. Biotechnol. 2011, 29, 24–26. [Google Scholar] [CrossRef] [PubMed]
  16. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  17. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed]
  18. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed]
  19. Ambrosio, A.B.; Nascimento, L.C.D.; Oliveira, B.V.; Teixeira, P.J.P.L.; A Tiburcio, R.; Thomazella, D.P.T.; Leme, A.F.P.; Carazzolle, M.F.; O Vidal, R.; Mieczkowski, P.; et al. Global analyses of Ceratocystis cacaofunesta mitochondria: From genome to proteome. BMC Genom. 2013, 14, 91. [Google Scholar]
  20. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef]
  21. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  22. Khaitovich, P.; Enard, W.; Lachmann, M.; Pääbo, S. Evolution of primate gene expression. Nat. Rev. Genet. 2006, 7, 693–702. [Google Scholar] [CrossRef]
  23. Feigin, C.Y.; Newton, A.H.; Doronina, L.; Schmitz, J.; Hipsley, C.A.; Mitchell, K.J.; Gower, G.; Llamas, B.; Soubrier, J.; Heider, T.N.; et al. Genome of the Tasmanian tiger provides insights into the evolution and demography of an extinct marsupial carnivore. Nat. Ecol. Evol. 2018, 2, 182–192. [Google Scholar]
  24. Boureux, A.; Vignal, E.; Faure, S.; Fort, P. Evolution of the Rho family of ras-like GTPases in eukaryotes. Mol. Biol. Evol. 2007, 24, 203–216. [Google Scholar] [CrossRef]
  25. Etienne-Manneville, S.; Hall, A. Rho GTPases in cell biology. Nature 2002, 420, 629–635. [Google Scholar] [CrossRef]
  26. Van Aelst, L.; D’Souza-Schorey, C. Rho GTPases and signaling networks. Genes Dev. 1997, 11, 2295–2322. [Google Scholar] [PubMed]
  27. Granger, B.L.; Lazarides, E. Synemin: A new high molecular weight protein associated with desmin and vimentin filaments in muscle. Cell 1980, 22, 727–738. [Google Scholar] [CrossRef] [PubMed]
  28. Sewduth, R.N.; Jaspard-Vinassa, B.; Peghaire, C.; Guillabert, A.; Franzl, N.; Larrieu-Lahargue, F.; Moreau, C.; Fruttiger, M.; Dufourcq, P.; Couffinhal, T.; et al. The ubiquitin ligase PDZRN3 is required for vascular morphogenesis through Wnt/planar cell polarity signalling. Nat. Commun. 2014, 5, 4832. [Google Scholar] [CrossRef] [PubMed]
  29. de Angelis, C.; Galdiero, M.; Pivonello, C.; Garifalos, F.; Menafra, D.; Cariati, F.; Salzano, C.; Galdiero, G.; Piscopo, M.; Vece, A.; et al. The role of vitamin D in male fertility: A focus on the testis. Rev. Endocr. Metab. Disord. 2017, 18, 285–305. [Google Scholar] [CrossRef] [PubMed]
  30. Goossens, S.; Van Roy, F. Cadherin-mediated cell-cell adhesion in the testis. Front. Biosci. 2005, 10, 398–419. [Google Scholar] [CrossRef]
  31. Nollet, F.; Kools, P.; van Roy, F. Phylogenetic analysis of the cadherin superfamily allows identification of six major subfamilies besides several solitary members. J. Mol. Biol. 2000, 299, 551–572. [Google Scholar] [PubMed]
  32. Johnson, K.J.; Patel, S.R.; Boekelheide, K. Multiple cadherin superfamily members with unique expression profiles are produced in rat testis. Endocrinology 2000, 141, 675–683. [Google Scholar] [CrossRef] [PubMed]
  33. Byers, S.W.; Sujarit, S.; Jegou, B.; Butz, S.; Hoschutzky, H.; Herrenknecht, K.; MacCalman, C.; Blaschuk, O.W. Cadherins and cadherin-associated molecules in the developing and maturing rat testis. Endocrinology 1994, 134, 630–639. [Google Scholar] [CrossRef]
  34. Kopera, I.A.; Bilinska, B.; Cheng, C.Y.; Mruk, D.D. Sertoli-germ cell junctions in the testis: A review of recent data. Philos. Trans. R Soc. Lond. B Biol. Sci. 2010, 365, 1593–1605. [Google Scholar] [CrossRef] [PubMed]
  35. Le Magueresse-Battistoni, B. Serine proteases and serine protease inhibitors in testicular physiology: The plasminogen activation system. Reproduction 2007, 134, 721–729. [Google Scholar] [CrossRef] [PubMed]
  36. Thierry-Mieg, D.; Thierry-Mieg, J. AceView: A comprehensive cDNA-supported gene and transcripts annotation. Genome Biol. 2006, 7, S12. [Google Scholar] [CrossRef] [PubMed]
  37. Ahn, J.; Park, Y.-J.; Chen, P.; Lee, T.J.; Jeon, Y.-J.; Croce, C.M.; Suh, Y.; Hwang, S.; Kwon, W.-S.; Pang, M.-G.; et al. Comparative expression profiling of testis-enriched genes regulated during the development of spermatogonial cells. PLoS ONE 2017, 12, e0175787. [Google Scholar] [CrossRef] [PubMed]
  38. Clement, T.; Geyer, C.; Willis, W.; Goulding, E.; Eddy, M. Expressed Actin-Like 7b (Actl7b) Is Required for Mouse Spermatid Morphogenesis and Male Fertility. Biol. Reprod. 2012, 87, 140. [Google Scholar] [CrossRef]
  39. Hosen, M.B.; Islam, M.R.; Begum, F.; Kabir, Y.; Howlader, M.Z. Oxidative stress induced sperm DNA damage, a possible reason for male infertility. Iran. J. Reprod. Med. 2015, 13, 525–532. [Google Scholar]
  40. Cardoso-Moreira, M.; Halbert, J.; Valloton, D.; Velten, B.; Chen, C.; Shao, Y.; Liechti, A.; Ascenção, K.; Rummel, C.; Ovchinnikova, S.; et al. Gene expression across mammalian organ development. Nature 2019, 571, 505–509. [Google Scholar] [CrossRef]
  41. Fumasoni, I.; Meani, N.; Rambaldi, D.; Scafetta, G.; Alcalay, M.; Ciccarelli, F.D. Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates. BMC Evol. Biol. 2007, 7, 187. [Google Scholar] [CrossRef] [PubMed]
  42. Vervoort, M.; Meulemeester, D.; Béhague, J.; Kerner, P. Evolution of Prdm Genes in Animals: Insights from Comparative Genomics. Mol. Biol. Evol. 2016, 33, 679–696. [Google Scholar] [CrossRef] [PubMed]
  43. Yamaji, M.; Seki, Y.; Kurimoto, K.; Yabuta, Y.; Yuasa, M.; Shigeta, M.; Yamanaka, K.; Ohinata, Y.; Saitou, M. Critical function of Prdm14 for the establishment of the germ cell lineage in mice. Nat. Genet. 2008, 40, 1016–1022. [Google Scholar] [CrossRef] [PubMed]
  44. Eom, G.H.; Kim, K.; Kim, S.-M.; Kee, H.J.; Kim, J.-Y.; Jin, H.M.; Kim, J.-R.; Kim, J.H.; Choe, N.; Kim, K.-B.; et al. Histone methyltransferase PRDM8 regulates mouse testis steroidogenesis. Biochem. Biophys. Res. Commun. 2009, 388, 131–136. [Google Scholar] [CrossRef] [PubMed]
  45. Oliver, P.L.; Goodstadt, L.; Bayes, J.J.; Birtle, Z.; Roach, K.C.; Phadnis, N.; Beatson, S.A.; Lunter, G.; Malik, H.S.; Ponting, C.P. Accelerated Evolution of the Prdm9 Speciation Gene across Diverse Metazoan Taxa. PLoS Genet. 2009, 5, e1000753. [Google Scholar] [CrossRef] [PubMed]
  46. Keller, A.D.; Maniatis, T. Identification and characterization of a novel repressor of beta-interferon gene expression. Genes Dev. 1991, 5, 868–879. [Google Scholar] [CrossRef]
  47. Satie, A.P.; Mazaud-Guittot, S.; Seif, I.; Mahé, D.; He, Z.; Jouve, G.; Jégou, B.; Dejucq-Rainsford, N. Excess type I IFN signaling in the mouse seminiferous tubules leads to germ cell loss and sterility. J. Biol. Chem. 2011, 286, 23280–23295. [Google Scholar] [CrossRef]
  48. Ohinata, Y.; Payer, B.; O’Carroll, D.; Ancelin, K.; Ono, Y.; Sano, M.; Barton, S.C.; Obukhanych, T.; Nussenzweig, M.; Tarakhovsky, A.; et al. Blimp1 is a critical determinant of the germ cell lineage in mice. Nature 2005, 436, 207–213. [Google Scholar] [CrossRef]
Animals 14 02316 g001
Figure 1. Acquisition/loss of testis-expressed genes in 10 species. The genes at each clade are listed in Table 1, Table 2, Tables S2 and S3, and the left and right of the ‘/’ symbols mean acquisition and loss of testis-expressed genes, respectively.
Figure 1. Acquisition/loss of testis-expressed genes in 10 species. The genes at each clade are listed in Table 1, Table 2, Tables S2 and S3, and the left and right of the ‘/’ symbols mean acquisition and loss of testis-expressed genes, respectively.
Animals 14 02316 g001
Animals 14 02316 g002aAnimals 14 02316 g002b
Figure 2. The phylogenetic trees, branch lengths, and dn/ds ratios for three therian-specific testis genes in 10 species. (A) ARHGAP28 genes (1833 bp). (B) SYNM genes (2235 bp). (C) PDZRN3 genes (2256 bp). The trees were constructed by the NJ method using the number of nucleotide differences. The average and standard deviation of branch lengths in six eutherians are shown above the black bold branch, and those in two marsupials are shown above the gray bold branch. The dn/ds ratio in six eutherian or two marsupial pairs is shown under the branch. ** and * mean p < 0.01 and p < 0.05, respectively, and are supported by Fisher’s exact test.
Figure 2. The phylogenetic trees, branch lengths, and dn/ds ratios for three therian-specific testis genes in 10 species. (A) ARHGAP28 genes (1833 bp). (B) SYNM genes (2235 bp). (C) PDZRN3 genes (2256 bp). The trees were constructed by the NJ method using the number of nucleotide differences. The average and standard deviation of branch lengths in six eutherians are shown above the black bold branch, and those in two marsupials are shown above the gray bold branch. The dn/ds ratio in six eutherian or two marsupial pairs is shown under the branch. ** and * mean p < 0.01 and p < 0.05, respectively, and are supported by Fisher’s exact test.
Animals 14 02316 g002aAnimals 14 02316 g002b
Table 1. Marsupial-specific testis-expressed genes.
Table 1. Marsupial-specific testis-expressed genes.
NIPAL1NIPA-like domain containing 1
LMOD3Leiomodin 3
FUT9Fucosyltransferase 9
TECTBTectorin beta
AMPD1Adenosine monophosphate deaminase 1
TYRP1Tyrosinase-related protein 1
PRDM1PR/SET domain 1
BCL11BB-cell CLL/lymphoma 11B
BEST3Bestrophin 3
IHHIndian hedgehog
IKBKEInhibitor of nuclear factor kappa B kinase subunit epsilon
DAB1DAB1, reelin adaptor protein
TRAF3IPTRAF3-interacting protein 3
RBM20RNA-binding motif protein 20
SERPINA10Serpin family A member 10
CDH7Cadherin 7
CDH19Cadherin 19
CDH20Cadherin 20
IGFBP1Insulin-like growth factor binding protein 1
ZNF750Zinc finger protein 750
THEMISThymocyte selection associated
NTSNeurotensin
Table 2. Therian-specific testis-expressed genes.
Table 2. Therian-specific testis-expressed genes.
ARHGAP28Rho GTPase-activating protein 28
CYP2R1Cytochrome P450 family 2 subfamily R member 1
ESRP2Epithelial splicing regulatory protein 2
ATP9AATPase phospholipid transporting 9A
MYO15AMyosin XVA
HSPA12AHeat shock protein family A (Hsp70) member 12A
TOGARAM2TOG array regulator of axonemal microtubules 2
PDZRN3PDZ domain containing ring finger 3
SLC7A2Solute carrier family 7 member 2
SYNMSynemin
CDH18Cadherin 18
GJA5Gap junction protein, alpha 5
SH3RF2SH3 domain containing ring finger 2
WFDC1WAP four-disulfide core domain 1
ELFN2Extracellular leucine rich repeat and fibronectin type III domain containing 2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Katsura, Y.; Shigenobu, S.; Satta, Y. Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria. Animals 2024, 14, 2316. https://doi.org/10.3390/ani14162316

AMA Style

Katsura Y, Shigenobu S, Satta Y. Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria. Animals. 2024; 14(16):2316. https://doi.org/10.3390/ani14162316

Chicago/Turabian Style

Katsura, Yukako, Shuji Shigenobu, and Yoko Satta. 2024. "Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria" Animals 14, no. 16: 2316. https://doi.org/10.3390/ani14162316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop