Male-specific region of the bovine Y chromosome is gene rich with a high transcriptomic activity in testis development

Comparative genomics illuminates karyotype and sex chromosome evolution of sharks

Chondrichthyes is an important lineage to reconstruct the evolutionary history of vertebrates. Here, we analyzed genome synteny for six chondrichthyan chromosome-level genomes. Our comparative analysis reveals a slow evolutionary rate of chromosomal changes, with infrequent but independent fusions observed in sharks, skates, and chimaeras. The chondrichthyan common ancestor had a proto-vertebrate-like karyotype, including the presence of 18 microchromosome pairs. The X chromosome is a conversed microchromosome shared by all sharks, suggesting a likely common origin of the sex chromosome at least 181 million years ago. We characterized the Y chromosomes of two sharks that are highly differentiated from the X except for a small young evolutionary stratum and a small pseudoautosomal region. We found that shark sex chromosomes lack global dosage compensation but that dosage-sensitive genes are locally compensated. Our study on shark chromosome evolution enhances our understanding of shark sex chromosomes and vertebrate chromosome evolution.

The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Synergistic enhancement of the mouse Pramex1 and Pramel1 in repressing retinoic acid (RA) signaling during gametogenesis

BACKGROUND

PRAME constitutes one of the largest multi-copy gene families in Eutherians, encoding cancer-testis antigens (CTAs) with leucine-rich repeats (LRR) domains, highly expressed in cancer cells and gametogenic germ cells. This study aims to elucidate genetic interactions between two members, Pramex1 and Pramel1, in the mouse Prame family during gametogenesis using a gene knockout approach.

RESULT

Single-gene knockout (sKO) of either Pramex1 or Pramel1 resulted in approximately 7% of abnormal seminiferous tubules, characterized by a Sertoli-cell only (SCO) phenotype, impacting sperm count and fecundity significantly. Remarkably, sKO female mice displayed normal reproductive functions. In contrast, Pramex1/Pramel1 double knockout (dKO) mice exhibited reduced fecundity in both sexes. In dKO females, ovarian primary follicle count decreased by 50% compared to sKO and WT mice, correlating with a 50% fecundity decrease. This suggested compensatory roles during oogenesis in Pramex1 or Pramel1 sKO females. Conversely, dKO males showed an 18% frequency of SCO tubules, increased apoptotic germ cells, and decreased undifferentiated spermatogonia compared to sKO and WT testes. Western blot analysis with PRAMEX1- or PRAMEL1-specific antibodies on sKO testes revealed compensatory upregulation of each protein (30-50%) in response to the other gene's deletion. Double KO males exhibited more severe defects in sperm count and litter size, surpassing Pramex1 and Pramel1 sKO accumulative effects, indicating a synergistic enhancement interaction during spermatogenesis. Additional experiments administering trans-retinoic acid (RA) and its inhibitor (WIN18,446) in sKO, dKO, and WT mice suggested that PRAMEX1 and PRAMEL1 synergistically repress the RA signaling pathway during spermatogenesis.

CONCLUSION

Data from sKO and dKO mice unveil a synergistic interaction via the RA signaling pathway between Pramex1 and Pramel1 genes during gametogenesis. This discovery sets the stage for investigating interactions among other members within the Prame family, advancing our understanding of multi-copy gene families involved in germ cell formation and function.

Genetic conflict and the origin of multigene families: implications for sex chromosome evolution

Sex chromosomes are havens for intragenomic conflicts. The absence of recombination between sex chromosomes creates the opportunity for the evolution of segregation distorters: selfish genetic elements that hijack different aspects of an individual's reproduction to increase their own transmission. Biased (non-Mendelian) segregation, however, often occurs at a detriment to their host's fitness, and therefore can trigger evolutionary arms races that can have major consequences for genome structure and regulation, gametogenesis, reproductive strategies and even speciation. Here, we review an emerging feature from comparative genomic and sex chromosome evolution studies suggesting that meiotic drive is pervasive: the recurrent evolution of paralogous sex-linked gene families. Sex chromosomes of several species independently acquire and co-amplify rapidly evolving gene families with spermatogenesis-related functions, consistent with a history of intragenomic conflict over transmission. We discuss Y chromosome features that might contribute to the tempo and mode of evolution of X/Y co-amplified gene families, as well as their implications for the evolution of complexity in the genome. Finally, we propose a framework that explores the conditions that might allow for recurrent bouts of fixation of drivers and suppressors, in a dosage-sensitive fashion, and therefore the co-amplification of multigene families on sex chromosomes.

Versatile generation of precise gene edits in bovines using SEGCPN

BACKGROUND

Gene knockout and knock-in have been widely performed in large farm animals based on genome editing systems. However, many types of precise gene editing, including targeted deletion, gene tagging, and large gene fragment replacement, remain a challenge in large farm animals.

RESULTS

Here, we established versatile self-excising gene-targeting technology in combination with programmable nucleases (SEGCPN) to efficiently generate various types of precise gene editing in bovine. First, we used this versatile method to successfully generate bovine embryos with point mutations and 11-bp deletions at the MSTN locus. Second, we successfully generated bulls with EGFP labeling at the SRY locus. Finally, we successfully generated humanized cows in which the endogenous 18-kb α-casein gene was replaced with a 2.6-kb human α-lactalbumin gene.

CONCLUSIONS

In summary, our new SEGCPN method offers unlimited possibilities for various types of precise gene editing in large animals for application both in agriculture and disease models.

Transcriptomic Analysis of the Developing Testis and Spermatogenesis in Qianbei Ma Goats

Reproductive competence in male mammals depends on testicular function. Testicular development and spermatogenesis in goats involve highly complex physiological processes. In this study, six testes were, respectively, obtained from each age group, immature (1 month), sexually mature (6 months) and physically mature (12 months old) Qianbei Ma goats. RNA-Seq was performed to assess testicular mRNA expression in Qianbei Ma goats at different developmental stages. Totally, 18 libraries were constructed to screen genes and pathways involved in testis development and spermatogenesis. Totally, 9724 upregulated and 4153 downregulated DEGs were found between immature (I) and sexually mature (S) samples; 7 upregulated and 3 downregulated DEGs were found between sexually mature (S) and physically mature (P) samples, and about 4% of the DEGs underwent alternative splicing events between I and S. Select genes were assessed by qRT-PCR, corroborating RNA-Seq findings. The detected genes have key roles in multiple developmental stages of goat testicular development and spermatogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine differentially expressed genes (DEGs). GO analysis revealed DEGs between S and P contributed to "reproduction process", "channel activity" and "cell periphery part" between I and S, and in "ion transport process", "channel activity" and "transporter complex part". KEGG analysis suggested the involvement of "glycerolipid metabolism", "steroid hormone biosynthesis" and "MAPK signaling pathway" in testis development and spermatogenesis. Genes including IGF1, TGFB1, TGFBR1 and EGFR may control the development of the testis from immature to sexually mature, which might be important candidate genes for the development of goat testis. The current study provides novel insights into goat testicular development and spermatogenesis.

Sexing of cattle embryos using RNA-sequencing data or polymerase chain reaction based on a complete sequence of cattle chromosome Y

