Ti-Cheng Chang

Copy number variation of PRAMEY across breeds and its association with male fertility in Holstein sires

Multi-copy gene families are especially prevalent in the male-specific region (MSY) of the mammalian Y chromosome. Copy number variations (CNV) of these Y-linked gene families have been shown to affect human and animal fertility. The PRAMEY (Preferentially expressed antigen in melanoma, Y-linked) gene family is a newly identified, bovid-specific Y-linked gene family, which codes for a cancer/testis antigen that is expressed predominantly in testis and various tumors. The PRAMEY gene family is believed to play an important role in spermatogenesis and male fertility in cattle. The objective of this study was to investigate the CNV of PRAMEY within and across breeds and to determine whether CNV was associated with reproductive traits in Holstein bulls. A quantitative real-time PCR method was applied to measure the copy number of PRAMEY among 460 bulls using a Y-linked single copy gene, DDX3Y (DEAD box polypeptide 3, Y-linked), as a reference. The median copy number of PRAMEY was 13, ranging from 2 to 31. Significant variations in PRAMEY copy number were observed among 15 breeds investigated. Holstein bulls had the lowest median copy number (12), whereas Limousin bulls possessed the highest median copy number (26). Furthermore, bulls in the taurine lineage (13) had a significantly lower median copy number than those bulls in the indicine lineage (20). Association analysis revealed that PRAMEY copy number was correlated negatively with scrotal circumference (SC), relative scrotal circumference (RLSC), percentage of normal sperm (PNS), and nonreturn rate (NRR), but had no significant association with postthaw motility (PTM), incubated motility (IM), percentage of intact acrosome (PIA), sire conception rate (SCR), or relative breeding efficiency (RBE). The data from this study indicate that CNV of the PRAMEY gene family is associated with male reproductive traits and may serve as a valuable marker for sire fertility selection at an early age in cattle.

Molecular Characterization of the DDX3Y Gene and Its Homologs in Cattle

DDX3Y (also known as DBY) is a member of the DEAD box protein family, which is involved in ATP-dependent RNA unwinding, needed in a variety of cellular processes including splicing, ribosome biogenesis and RNA degradation. In the human, DDX3Y is located in the AZFa interval in the Y chromosome. Deletion of the AZFa region has been shown to disrupt spermatogenesis, causing subfertility and infertility in otherwise healthy men. Here, we report the characterization of the bovine (b) DDX3Y gene and its homologs DDX3X and PL10. We found 2 transcripts for the bDDX3Y (bDDX3Y-L and -S), which correspond to the long and short transcripts of the human DDX3Y and mouse Ddx3y gene. The 2 transcripts are identical except for a 3-bp (AGT) insertion at the position of nt 2025 and an expanded 3′UTR (nt 2155–2769) in bDDX3Y-L. The bDDX3Y-S encodes a peptide of 660 amino acids (aa), while the bDDX3Y-L encodes a peptide of 661 aa as the result of an additional serine (S) insertion at the position of aa 634. Both bDDX3Y isoforms contain the conserved DEAD-box motif. The bDDX3Y is composed of 17 exons. The homologous gene on the X chromosome, bDDX3X, is highly conserved to the Y-copy at mRNA (83%) and protein (88%) levels as well as in the genomic structure. The autosomal copy, bPL10, mapped on BTA15, is a processed pseudogene with a similarity of 88.1% to bDDX3Y and 93.7% to bDDX3X mRNA, suggesting that PL10 is a retroposon of DDX3X. RT-PCR analyses showed that bDDX3Y-L, -S, bDDX3X and bPL10 were all widely expressed with predominant expression in testis and brain. Testicular section in situ hybridization revealed that sense and anti-sense RNAs of bDDX3Y-L, -S, and bDDX3X were expressed in interstitial cells. These results together with the finding that the pseudogene bPL10 is transcriptionally active in this study provide a basis for further investigating the DDX3 gene function in spermatogenesis, male fertility and gene evolution in mammals.

Copy number variations of the extensively amplified Y-linked genes, HSFY and ZNF280BY, in cattle and their association with male reproductive traits in Holstein bulls.