When necessary, RNA-sequencing data or polymerase chain reaction (PCR) assays can be used to determine the presence of the chromosome Y (ChrY) in samples. This information allows for biological variation due to sexual dimorphism to be studied. A prime example is when researchers conduct RNA-sequencing of single embryos, or conceptuses, prior to the development of gonads. A recent publication of a complete sequence of the ChrY has removed limitations for the development of these procedures in cattle, otherwise imposed by the absence of a ChrY in the reference genome. Using the sequence of the cattle ChrY and transcriptome data, we conducted a systematic search for genes in the ChrY that are exclusively expressed in male tissues. The genes ENSBIXG00000029763, ENSBIXG00000029774, ENSBIXG00000029788, and ENSBIXG00000029892 were consistently expressed across male tissues and lowly expressed or absent in female samples. We observed that the cumulative values of counts per million were 2688-fold greater in males than the equivalent values in female samples. Thus, we deemed these genes suitable for the sexing of samples using RNA-sequencing data. We successfully used this set of genes to infer the sex of 22 cattle blastocysts (8 females and 14 males). Additionally, the completed sequence of the cattle ChrY has segments in the male-specific region that are not repeated. We designed a pair of oligonucleotides that targets one of these non-repeated regions in the male-specific sequence of the ChrY. Using this pair of oligonucleotides, in a multiplexed PCR assay with oligonucleotides that anneal to an autosome chromosome, we accurately identified the sex of cattle blastocysts. We developed efficient procedures for the sexing of samples in cattle using either transcriptome data or their DNA. The procedures using RNA-sequencing will greatly benefit researchers who work with samples limited in cell numbers which are only sufficient to produce transcriptome data. The oligonucleotides used for the accurate sexing of samples using PCR are transferable to other cattle tissue samples.

Role of the bovine PRAMEY protein in sperm function during in vitro fertilization (IVF)

Preferentially expressed antigen in melanoma (PRAME) is a cancer/testis antigen (CTA) that is predominantly expressed in normal male gonad tissues and a variety of tumors. PRAME proteins are present in the acrosome and sperm tail, but their role in sperm function is unknown. The objective of this study was to examine the function of the bovine Y-linked PRAME (PRAMEY) during spermatozoal capacitation, the acrosome reaction (AR), and fertilization. Freshly ejaculated spermatozoa were induced to capacitate and undergo AR in vitro. Western blotting results revealed a decrease in the PRAMEY protein in capacitated spermatozoa, and the release of the PRAMEY protein from the acrosome during the AR, suggesting its involvement in sperm capacitation and AR. IVF was performed using in vitro matured bovine oocytes and cauda epididymal spermatozoa either treated with PRAMEY antibody, rabbit IgG, or DPBS. Sperm-egg binding and early embryos were examined at 6 and 45 h post IVF, respectively. The number of spermatozoa that bound per oocyte was nearly two-fold greater in the PRAMEY antibody treatment group (34.4) when compared to both the rabbit IgG (17.6) and DPBS (18.1) controls (P < 0.01). Polyspermy rate in the antibody-treated group (18.9%) was three-fold greater than the rabbit IgG control (6.0%) (P < 0.01). The results indicate that PRAMEY may play a role in anti-polyspermy defense. This study thus provides the initial evidence for the involvement of the PRAME protein family in sperm function and fertilization.

Identification of Y-SNPs within ovine MSY region and their association with testicular size

The screening of genomic variations within the male-specific region of the mammalian Y chromosome (MSY) is one of the most effective ways to investigate paternal evolutionary history and identify molecular markers related to male fertility. The current study was to identify single nucleotide polymorphisms (SNPs) within single-copy genes of the ovine MSY, and confirm whether they are associated with testicular size. A total of 21 Y-specific gene fragments were successfully amplified to screen Y-SNPs in 956 rams across nine sheep breeds. Three Y-SNPs, including SRY16: g.88 A > G in South African Mutton Merino sheep, ZFY16: g.146 C > T in Suffolk and South African Mutton Merino sheep, and EIF2S3Y2: g.77 C > G in Hu and Tan sheep, were identified using DNA-pooled sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) methods. The investigation of the global distribution for three Y-SNPs showed that the C allele of ZFY16: g.146 C > T co-segregated with haplogroup y-HC, and the C/G allele of EIF2S3Y2: g.77 C > G co-segregated with haplogroup y-HA/y-HB1 in Hu sheep according to data mining from a previous study. In addition, association analysis revealed that ZFY16: g.146 C > T had a significant effect on yearling scrotal circumference in Suffolk sheep, and EIF2S3Y2: g.77 C > G was significantly associated with testicular and epididymis weight in Hu sheep (P ≤ 0.05). The current study concluded that Y-SNPs were associated with testicular size in specific sheep, which provides valuable candidate makers for selecting elite rams at an early age.

Association between the morphokinetics of in-vitro-derived bovine embryos and the transcriptomic profile of the derived blastocysts

The time-lapse system is a non-invasive method that enables a continuous evaluation through embryo development. Here, we examined the association between the morphokinetics of the developing embryo and the transcriptomic profile of the formed blastocysts. Bovine oocytes were matured and fertilized in vitro; then, the putative zygotes were cultured in an incubator equipped with a time-lapse system. Based on the first-cleavage pattern, embryos were categorized as normal or abnormal (68.5±2.2 and 31.6±2.3%, respectively; P<0.001). A cleaved embryo was defined as normal when it first cleaved into two equal blastomeres; it was classified as synchronous or asynchronous according to its subsequent cleavages. An abnormal pattern was defined as direct, unequal, or reverse cleavage. Direct cleavage was classified as division from one cell directly into three or more blastomeres; unequal cleavage was classified as division that resulted in asymmetrically sized blastomeres; and reverse cleavage of the first division was classified as reduced number of blastomeres from two to one. Of the normally cleaving embryos, 60.2±3.1% underwent synchronous cleavage into 4, 8, and 16 blastomeres, and 39.7±3.1% cleaved asynchronously (P<0.001). The blastocyte formation rate was lower for the synchronously vs. the asynchronously cleaved embryos (P<0.03). The abnormally cleaved embryos showed low competence to develop to blastocysts, relative to the normally cleaved embryos (P<0.001). Microarray analysis revealed 895 and 643 differentially expressed genes in blastocysts that developed from synchronously and asynchronously cleaved embryos, respectively, relative to those that developed from directly cleaved embryos. The genes were related to the cell cycle, cell differentiation, metabolism, and apoptosis. About 180 differentially expressed genes were found between the synchronously vs. the asynchronously cleaved embryos, related to metabolism and the apoptosis mechanism. We provide the first evidence indicating that an embryo's morphokinetics is associated with the transcriptome profile of the derived blastocyst, which might be practically relevant for the embryo transfer program.

PRAMEY: A Bovid-Specific Y-Chromosome Multicopy Gene Is Highly Related to Postnatal Testicular Growth in Hu Sheep

PRAMEY (preferentially expressed antigen in melanoma, Y-linked) belongs to the cancer-testis antigens (CTAs) gene family and is predominantly expressed in testis, playing important roles in spermatogenesis and testicular development. This study cloned the full-length cDNA sequence of ovine PRAMEY using the rapid amplification of cDNA ends (RACE) method and analyzed the expression profile and copy number variation (CNV) of PRAMEY using quantitative real-time PCR (qPCR). The results revealed that the PRAMEY cDNA was 2099 bp in length with an open reading frame (ORF) of 1536 bp encoding 511 amino acids. PRAMEY was predominantly expressed in the testis and significantly upregulated during postnatal testicular development. The median copy number (MCN) of PRAMEY was 4, varying from 2 to 25 in 710 rams across eight sheep breeds. There was no significant correlation between the CNV of PRAMEY and testicular size, while a significant positive correlation was observed between the mRNA expression and testicular size in Hu sheep. The current study suggests that the expression levels of PRAMEY were closely associated with testicular size, indicating that PRAMEY may play an important role in testicular growth.

Copy number variation of ZNF280BY across eight sheep breeds and its association with testicular size of Hu sheep

ZNF280BY, a bovid-specific Y chromosome gene, was firstly found to be highly expressed in bovine testis, indicating it may play important roles in testicular development and male fertility. In this study, we firstly cloned the full-length cDNA of ovine ZNF280BY containing 1993 bp, and with a 1632 bp open reading frame. ZNF280BY was predominantly expressed in the testis, and its expression level was significantly higher in large testis than in small testis in Hu sheep at 6 months of age. In addition, the expression level of ZNF280BY significantly increased during testicular development, showing the highest expression level at 12 months of age. ZNF280BY showed copy number variation (CNV) in 723 rams from eight sheep breeds, ranging from 17 to 514 copies, with a median copy number of 188. Pearson correlation analysis showed that the CNV of ZNF280BY was negatively correlated with testis size in Hu sheep. Furthermore, its mRNA expression level in testis had no significant correlation with the CNV but was significantly correlated with testis size. This study concluded that the expression of ZNF280BY was closely related to testicular development, and the CNV of ZNF280BY could be used as an important genetic marker to evaluate the ram reproductive capacity at an early stage in Hu sheep.

The Dynamic of PRAMEY Isoforms in Testis and Epididymis Suggests Their Involvement in Spermatozoa Maturation

The preferentially expressed antigen in melanoma, Y-linked (PRAMEY) is a cancer/testis antigen expressed predominantly in bovine spermatogenic cells, playing an important role in germ cell formation. To better understand PRAMEY’s function during spermatogenesis, we studied the dynamics of PRAMEY isoforms by Western blotting (WB) with PRAMEY-specific antibodies. The PRAMEY protein was assessed in the bovine testicular and epididymal spermatozoa, fluid and tissues, and as well as in ejaculated semen. The protein was further examined, at a subcellular level in sperm head and tail, as well as in the subcellular components, including the cytosol, nucleus, membrane, and mitochondria. RNA expression of PRAMEY was also evaluated in testis and epididymal tissues. Our WB results confirmed the previously reported four isoforms of PRAMEY (58, 30, 26, and 13 kDa) in the bovine testis and spermatozoa. We found that testicular spermatozoa expressed the 58 and 30 kDa isoforms. As spermatozoa migrated to the epididymis, they expressed two additional isoforms, 26 and 13 kDa. Similarly, the 58 and 30 kDa isoforms were detected only in the testis fluid, while all four isoforms were detected in fluid from the cauda epididymis. Tissue evaluation indicated a significantly higher expression of the 58 and 13 kDa isoforms in the cauda tissue when compared to both the testis and caput tissue (p < 0.05). These results indicated that testis samples (spermatozoa, fluid, and tissue) expressed predominantly the 58 and 30 kDa PRAMEY isoforms, suggesting their involvement in spermatogenesis. In contrast, the 26 kDa isoform was specific to epididymal sperm and the 13 kDa isoform was marked in samples derived from the cauda epididymis, suggesting their involvement in sperm maturation. Results from the sperm head and tail experiments indicated that the 13 kDa isoform increased 4-fold in sperm tails from caput to cauda, suggesting this isoform may have a significant role in tail function. Additionally, the 13 kDa isoform increased significantly (p < 0.05) in the cytosol during epididymal passage and tended to increase in other subcellular components. The expression of PRAMEY in the sperm subcellular components during epididymal maturation suggests the involvement of PRAMEY, especially the 13 kDa isoform, in sperm motility.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

The selection of testis size can improve the reproductive capacity of livestock used for artificial insemination and has been considered as an important strategy for accelerating the breeding process. Although much work has been done to investigate the mechanisms of testis development in various species, there is little information available in regard to the differences in transcriptomic profiling of sheep testes at different developmental stages. In this work, we aimed to identify differentially expressed genes (DEGs) by RNA-Seq in sheep during different growth stages, including 0 month old (infant, M0), 3 months old (puberty, M3), 6 months old (sexual maturity, M6) and 12 months old (body maturity, M12). A total of 4,606 (2,381 up and 2,225 down), 7,500 (4,368 up and 3,132 down), 15 (8 up and seven down) DEGs were identified in M3_vs_M0, M6_vs_M3, and M12_vs_M6 comparison, respectively. Of which, a number of genes were continuously up-regulated and down-regulated with testicular development, including ODF3, ZPBP1, PKDREJ, MYBL1, PDGFA, IGF1, LH, INSL3, VIM, AMH, INHBA, COL1A1, COL1A2, and INHA. GO analysis illustrated that DEGs were mainly involved in testis development and spermatogenesis. KEGG analysis identified several important pathways and verified several reproduction-associated DEGs such as COL1A1, COL1A2, PDGFA, and IGF1. In addition, two gene modules highly associated with testis development and core genes with testis size were identified using weighted gene co-expression network analysis (WGCNA), including hub genes positively associated with testis size such as RANBP9, DNAH17, SPATA4, CIB4 and SPEM1, and those negatively associated with testis size such as CD81, CSK, PDGFA, VIM, and INHBA. This study comprehensively identified key genes related to sheep testicular development, which may provide potential insights for understanding male fertility and better guide in animal breeding.

A mutation modulating DDX3Y gene expression cosegregates with the major Y-chromosomal haplogroups and with testis size in Hu sheep

Variations in the Y-chromosome are usually correlated with male-specific traits. However, this condition has been described only sporadically, even in human genetics. The present study was conducted to clone the full-length gene sequence of ovine DEAD-box helicase 3, Y-linked (DDX3Y), and investigate the effect of the expression and variation of DDX3Y on the reproductive traits of Hu sheep. Consequently, we identified the full coding sequence and genomic sequence of ovine DDX3Y. Quantitative PCR (qPCR) analysis showed that ovine DDX3Y was highly expressed in testis, and the expression level increased during testicular development. Furthermore, individuals with larger testis at 6 months expressed significantly more DDX3Y mRNA in the testis than individuals with smaller testis. Notably, a novel SNP (g. 12657 C>A) in the 3' untranslated region was identified in Hu sheep and Tan sheep according to the investigation of the full DDX3Y genomic sequence of 1069 individuals from nine sheep breeds. Association analysis revealed that the SNP was significantly related to testis size in Hu sheep. Meanwhile, Hu rams with the derived C allele showed significantly higher expression levels of DDX3Y in testis than those with the ancestral A allele. In addition, data mining in a previous study showed that the C allele cosegregated with the globally major Y-chromosomal haplogroups y-HA and y-HC, and the A allele is found in all rams with haplogroups y-HB1, y-HB2 and y-HD. This study suggests that the association of the Y-chromosomal haplogroups with testis size in Hu sheep can be extrapolated to the sheep population worldwide.

A near complete genome for goat genetic and genomic research

Background

Goat, one of the first domesticated livestock, is a worldwide important species both culturally and economically. The current goat reference genome, known as ARS1, is reported as the first nonhuman genome assembly using 69× PacBio sequencing. However, ARS1 suffers from incomplete X chromosome and highly fragmented Y chromosome scaffolds.

Results

Here, we present a very high-quality de novo genome assembly, Saanen_v1, from a male Saanen dairy goat, with the first goat Y chromosome scaffold based on 117× PacBio long-read sequencing and 118× Hi-C data. Saanen_v1 displays a high level of completeness thanks to the presence of centromeric and telomeric repeats at the proximal and distal ends of two-thirds of the autosomes, and a much reduced number of gaps (169 vs. 773). The completeness and accuracy of the Saanen_v1 genome assembly are also evidenced by more assembled sequences on the chromosomes (2.63 Gb for Saanen_v1 vs. 2.58 Gb for ARS1), a slightly increased mapping ratio for transcriptomic data, and more genes anchored to chromosomes. The eight putative large assembly errors (1 to ~ 7 Mb each) found in ARS1 were amended, and for the first time, the substitution rate of this ruminant Y chromosome was estimated. Furthermore, sequence improvement in Saanen_v1, compared with ARS1, enables us to assign the likely correct positions for 4.4% of the single nucleotide polymorphism (SNP) probes in the widely used GoatSNP50 chip.