BackgroundRecent transcriptomic analysis of the bovine Y chromosome revealed at least six multi-copy protein coding gene families, including TSPY, HSFY and ZNF280BY, on the male-specific region (MSY). Previous studies indicated that the copy number variations (CNVs) of the human and bovine TSPY were associated with male fertility in men and cattle. However, the relationship between CNVs of the bovine Y-linked HSFY and ZNF280BY gene families and bull fertility has not been investigated.ResultsWe investigated the copy number (CN) of the bovine HSFY and ZNF280BY in a total of 460 bulls from 15 breeds using a quantitative PCR approach. We observed CNVs for both gene families within and between cattle breeds. The median copy number (MCN) of HSFY among all bulls was 197, ranging from 21 to 308. The MCN of ZNF280BY was 236, varying from 28 to 380. Furthermore, bulls in the Bos taurus (BTA) lineage had a significantly higher MCN (202) of HSFY than bulls in the Bos indicus (BIN) lineage (178), while taurine bulls had a significantly lower MCN (231) of ZNF280BY than indicine bulls (284). In addition, the CN of ZNF280BY was positively correlated to that of HSFY on the BTAY. Association analysis revealed that the CNVs of both HSFY and ZNF280BY were correlated negatively with testis size, while positively with sire conception rate.ConclusionThe bovine HSFY and ZNF280BY gene families have extensively expanded on the Y chromosome during evolution. The CN of both gene families varies significantly among individuals and cattle breeds. These variations were associated with testis size and bull fertility in Holstein, suggesting that the CNVs of HSFY and ZNF280BY may serve as valuable makers for male fertility selection in cattle.

The molecular evolution of PL10 homologs

BackgroundPL10 homologs exist in a wide range of eukaryotes from yeast, plants to animals. They share a DEAD motif and belong to the DEAD-box polypeptide 3 (DDX3) subfamily with a major role in RNA metabolism. The lineage-specific expression patterns and various genomic structures and locations of PL10 homologs indicate these homologs have an interesting evolutionary history.ResultsPhylogenetic analyses revealed that, in addition to the sex chromosome-linked PL10 homologs, DDX3X and DDX3Y, a single autosomal PL10 putative homologous sequence is present in each genome of the studied non-rodent eutheria. These autosomal homologous sequences originated from the retroposition of DDX3X but were pseudogenized during the evolution. In rodents, besides Ddx3x and Ddx3y, we found not only Pl10 but another autosomal homologous region, both of which also originated from the Ddx3x retroposition. These retropositions occurred after the divergence of eutheria and opossum. In contrast, an additional X putative homologous sequence was detected in primates and originated from the transposition of DDX3Y. The evolution of PL10 homologs was under positive selection and the elevated Ka/Ks ratios were observed in the eutherian lineages for DDX3Y but not PL10 and DDX3X, suggesting relaxed selective constraints on DDX3Y. Contrary to the highly conserved domains, several sites with relaxed selective constraints flanking the domains in the mammalian PL10 homologs may play roles in enhancing the gene function in a lineage-specific manner.ConclusionThe eutherian DDX3X/DDX3Y in the X/Y-added region originated from the translocation of the ancient PL10 ortholog on the ancestral autosome, whereas the eutherian PL10 was retroposed from DDX3X. In addition to the functional PL10/DDX3X/DDX3Y, conserved homologous regions on the autosomes and X chromosome are present. The autosomal homologs were also derived from DDX3X, whereas the additional X-homologs were derived from DDX3Y. These homologs were apparently pseudogenized but may still be active transcriptionally. The evolution of PL10 homologs was positively selected.

Phylogeography and Domestication of Chinese Swamp Buffalo

To further probe into whether swamp buffaloes were domesticated once or multiple times in China, this survey examined the mitochondrial DNA (mtDNA) Control Region (D-loop) diversity of 471 individuals representing 22 populations of 455 Chinese swamp buffaloes and 16 river buffaloes. Phylogenetic analysis revealed that Chinese swamp buffaloes could be divided into two distinct lineages, A and B, which were defined previously. Of the two lineages, lineage A was predominant across all populations. For predominant lineage A, Southwestern buffalo populations possess the highest genetic diversity among the three hypothesized domestication centers (Southeastern, Central, and Southwestern China), suggesting Southwestern China as the most likely location for the domestication of lineage A. However, a complex pattern of diversity is detected for the lineage B, preventing the unambiguous pinpointing of the exact place of domestication center and suggesting the presence of a long-term, strong gene flow among swamp buffalo populations caused by extensive migrations of buffaloes and frequent human movements along the Yangtze River throughout history. Our current study suggests that Southwestern China is the most likely domestication center for lineage A, and may have been a primary center of swamp buffalo domestication. More archaeological and genetic evidence is needed to show the process of domestication.