Conclusions

The updated goat genome assembly including both sex chromosomes (X and Y) and the autosomes with high-resolution quality will serve as a valuable resource for goat genetic research and applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00668-5.

Molecular fingerprint of female bovine embryos produced in vitro with high competence to establish and maintain pregnancy†

The objective was to identify the transcriptomic profile of in vitro-derived embryos with high competence to establish and maintain gestation. Embryos produced with X-sorted sperm were cultured from day 5 to day 7 in serum-free medium containing 10 ng/ml recombinant bovine colony-stimulating factor 2 (CSF2) or vehicle. The CSF2 was administered because this molecule can increase blastocyst competence for survival after embryo transfer. Blastocysts were harvested on day 7 of culture and manually bisected. One demi-embryo from a single blastocyst was transferred into a synchronized recipient and the other half was used for RNA-seq analysis. Using P < 0.01 and a fold change >2-fold or <0.5 fold as cutoffs, there were 617 differentially expressed genes (DEG) between embryos that survived to day 30 of gestation vs those that did not, 470 DEG between embryos that survived to day 60 and those that did not, 432 DEG between embryos that maintained pregnancy from day 30 to day 60 vs those where pregnancy failed after day 30, and 635 DEG regulated by CSF2. Pathways and ontologies in which DEG were overrepresented included many related to cellular responses to stress and cell survival. It was concluded that gene expression in the blastocyst is different between embryos that are competent to establish and maintain pregnancy vs those that are not. The relationship between expression of genes related to cell stress and subsequent embryonic survival probably reflects cellular perturbations caused by embryonic development taking place in the artificial environment associated with cell culture.

Gene expression profiles of bovine genital ridges during sex determination and early differentiation of the gonads†

Most current knowledge of sex determination in mammals has emerged from mouse and human studies. To investigate the molecular regulation of the sex determination process in cattle, we used an RNA sequencing strategy to analyze the transcriptome landscape of male and female bovine fetal gonads collected in vivo at key developmental stages: before, during, and after SRY gene activation on fetal days D35 (bipotential gonad formation), D39 (peak SRY expression), and D43 (early gonad differentiation). Differentially expressed genes (DEGs) were identified in male vs. female germinal ridges and among group genes showing similar expression profiles during the three periods. There were 143, 96, and 658 DEG between males and female fetuses at D35, D39, and D43, respectively. On D35, genes upregulated in females were enriched in translation, nuclear export, RNA localization, and mRNA splicing events, whereas those upregulated in males were enriched in cell proliferation regulation and male sex determination terms. In time-course experiments, 767 DEGs in males and 545 DEGs in females were identified between D35 vs. D39, and 3157 DEGs in males and 2008 in females were identified between D39 vs. D43. Results highlight unique aspects of sex determination in cattle, such as the expression of several Y chromosome genes (absent in mice and humans) before SRY expression and an abrupt increase in the nuclear expression of SOX10 (instead of SOX9 expression in the Sertoli cell cytoplasm as observed in mice) during male determination and early differentiation.

Efficient TALEN-mediated gene knockin at the bovine Y chromosome and generation of a sex-reversal bovine

Functional elucidation of bovine Y-chromosome genes requires available genome editing technologies. Meanwhile, it has yet to be proven whether the bovine Sry gene is the main or single factor involved in the development of the male phenotype in bovine. Here, we efficiently knocked out four Y-linked genes (Sry, ZFY, DDX3Y, and EIF2S3Y) in bovine fetal fibroblasts (BFFs) with transcription activator-like effector nucleases (TALENs) individually. Furthermore, we used TALEN-mediated gene knockin at the Sry gene and generated a sex-reversal bovine by somatic cell nuclear transfer (SCNT). The resulting bovine had only one ovary and was sterile. We demonstrate, for the first time, that the Sry gene is an important sex-determining gene in bovine. Our method lays a solid foundation for detecting the biology of the bovine Y chromosome, as it may provide an alternative biological model system for the study of mammalian sex determination, and new methods for the practical application in agricultural, especially for sex predetermination.

The PRAME family of cancer testis antigens is essential for germline development and gametogenesis†

Preferentially expressed antigen in melanoma (PRAME) belongs to a group of cancer/testis antigens that are predominately expressed in the testis and a variety of tumors, and are involved in immunity and reproduction. Much of the attention on PRAME has centered on cancer biology as PRAME is a prognostic biomarker for a wide range of cancers and a potential immunotherapeutic target. Less information is available about the PRAME family's function (s) during gametogenesis and in the overall reproduction process. Here, we review the current knowledge of the PRAME gene family and its function in germline development and gametogenesis. Members of the PRAME family are leucine rich repeat proteins, localized in nucleus and cytoplasm, with multifaceted roles in germ cells. As transcriptional regulators, the PRAME family proteins are involved in germline development, particularly in the maintenance of embryonic stem cell pluripotency, development of primordial germ cells, and differentiation/proliferation of spermatogenic and oogenic cells. The PRAME family proteins are also enriched in cytoplasmic organelles, such as rough endoplasmic reticulum, Golgi vesicle, germinal granules, centrioles, and play a role in the formation of the acrosome and sperm tail during spermiogenesis. The PRAME gene family remains transcriptionally active in the germline throughout the entire life cycle and is essential for gametogenesis, with some members specific to either male or female germ cells, while others are involved in both male and female gametogenesis. A potential molecular mechanism that underlies the function of PRAME, and is shared by gametogenesis and oncogenesis is also discussed.

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Effect of the X chromosome in genomic evaluations of reproductive traits in beef cattle

The aim of this study was to evaluate whether there are predictive advantages for breeding values with inclusion of X chromosome genomic markers for reproductive (occurrence of early pregnancy - P16 and age at first calving - AFC) and andrological (scrotal circumference -SC) variables in beef cattle. There were 3263 genotypes of females and males evaluated. There were breeding value estimates for SC, AFC and P16 considering two scenarios: 1) only autosomal markers or 2) autosomal and X chromosome markers. To evaluate effects of inclusion of X chromosome markers on selection, responses to selection were compared including or not including genomic marker information from the X chromosome. There were greater heritability estimates for SC (0.40 and 0.31), AFC (0.11 and 0.09) and P16 (0.43 and 0.38) when analyses included, compared with not including, genomic marker information from the X chromosome. When selection is based on results from analyses that did not include information for the X chromosome, there was about a 7 % lesser mean genomic breeding value for the SC traits for selected animals. For P16, there was an approximate 4% lesser breeding value without inclusion of genomic marker information from the X chromosome, while this inclusion did not have as great an effect on the breeding value for AFC. There was an average predictive correlation of 0.79, 0.98 and 0.84 for SC, AFC and P16, respectively. These estimates indicate inclusion of the X chromosome genomic marker information in the analysis can improve prediction of genomic breeding values, especially for SC.

Sequence Transpositions Restore Genes on the Highly Degenerated W Chromosomes of Songbirds

The female-specific W chromosomes of most Neognathae birds are highly degenerated and gene-poor. Previous studies have demonstrated that the gene repertoires of the Neognathae bird W chromosomes, despite being in small numbers, are conserved across bird species, likely due to purifying selection maintaining the regulatory and dosage-sensitive genes. Here we report the discovery of DNA-based sequence duplications from the Z to the W chromosome in birds-of-paradise (Paradisaeidae, Passeriformes), through sequence transposition. The original transposition involved nine genes, but only two of them (ANXA1 and ALDH1A1) survived on the W chromosomes. Both ANXA1 and ALDH1A1 are predicted to be dosage-sensitive, and the expression of ANXA1 is restricted to ovaries in all the investigated birds. These analyses suggest the newly transposed gene onto the W chromosomes can be favored for their role in restoring dosage imbalance or through female-specific selection. After examining seven additional songbird genomes, we further identified five other transposed genes on the W chromosomes of Darwin's finches and one in the great tit, expanding the observation of the Z-to-W transpositions to a larger range of bird species, but not all transposed genes exhibit dosage-sensitivity or ovary-biased expression We demonstrate a new mechanism by which the highly degenerated W chromosomes of songbirds can acquire genes from the homologous Z chromosomes, but further functional investigations are needed to validate the evolutionary forces underlying the transpositions.

Dynamic evolution of great ape Y chromosomes

The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered.

The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.

Patterns of DNA variation between the autosomes, the X chromosome and the Y chromosome in Bos taurus genome

The new ARS-UCD1.2 assembly of the bovine genome has considerable improvements over the previous assembly and thus more accurate identification of patterns of genetic variation can be achieved with it. We explored differences in genetic variation between autosomes, the X chromosome, and the Y chromosome. In particular, variant densities, annotations, lengths (only for InDels), nucleotide divergence, and Tajima’s D statistics between chromosomes were considered. Whole-genome DNA sequences of 217 individuals representing different cattle breeds were examined. The analysis included the alignment to the new reference genome and variant identification. 23,655,295 SNPs and 3,758,781 InDels were detected. In contrast to autosomes, both sex chromosomes had negative values of Tajima’s D and lower nucleotide divergence. That implies a correlation between nucleotide diversity and recombination rate, which is obviously reduced for sex chromosomes. Moreover, the accumulation of nonsynonymous mutations on the Y chromosome could be associated with loss of recombination. Also, the relatively lower effective population size for sex chromosomes leads to a lower expected density of variants.

Genes associated with survival of female bovine blastocysts produced in vivo

The objective was to characterize the transcriptome profile of in vivo-derived female embryos competent to establish and maintain gestation. Blastocysts from superovulated heifers were bisected to generate two demi-embryos. One demi-embryo was transferred into a synchronized recipient and the other part was used for RNA-seq analysis. Data on transcript abundance was analyzed for 4 demi-embryos that established and maintained pregnancy to day 60 (designated as PP) and 3 that did not result in a pregnancy at day 30 (designated as NP). Using a false discovery rate of P < 0.10 as cutoff, a total of 155 genes were differentially expressed between PP and NP embryos, of which 73 genes were upregulated and 82 genes were downregulated in the PP group. The functional cluster with the greatest enrichment score for embryos that survived, representing 28 genes (48% of the annotated genes), was related to membrane proteins, particularly those related to olfaction and neural development and function. The functional cluster with the greatest enrichment score for downregulated genes in embryos that survived included terms related to oxidative phosphorylation, mitochondrial function, and transmembrane proteins. In conclusion, competence of in vivo-derived female bovine embryos to survive after transfer is associated with increased expression of genes encoding transmembrane proteins, perhaps indicative of differentiation of the inner cell mass to epiblast, and decreased expression of genes involved in oxidative phosphorylation, perhaps indicative of reduced metabolic activity.

Highly efficient synchronization of sheep skin fibroblasts at G2/M phase and isolation of sheep Y chromosomes by flow cytometric sorting

At present, based on whole genome sequencing, sequences and genes annotation of the sheep (Ovis aries) Y chromosome are still absent. The isolation of Y chromosomes followed by sequencing has been approved as an effective approach to analyze this complex chromosome in other species. In this study, we established a highly efficient synchronization method for G2/M phase of sheep fibroblasts, which was successfully applied to flow-sorting chromosomes of sheep, with a focus on isolation and sequencing of the ovine Y chromosome. The isolated (~80,000) Y chromosomes were verified by fluorescence quantitative real-time polymerase chain reaction, further confirmed by fluorescence in situ hybridization, and amplified by the MALBAC method before next-generation sequencing. The sequence results indicated that 68.90% of reads were Y chromosome-related sequences as they are homologous to the bovine Y chromosome. The remaining 31.1% of reads were aligned to the sheep reference genome, including 13.57% reads to chromosome X and 6.68% to chromosome 17. Importantly, the paired-end reads that are properly aligned to the bovine Y sequence assembly accounted for 46.49%, indicating the success in the ovine Y chromosome isolation and the high quality of the Y chromosome sequences. This study not only set up a foundation for future sequencing, assembly and annotation of the ovine Y chromosome, but also provide a validated approach to overcoming difficulties in sequencing Y chromosome in other mammalian species.

Gene mapping and genomic prediction of bull fertility using sex chromosome markers

Service sire has been recognized as an important factor affecting dairy herd fertility. Our group has reported promising results on gene mapping and genomic prediction of dairy bull fertility using autosomal SNP markers. Little is known, however, about the genetic contribution of sex chromosomes, which are enriched in genes related to sexual development and reproduction. As such, the main goal of this study was to investigate the effect of SNP markers on X and Y chromosomes (BTAX and BTAY, respectively) on sire conception rate (SCR) in US Holstein bulls. The analysis included a total of 5,014 bulls with SCR records and genotypes for roughly 291k SNP located on the autosomes, 1.5k SNP located on the pseudoautosomal region (PAR), 13.7k BTAX-specific SNP, and 24 BTAY-specific SNP. We first performed genomic scans of the sex chromosomes, and then we evaluated the genomic prediction of SCR including BTAX SNP markers in the predictive models. Two markers located on PAR and 3 markers located on the X-specific region showed significant associations with sire fertility. Interestingly, these regions harbor genes, such as FAM9B, TBL1X, and PIH1D3, that are directly implicated in testosterone concentration, spermatogenesis, and sperm motility. On the other hand, BTAY showed very low genetic variability, and none of the segregating markers were associated with SCR. Notably, model predictive ability was largely improved by including BTAX markers. Indeed, the combination of autosomal with BTAX SNP delivered predictive correlations around 0.343, representing an increase in accuracy of about 7.5% compared with the standard whole autosomal genome approach. Overall, this study provides evidence of the importance of both PAR and X-specific regions in male fertility in dairy cattle. These findings may help to improve conception rates in dairy herds through accurate genome-guided decisions on bull fertility.

Deletion of the mouse X-linked Prame gene causes germ cell reduction in spermatogenesis

Preferentially expressed antigen in melanoma (PRAME) is cancer/testis antigen and a transcriptional repressor, inhibiting the signaling of retinoic acid through the retinoic acid receptor (RAR) for promoting cell proliferation and preventing cell apoptosis in cancer cells. The role of PRAME in testis and germline is unknown. We report here the generation and characterization of an X-linked Prame conditional knockout (cKO) mouse. Although fertile, the testis size (p < .01) and sperm count (p < .05) of the Prame cKO mice were significantly reduced by 12% at 4 months of age compared with the Prame floxed mice. Histological, immunofluorescence with germ cell-specific markers and terminal deoxynucleotidyl transferase dUTP nick end labeling analyses of testis cross-sections at postnatal day 7 (P7), P14, P21, P35, P120, and P365 indicated a significant increase in apoptotic germ cells at P7 and P14 and an increase in abnormal seminiferous tubules at P21 and P35. Germ cells were gradually lost resulting in two different phenotypes in the Prame cKO testes: Sertoli-cell-only for some of the affected tubules in young mice (at P35) and germ cell arrest at spermatogonia stage for other affected tubules in mature mice. Both phenotypes were a consequence of disruption in RAR signaling pathway by the depletion of Prame at a different time point during the first and subsequent rounds of spermatogenesis. The results suggest that Prame plays a minor, but important role in spermatogenesis and different paralogs in the Prame gene family may be functionally and partially redundant.