Differential Expression of PRAMEL1, a Cancer/Testis Antigen, during Spermatogenesis in the Mouse

PRAME belongs to a group of cancer/testis antigens (CTAs) that are characterized by their restricted expression in normal gametogenic tissues and a variety of tumors. The PRAME family is one of the most amplified gene families in the mouse and other mammalian genomes. Members of the PRAME gene family encode leucine-rich repeat (LRR) proteins functioning as transcription regulators in cancer cells. However, the role of PRAME in normal gonads is unknown. The objective of this study is to characterize the temporal and spatial expression of the mouse Pramel1 gene, and to determine the cellular localization of the PRAMEL1 protein during the mouse spermatogenesis. Our results indicated that the mouse Pramel1 was expressed in testis only. The mRNA and protein expression level was low in the newborn testes, and gradually increased from 1- to 3-week-old testes, and then remained constant after three weeks of age. Immunofluorescent staining on testis sections with the mouse PRAMEL1 antibody revealed that PRAMEL1 was localized in the cytoplasm of spermatocytes and the acrosomal region of round, elongating and elongated spermatids. Further analyses on the testis squash preparation and spermatozoa at a subcellular level indicated that the protein localization patterns of PRAMEL1 were coordinated with morphological alterations during acrosome formation in spermatids, and were significantly different in connecting piece, middle piece and principal piece of the flagellum between testicular and epididymal spermatozoa. Collectively, our results suggest that PRAMEL1 may play a role in acrosome biogenesis and sperm motility.

Regional Selection Acting on the OFD1 Gene Family

The OFD1 (oral-facial-digital, type 1) gene is implicated in several developmental disorders in humans. The X-linked OFD1 (OFD1X) is conserved in Eutheria. Knowledge about the Y-linked paralog (OFD1Y) is limited. In this study, we identified an OFD1Y on the bovine Y chromosome, which is expressed differentially from the bovine OFD1X. Phylogenetic analysis indicated that: a) the eutherian OFD1X and OFD1Y were derived from the pair of ancestral autosomes during sex chromosome evolution; b) the autosomal OFD1 pseudogenes, present in Catarrhini and Murinae, were derived from retropositions of OFD1X after the divergence of primates and rodents; and c) the presence of OFD1Y in the ampliconic region of the primate Y chromosome is an indication that the expansion of the ampliconic region may initiate from the X-degenerated sequence. In addition, we found that different regions of OFD1/OFD1X/OFD1Y are under differential selection pressures. The C-terminal half of OFD1 is under relaxed selection with an elevated Ka/Ks ratio and clustered positively selected sites, whereas the N-terminal half is under stronger constraints. This study provides some insights into why the OFD1X gene causes OFD1 (male-lethal X-linked dominant) and SGBS2 & JSRDs (X-linked recessive) syndromes in humans, and reveals the origin and evolution of the OFD1 family, which will facilitate further clinical investigation of the OFD1-related syndromes.

A high-resolution comparative map of porcine chromosome 4 (SSC4)

We used the IMNpRH2(12,000-rad) RH and IMpRH(7,000-rad) panels to integrate 2019 transcriptome (RNA-seq)-generated contigs with markers from the porcine genetic and radiation hybrid (RH) maps and bacterial artificial chromosome finger-printed contigs, into 1) parallel framework maps (LOD ≥ 10) on both panels for swine chromosome (SSC) 4, and 2) a high-resolution comparative map of SSC4, thus and human chromosomes (HSA) 1 and 8. A total of 573 loci were anchored and ordered on SSC4 closing gaps identified in the porcine sequence assembly Sscrofa9. Alignment of the SSC4 RH with the genetic map identified five microsatellites incorrectly mapped around the centromeric region in the genetic map. Further alignment of the RH and comparative maps with the genome sequence identified four additional regions of discrepancy that are also suggestive of errors in assembly, three of which were resolved through conserved synteny with blocks on HSA1 and HSA8.