Copy Number Variations of Four Y-Linked Genes in Swamp Buffaloes

Simple Summary

The amplification of the male-specific region of the Y chromosome was a unique phenomenon during mammalian sex chromosome evolution. The Y-linked copy number variations of many species have been confirmed. However, the Y-linked copy number variations (CNVs) in water buffalo are still unknown. In this study, we investigated the copy number variations of four Y-linked genes (SRY, UTY, DBY, and OFD1Y) in buffalo. Our results showed that UTY was a single-copy gene in buffalo, while DBY, OFD1Y, and SRY exhibited copy number variations in buffalo.

Abstract

Copy number variation (CNV), a significant source of genetic diversity in the mammalian Y chromosome, is associated with the development of many complex phenotypes, such as spermatogenesis and male fertility. The contribution of Y-linked CNVs has been studied in various species, however, water buffalo has not been explored in this area and the genetic information still remains unknown. The aim of the current study was to investigate the CNVs of four Y-linked genes, including, sex determining Region of Y-Chromosome (SRY), ubiquitously transcribed tetratricopeptide repeat gene protein on the chromosome Y (UTY), DEAD-box helicase 3 Y-linked (DDX3Y, also known as DBY), and oral-facial-digital syndrome 1 Y-linked (OFD1Y) in 254 swamp buffaloes from 15 populations distributed across China, Vietnam, and Laos using quantitative real-time PCR (qPCR). Our results revealed the prevalence of a single-copy UTY gene in buffaloes. The DBY and OFD1Y represented CNVs among and within different buffalo breeds. The SRY showed CNVs only in Vietnamese and Laotian buffaloes. In conclusion, this study indicated that DBY, OFD1Y, and SRY showed CNVs, while the UTY was a single-copy gene in swamp buffaloes.

New insights into mammalian sex chromosome structure and evolution using high-quality sequences from bovine X and Y chromosomes

BACKGROUND

Mammalian X chromosomes are mainly euchromatic with a similar size and structure among species whereas Y chromosomes are smaller, have undergone substantial evolutionary changes and accumulated male specific genes and genes involved in sex determination. The pseudoautosomal region (PAR) is conserved on the X and Y and pair during meiosis. The structure, evolution and function of mammalian sex chromosomes, particularly the Y chromsome, is still poorly understood because few species have high quality sex chromosome assemblies.

RESULTS

Here we report the first bovine sex chromosome assemblies that include the complete PAR spanning 6.84 Mb and three Y chromosome X-degenerate (X-d) regions. The PAR comprises 31 genes, including genes that are missing from the X chromosome in current cattle, sheep and goat reference genomes. Twenty-nine PAR genes are single-copy genes and two are multi-copy gene families, OBP, which has 3 copies and BDA20, which has 4 copies. The Y chromosome X-d1, 2a and 2b regions contain 11, 2 and 2 gametologs, respectively.

CONCLUSIONS

The ruminant PAR comprises 31 genes and is similar to the PAR of pig and dog but extends further than those of human and horse. Differences in the pseudoautosomal boundaries are consistent with evolutionary divergence times. A bovidae-specific expansion of members of the lipocalin gene family in the PAR reported here, may affect immune-modulation and anti-inflammatory responses in ruminants. Comparison of the X-d regions of Y chromosomes across species revealed that five of the X-Y gametologs, which are known to be global regulators of gene activity and candidate sexual dimorphism genes, are conserved.

Comparison of mRNA expression from Y-chromosome X-degenerate region genes in taurine cattle, yaks and interspecific hybrid bulls

The yattle (dzo) is an interspecific hybrid of the taurine cattle (Bos taurus) and the domestic yak (Bos grunniens). F₁ hybrid yattle bulls are sterile due to spermatogenic arrest and have misregulation of spermatogenesis genes in the testes. However, the expression pattern of Y chromosome-linked genes in cattle, yaks and yattle testes is still unknown. In this study, we analyzed the mRNA expression pattern of 10 genes known to be present as single copies in the X-degenerate region of the bovine male-specific region of the Y chromosome. Using male-specific primers and reverse transcription quantitative PCR, the ubiquitously transcribed tetratricopeptide repeat gene, Y-linked (UTY), oral-facial-digital syndrome 1, Y-linked (OFD1Y) and ubiquitin specific peptidase 9, Y-linked (USP9Y) genes were ubiquitously expressed and significantly more highly expressed in yattle than in cattle and yaks testes (P < 0.001). RNA binding motif protein, Y-linked (RBMY) had testes-specific expression, and eukaryotic translation initiation factor 1A, Y-linked (EIF1AY) was expressed mainly in testis, whereas yattle and cattle did not show significant differences with respect to the expression of RBMY and EIF1AY. Thus, based on the model of yattle bull sterility, the high expression of UTY, OFD1Y and USP9Y may be associated with yattle infertility.

Gene expression dynamics during the gonocyte to spermatogonia transition and spermatogenesis in the domestic yak

Background

Spermatogenesis is a cellular differentiation process that includes three major events: mitosis of spermatogonia, meiosis of spermatocytes and spermiogenesis. Steady-state spermatogenesis relies on functions of spermatogonial stem cells (SSCs). Establishing and maintaining a foundational SSC pool is essential for continued spermatogenesis in mammals. Currently, our knowledge about SSC and spermatogenesis is severely limited in domestic animals.

Results

In the present study, we examined transcriptomes of testes from domestic yaks at four different stages (3, 5, 8 and 24 months of age) and attempted to identify genes that are associated with key developmental events of spermatogenesis. Histological analyses showed that the most advanced germ cells within seminiferous tubules of testes from 3, 5, 8 and 24 months old yaks were gonocytes, spermatogonia, spermatocytes and elongated spermatids, respectively. RNA-sequencing (RNA-seq) analyses revealed that 11904, 4381 and 2459 genes were differentially expressed during the gonocyte to spermatogonia transition, the mitosis to meiosis transition and the meiosis to post-meiosis transition. Further analyses identified a list of candidate genes than may regulate these important cellular processes. CXCR4, a previously identified SSC niche factor in mouse, was one of the up-regulated genes in the 5 months old yak testis. Results of immunohistochemical staining confirmed that CXCR4 was exclusively expressed in gonocytes and a subpopulation of spermatogonia in the yak testis.

Conclusions

Together, these findings demonstrated histological changes of postnatal testis development in the domestic yak. During development of spermatogonial lineage, meiotic and haploid germ cells are supported by dynamic transcriptional regulation of gene expression. Our transcriptomic analyses provided a list of candidate genes that potentially play crucial roles in directing the establishment of SSC and spermatogenesis in yak.

Electronic supplementary material

The online version of this article (10.1186/s40104-019-0360-7) contains supplementary material, which is available to authorized users.

Synteny search identifies carnivore Y chromosome for evolution of male specific genes

The explosive accumulation of mammalian genomes has provided a valuable resource to characterize the evolution of the Y chromosome. Unexpectedly, the Y-chromosome sequence has been characterized in only a small handful of species, with the majority being model organisms. Thus, identification of Y-linked scaffolds from unordered genome sequences is becoming more important. Here, we used a syntenic-based approach to generate the scaffolds of the male-specific region of the Y chromosome (MSY) from the genome sequence of 6 male carnivore species. Our results identified 14, 15, 9, 28, 14 and 11 Y-linked scaffolds in polar bears, pacific walruses, red pandas, cheetahs, ferrets and tigers, covering 1.55 Mbp, 2.62 Mbp, 964 Kb, 1.75 Mb, 2.17 Mbp and 1.84 Mb MSY, respectively. All the candidate Y-linked scaffolds in 3 selected species (red pandas, polar bears and tigers) were successfully verified using polymerase chain reaction. We re-annotated 8 carnivore MSYs including these 6 Y-linked scaffolds and domestic dog and cat MSY; a total of 11 orthologous genes conserved in at least 7 of the 8 carnivores were identified. These 11 Y-linked genes have significantly higher evolutionary rates compared with their X-linked counterparts, indicating less purifying selection for MSY genes. Taken together, our study shows that the approach of synteny search is a reliable and easily affordable strategy to identify Y-linked scaffolds from unordered carnivore genomes and provides a preliminary evolutionary study for carnivore MSY genes.

Copy number variation of ZNF280AY across 21 cattle breeds and its association with the reproductive traits of Holstein and Simmental bulls

The mammalian Y chromosome gene families in the ampliconic region are expressed predominantly or exclusively in the testis, and their copy number variations (CNV) are significantly associated with male reproductive traits, suggesting they have important roles in spermatogenesis and testicular development. ZNF280AY (zinc finger protein 280A, Y-linked) is a member of the zinc finger protein family and has been identified as a bovid-specific Y-chromosome gene. The current study applied a reliable quantitative real-time PCR method to estimate the CNV of ZNF280AY in 715 bulls across 21 cattle breeds and to further investigate the association of the CNV of ZNF280AY with bull reproductive traits and ZNF280AY mRNA expression levels in adult testis. The results revealed that the median copy number of ZNF280AY was 47, and the copy number varied from 11 to 154, showing significant CNV between and within the investigated cattle breeds. In addition, all 715 bulls were classified into Y1, Y2, and Y3 lineage groups based on a rapid genotyping method described previously. Pairwise comparisons indicated that bulls belonging to the Y1 lineage had a significantly lower median copy number (40) than bulls belonging to the Y2 (52) and Y3 lineages (57). Association analysis revealed that the CNV of ZNF280AY was correlated negatively with the percentage of normal sperm and sperm concentration in Holstein bulls, whereas no significant correlation was observed with ejaculation volume, total sperm count, sperm motility, postthaw motility (PTM), and scrotal circumference in Holstein and Simmental bulls. Furthermore, no correlation was observed between ZNF280AY copy number and ZNF280AY mRNA expression levels in the testis. The current study suggests that the CNV of the ZNF280AY gene family is associated with male reproductive traits and may serve as a valuable marker for early bull fertility selection in Holstein breeding programs.

Comparison of Y-chromosome-linked TSPY, TSPY2, and PRAMEY genes in Taurus cattle, yaks, and interspecific hybrid bulls

Domestic yaks (Bos grunniens) and domestic Taurus cattle (Bos taurus) are closely related. An interesting phenomenon in interspecific crossings is male sterility in the F₁ hybrid (yattle) and F₂ backcross, with no late meiotic cells or spermatids in the seminiferous tubules. The mammalian Y chromosome is crucial for spermatogenesis and male fertility. This study investigated the copy number variations and mRNA of Y-transitional region genes TSPY2 (testis specific protein, Y-linked 2 and testis-specific Y-encoded protein 3-like) and PRAMEY (preferentially expressed antigen in melanoma, Y-linked), and Y-ampliconic region genes TSPY (testis-specific Y-encoded protein 1-like), ZNF280BY (zinc finger protein 280B, Y-linked) and HSFY (heat-shock transcription factor, Y-linked) in mature testes from Taurus cattle, yaks, and yattle. Phylogenetic trees divided 33 copies of TSPY into major 2 types (TSPY-T1 and TSPY-T2), 19 copies of TSPY2 into 2 types (TSPY2-T1 and T2), and 8 copies of PRAMEY into 4 types (PRAMEY-T1 to T4). Searching by the Basic Local Alignment Search Tool of the TSPY2 coding sequences in GenBank revealed that TSPY2 was conserved in Bovidae. The TSPY2-T2 sequences were absent, whereas PRAMEY-T2 and PRAMEY-T4 were amplified on the yak Y chromosome. The average copy numbers of TSPY-T2 and ZNF280BY were significantly different between cattle and yaks. The TSPY-T2, TSPY2, PRAMEY, ZNF280BY, and HSFY genes were uniquely or predominantly expressed in testes. Reverse-transcription quantitative PCR showed that the TSPY-T2, PRAMEY-T2, HSFY, ZNF280BY, protamine 1 (PRM1), and protamine 2 (PRM2) genes were almost not expressed in yattle. The PRM1 and PRM2 genes are used as positive markers for spermatozoa. Thus, our results showed that the genomic structure of the Y-transitional and Y-ampliconic region differed between Taurus cattle and yaks. Dysregulated expression of Y-ampliconic region genes TSPY-T2, HSPY, ZNF280BY, and Y-transitional region gene PRAMEY-T2 may be associated with hybrid male sterility in yattle.

The evolutionary history of grey wolf Y chromosomes

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Mammalian Sex Chromosome Structure, Gene Content, and Function in Male Fertility

Mammalian sex chromosomes evolved from an ordinary pair of autosomes. The X chromosome is highly conserved, whereas the Y chromosome varies among species in size, structure, and gene content. Unlike autosomes that contain randomly mixed collections of genes, the sex chromosomes are enriched in testis-biased genes related to sexual development and reproduction, particularly in spermatogenesis and male fertility. This review focuses on how sex chromosome dosage compensation takes place and why meiotic sex chromosome inactivation occurs during spermatogenesis. Furthermore, the review also emphasizes how testis-biased genes are enriched on the sex chromosomes and their functions in male fertility. It is concluded that sex chromosomes are critical to sexual development and male fertility; however, our understanding of how sex chromosome genes direct sexual development and fertility has been hampered by the structural complexities of the sex chromosomes and by the multicopy nature of the testis gene families that also play a role in immunity, cancer development, and brain function.

Spermatogenesis and the Evolution of Mammalian Sex Chromosomes

Developmental constraint and sexual conflict shape the evolution of heteromorphic sex chromosomes. These contrasting forces are perhaps strongest during spermatogenesis in species with XY males. In this review, we consider how the unique regulatory environment and selective pressures of spermatogenesis interact to impact sex chromosome evolution in mammals. We explore how each developmental phase of spermatogenesis influences sex chromosome gene content, structure, and rate of molecular evolution, and how these attributes may contribute to speciation. We argue that a developmental context is fundamental to understanding sex chromosome evolution and that an evolutionary perspective can shed new light on our understanding of sperm development.

Eight Y chromosome genes show copy number variations in horses

Abstract. Copy number variations (CNVs), which represent a significant source of genetic diversity on the Y chromosome in mammals, have been shown to be associated with the development of many complex phenotypes, such as reproduction and male fertility. The occurrence of CNVs has been confirmed on the Y chromosome in horses. However, the copy numbers (CNs) of Equus caballus Y chromosome (ECAY) genes are largely unknown. To demonstrate the copy number variations of Y chromosome genes in horses, the quantitative real-time polymerase chain reaction (qPCR) method was applied to measure the CNVs of the eukaryotic translation initiation factor 1A Y (EIF1AY), equine testis-specific transcript on Y 1 (ETSTY1), equine testis-specific transcript on Y 4 (ETSTY4), equine testis-specific transcript on Y 5 (ETSTY5), equine transcript Y4 (ETY4), ubiquitin activating enzyme Y (UBE1Y), sex determining region Y (SRY), and inverted repeat 2 Y (YIR2) across 14 Chinese domestic horse breeds in this study. Our results revealed that these eight genes were multi-copy; furthermore, some of the well acknowledged single-copy genes such as SRY and EIF1AY were found to be multi-copy in this research. The median copy numbers (MCNs) varied among different breeds for the same gene. The CNVs of Y chromosome genes showed different distribution patterns among Chinese horse breeds, indicating the impact of natural selection on copy numbers. Our results will provide fundamental information for future functional studies.

Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia

Cattle domestication and the complex histories of East Asian cattle breeds warrant further investigation. Through analysing the genomes of 49 modern breeds and eight East Asian ancient samples, worldwide cattle are consistently classified into five continental groups based on Y-chromosome haplotypes and autosomal variants. We find that East Asian cattle populations are mainly composed of three distinct ancestries, including an earlier East Asian taurine ancestry that reached China at least ~3.9 kya, a later introduced Eurasian taurine ancestry, and a novel Chinese indicine ancestry that diverged from Indian indicine approximately 36.6–49.6 kya. We also report historic introgression events that helped domestic cattle from southern China and the Tibetan Plateau achieve rapid adaptation by acquiring ~2.93% and ~1.22% of their genomes from banteng and yak, respectively. Our findings provide new insights into the evolutionary history of cattle and the importance of introgression in adaptation of cattle to new environmental challenges in East Asia.

There are various indigenous cattle breeds in East Asia which have a complex history. Here, the authors analyse the genomes of 49 modern breeds and eight ancient samples and identify three distinct ancestries and multiple adaptive introgressions from other bovine species.

Review: Genomics of bull fertility

Fertility is one of the most economically important traits in both beef and dairy cattle production; however, only female fertility is typically subjected to selection. Male and female fertility have only a small positive genetic correlation which is likely due to the existence of a relatively small number of genetic variants within each breed that cause embryonic and developmental losses. Genomic tools have been developed that allow the identification of lethal recessive loci based upon marker haplotypes. Selection against haplotypes harbouring lethal alleles in conjunction with selection to improve female fertility will result in an improvement in male fertility. Genomic selection has resulted in a two to fourfold increase in the rate of genetic improvement of most dairy traits in US Holstein cattle, including female fertility. Considering the rapidly increasing rate of adoption of high-throughput single nucleotide polymorphism genotyping in both the US dairy and beef industries, genomic selection should be the most effective of all currently available approaches to improve male fertility. However, male fertility phenotypes are not routinely recorded in natural service mating systems and when artificial insemination is used, semen doses may be titrated to lower post-thaw progressively motile sperm numbers for high-merit and high-demand bulls. Standardization of sperm dosages across bull studs for semen distributed from young bulls would allow the capture of sire conception rate phenotypes for young bulls that could be used to generate predictions of genetic merit for male fertility in both males and females. These data would allow genomic selection to be implemented for male fertility in addition to female fertility within the US dairy industry. While the rate of use of artificial insemination is much lower within the US beef industry, the adoption of sexed semen in the dairy industry has allowed dairy herds to select cows from which heifer replacements are produced and cows that are used to produce terminal crossbred bull calves sired by beef breed bulls. Capture of sire conception rate phenotypes in dairy herds utilizing sexed semen will contribute data enabling genomic selection for male fertility in beef cattle breeds. As the commercial sector of the beef industry increasingly adopts fixed-time artificial insemination, sire conception rate phenotypes can be captured to facilitate the development of estimates of genetic merit for male fertility within US beef breeds.

Y-chromosome polymorphisms of the domestic Bactrian camel in China

Single-nucleotide polymorphisms (SNPs), microsatellites and copy number variation (CNV) were studied on the Y chromosome to understand the paternal origin and phylogenetic relationships for resource protection, rational development and utilization of the domestic Bactrian camel in China. Our sample set consisted of 94 Chinese domestic Bactrian camels from four regions (Inner Mongolia, Gansu, Qinghai and Xinjiang), we screened 29 Y-chromosome-specific loci for SNPs, analysed 40 bovine-derived microsatellite loci and measured CNVs of HSFY and SRY through Sanger sequencing, automated fluorescence-based microsatellite analysis and quantitative real-time PCR, respectively. A multicopy gene, SRY, was first found, and sequence variation was only detected in SRY in a screen of 29 loci in 13 DNA pools of individual camels. In addition, a TG repeat in the USP9Y gene was identified as the first polymorphic microsatellite in the camel Y chromosome, whereas microsatellite based on bovine sequences were not detected. The frequency of each allele varied among different populations. For the Nanjiang, Hexi and Alashan populations, a 243-bp allele was found. For the Sunite population, 241-bp, 243-bp and 247-bp alleles were detected, and the frequencies of these alleles were 22.2%, 44.5% and 33.3%, respectively; 241-bp and 243-bp alleles were found in other populations. Finally, CNVs in two Y-chromosomal genes were detected; CNV for HSFY and SRY ranged from 1 to 3 and from 1 to 9, respectively.

Genome-wide analysis of long non-coding RNAs and their role in postnatal porcine testis development

A comprehensive and systematic understanding of the roles of lncRNAs in the postnatal development of the pig testis has still not been achieved. In the present study, we obtained more than one billion clean reads and identified 15,528 lncRNA transcripts; these transcripts included 5032 known and 10,496 novel porcine lncRNA transcripts and corresponded to 10,041 lncRNA genes. Pairwise comparisons identified 449 known and 324 novel lncRNAs that showed differential expression patterns. GO and KEGG pathway enrichment analyses revealed that the targeted genes were involved in metabolic pathways regulating testis development and spermatogenesis, such as the TGF-beta pathway, the PI3K-Akt pathway, the Wnt/β-catenin pathway, and the AMPK pathway. Using this information, we predicted some lncRNAs and coding gene pairs were predicted that may function in testis development and spermatogenesis; these are listed in detail. This study has provided the most comprehensive catalog to date of lncRNAs in the postnatal pig testis and will aid our understanding of their functional roles in testis development and spermatogenesis.

Abundance and Characterization of Perfect Microsatellites on the Cattle Y Chromosome

Microsatellites or simple sequence repeats (SSRs) are found in most organisms and play an important role in genomic organization and function. To characterize the abundance of SSRs (1-6 base-pairs [bp]) on the cattle Y chromsome, the relative frequency and density of perfect or uninterrupted SSRs based on the published Y chromosome sequence were examined. A total of 17,273 perfect SSRs were found, with total length of 324.78 kb, indicating that approximately 0.75% of the cattle Y chromosome sequence (43.30 Mb) comprises perfect SSRs, with an average length of 18.80 bp. The relative frequency and density were 398.92 loci/Mb and 7500.62 bp/Mb, respectively. The proportions of the six classes of perfect SSRs were highly variable on the cattle Y chromosome. Mononucleotide repeats had a total number of 8073 (46.74%) and an average length of 15.45 bp, and were the most abundant SSRs class, while the percentages of di-, tetra-, tri-, penta-, and hexa-nucleotide repeats were 22.86%, 11.98%, 11.58%, 6.65%, and 0.19%, respectively. Different classes of SSRs varied in their repeat number, with the highest being 42 for dinucleotides. Results reveal that repeat categories A, AC, AT, AAC, AGC, GTTT, CTTT, ATTT, and AACTG predominate on the Y chromosome. This study provides insight into the organization of cattle Y chromosome repetitive DNA, as well as information useful for developing more polymorphic cattle Y-chromosome-specific SSRs.