Testis transcriptome profiling identified lncRNAs involved in spermatogenic arrest of cattleyak

Transcriptome analysis reveals dysregulated gene expression networks in Sertoli cells of cattle-yak hybrids

Cattle-yak, the hybrid offspring of yak and taurine cattle, exhibits male sterility with normal female fertility. Spermatogenesis is arrested in adult cattle-yak, and apoptosis is elevated in spermatogenic cells. Currently, the mechanisms underlying these defects remain elusive. Sertoli cells are the only somatic cells that directly interact with spermatogenic cells in the seminiferous tubules and play essential roles in spermatogenesis. The present study was designed to investigate gene expression signatures and potential roles of Sertoli cells in hybrid sterility in cattle-yak. Immunohistochemical analysis showed that the 5 mC and 5hmC signals in Sertoli cells of cattle-yaks were significantly different from those of age-matched yaks (P < 0.05). Transcriptome profiling of isolated Sertoli cells identified 402 differentially expressed genes (DEGs) between cattle-yaks and yaks. Notably, niche factor glial cell derived neurotrophic factor (GDNF) was upregulated, and genes involved in retinoic acid (RA) biogenesis were changed in Sertoli cells of cattle-yak, suggesting possible impairments of spermatogonial fate decisions. Further studies showed that the numbers of proliferative gonocytes and undifferentiated spermatogonia in cattle-yak were significantly higher than those in yak (P < 0.01). Exogenous GDNF significantly promoted the proliferation of UCHL1-positive spermatogonia in yaks. Therefore, we concluded that altered GDNF expression and RA signaling impacted the fate decisions of undifferentiated spermatogonia in cattle-yak. Together, these findings highlight the role of Sertoli cells and their derived factors in hybrid sterility.

Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

A model based on immune-related lncRNA pairs and its potential prognostic value in immunotherapy for melanoma

A model based on long non-coding RNA (lncRNA) pairs independent of expression quantification was constructed to evaluate prognosis melanoma and response to immunotherapy in melanoma. RNA sequencing data and clinical information were retrieved and downloaded from The Cancer Genome Atlas and the Genotype-Tissue Expression databases. We identified differentially expressed immune-related lncRNAs (DEirlncRNAs), matched them, and used least absolute shrinkage and selection operator and Cox regression to construct predictive models. The optimal cutoff value of the model was determined using a receiver operating characteristic curve and used to categorize melanoma cases into high-risk and low-risk groups. The predictive efficacy of the model with respect to prognosis was compared with that of clinical data and ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data). Then, we analyzed the correlations of risk score with clinical characteristics, immune cell invasion, anti-tumor, and tumor-promoting activities. Differences in survival, degree of immune cell infiltration, and intensity of anti-tumor and tumor-promoting activities were also evaluated in the high- and low-risk groups. A model based on 21 DEirlncRNA pairs was established. Compared with ESTIMATE score and clinical data, this model could better predict outcomes of melanoma patients. Follow-up analysis of the model's effectiveness showed that patients in the high-risk group had poorer prognosis and were less likely to benefit from immunotherapy compared with those in the low-risk group. Moreover, there were differences in tumor-infiltrating immune cells between the high-risk and low-risk groups. By pairing the DEirlncRNA, we constructed a model to evaluate the prognosis of cutaneous melanoma independent of a specific level of lncRNA expression.

Whole transcriptome sequencing reveals core genes related to spermatogenesis in Bactrian camels

Bactrian camels survive and reproduce better in extreme climatic conditions than other domestic animals can. However, the reproductive efficiency of camels under their natural pastoral conditions is low. Several factors affect mammalian reproductive performance, including testicular development, semen quality, libido, and mating ability. Testis is a main reproductive organ of the male and is responsible for producing spermatozoa and hormones. However, our understanding of the expression patterns of the genes in camel testis is minimal. Thus, we performed total RNA-sequencing to investigate the gene expression pattern. As a result, 1,538 differential expressed mRNAs (DEmRNAs), 702 differential expressed long non-coding RNAs (DElncRNAs), and 61 differential expressed microRNAs (DEmiRNAs) were identified between pubertal and adult Bactrian camel testes. Then the genomic features, length distribution, and other characteristics of the lncRNAs and mRNAs in the Bactrian camel testis were investigated. Target genes of the DEmiRNAs and DEmRNAs were further subjected to gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Genes, such as AMHR2, FGF1, ACTL7A, GATA4, WNT4, ID2, LAMA1, IGF1, INHBB, and TLR2, were mainly involved in the TGF-β, PI3K-AKT, Wnt, GnRH, and Hippo signaling pathways which relate to spermatogenesis. Some of the DEmiRNAs were predicted to be associated with numerous DElncRNAs and DEmRNAs through competing endogenous RNA (ceRNA) regulatory network. At last, the candidate genes were validated by RT-qPCR, dual fluorescent reporter gene, and a fluorescence in situ hybridization (FISH) assay. This research provides high-throughput RNA sequencing data of the testes of Bactrian camels across different developmental stages. It lays the foundation for further investigations on lncRNAs, miRNAs, and mRNAs that involved in Bactrian camel spermatogenesis.

Identification of genomic regions associated with total and progressive sperm motility in Italian Holstein bulls

Sperm motility is directly related to the ability of sperm to move through the female reproductive tract to reach the ovum. Sperm motility is a complex trait that is influenced by environmental and genetic factors and is associated with male fertility, oocyte penetration rate, and reproductive success of cattle. In this study we carried out a GWAS in Italian Holstein bulls to identify candidate regions and genes associated with variations in progressive and total motility (PM and TM, respectively). After quality control, the final data set consisted of 5,960 records from 949 bulls having semen collected in 10 artificial insemination stations and genotyped at 412,737 SNPs (call rate >95%; minor allele frequency >5%). (Co)variance components were estimated using single trait mixed models, and associations between SNPs and phenotypes were assessed using a genomic BLUP approach. Ten windows that explained the greatest percentage of genetic variance were located on Bos taurus autosomes 1, 2, 4, 6, 7, 23, and 26 for TM and Bos taurus autosomes 1, 2, 4, 6, 8, 16, 23, and 26 for PM. A total of 150 genes for TM and 72 genes for PM were identified within these genomic regions. Gene Ontology enrichment analyses identified significant Gene Ontology terms involved with energy homeostasis, membrane functions, sperm-egg interactions, protection against oxidative stress, olfactory receptors, and immune system. There was significant enrichment of quantitative trait loci for fertility, calving ease, immune response, feed intake, and carcass weight within the candidate windows. These results contribute to understanding the architecture of the genetic control of sperm motility and may aid in the development of strategies to identify subfertile bulls and improve reproductive success.

Analysis of Chromatin Openness in Testicle Tissue of Yak and Cattle-Yak

Cattle-yak, a crossbreed of yak and cattle, which can exhibit obvious heterosis and can adapt to the harsh environmental conditions of the Qinghai Tibet Plateau (QTP). However, F1 cattle-yak were found to be sterile because they were unable to produce sperm, which adversely restricted the fixation of heterosis. Many prior attempts have been made to decipher the mechanism underlying the spermatogenesis stagnation of cattle-yak. However, the open chromatin region (OCR) map of yak and cattle-yak testes has not been generated yet. Here, we have analyzed the OCRs landscape of testicular tissues of cattle-yak and yaks by performing ATAC-seq technology. The OCRs of cattle-yak and yak testes displayed similar genome distribution and showed priority in intergenic regions, introns and promoters. The pathway enrichment analysis indicated that the differential OCRs-related genes were involved in spermatogenesis, involving the cell cycle, as well as Hippo, mTOR, MAPK, Notch, and Wnt signaling pathways. The integration of ATAC-seq and mRNA-seq indicated that the majority of the gene expression levels were positively correlated with chromatin openness. At the same time, we have identified a number of transcription factors (TFs) related to spermatogenesis and the differential expression of these TFs may contribute to the spermatogenesis stagnation of the cattle-yak. Overall, the findings of this study provide valuable information for advancing the research related to yak crossbreeding improvement and sperm production stagnation of cattle-yak.

First detection and complete genome analysis of porcine circovirus-like virus P1 and porcine circovirus-2 in yak in China

Porcine circovirus-like virus P1, like porcine circovirus type 2 (PCV2), is a potential pathogen of post-weaning multisystemic wasting syndrome in swine. Yaks are a valuable species and an iconic symbol of the Tibet Plateau which is the highest and largest plateau in the world. In this study, a total of 105 yak diarrheal samples, collected from 13 farms in Linzhi in the Tibet Plateau from January 2019 to December 2021, that were screened for P1 and PCV2 by polymerase chain reaction, 10.48% (n = 11) were positive for P1, 4.76% (n = 5) for PCV2, and 5.71% (n = 6) were positive for coinfection of P1 and PCV2. In addition, the whole genomes of eight P1 strains and eight PCV2 strains were sequenced. Alignment of deduced amino acid sequences of P1 ORF1 and PCV2 ORF2 gene revealed that ON012566 had one unique amino acid mutation at residues 137 (T to P). This mutation has important implication for the study of virus virulence, tissue tropism, and immune response. Phylogenetic analysis shows that the yak-origin P1 strains in this study with cattle-origin P1 reference strains were grouped into one cluster. The yak-origin PCV2 (ON012566) and a buffalo-origin PCV2 (KM116514) reference strain clustered in the same branch in the PCV2b regions. Meanwhile, the remaining PCV2 strains and buffalo-origin PCV2 reference strain (ON012565) clustered in the PCV2d regions. To summarize, to our knowledge, this is the first report on the molecular prevalence and genome characteristics of P1 and PCV2 in yaks in the world and will contribute to further study of the molecular epidemiology, source, and evolution of P1 and PCV2 strains.

Proteome Informatics in Tibetan Sheep (Ovis aries) Testes Suggest the Crucial Proteins Related to Development and Functionality

Testis has an indispensable function in male reproduction of domestic animals. Tibetan sheep (Ovis aries) is a locally adapted breed of sheep raised in the Qinghai-Tibet Plateau, with outsized roles in providing the livelihood for millions of residents. Nevertheless, less is known on how protein expression and their functional roles in developmental testes of such breed limit their use in breeding efforts. In this study, we obtained comprehensive protein profiles from testes of Tibetan sheep at three developmental stages (including pre-puberty, post-puberty, and adulthood) using data-independent acquisition-based proteomic strategy to quantitatively identify the differentially abundant proteins (DAPs) associated with testicular development and function and to unravel the molecular basis of spermatogenesis. A total of 6,221 proteins were differentially expressed in an age-dependent manner. The reliability of the gene expression abundance was corroborated by quantitative PCR and targeted parallel reaction monitoring. These DAPs were significantly enriched to biological processes concerning spermatid development and sperm deformation, mitosis, glycolytic process, cell-cell/extracellular matrix (ECM) junctions, cell proliferation, apoptosis, and migration and to the pathways including, developmental process and sexual reproduction-related (such as VEGF, estrogen, insulin, GnRH, Hippo, PI3K-Akt, mTOR, MAPK, and AMPK), and testicular cell events-related pathways (such as tight/gap/adherens junctions, ECM-receptor interaction, regulation of actin cytoskeleton, glycolysis, cell cycle, and meiosis). Based on these bioinformatics analysis, we constructed four protein–protein interaction network, among which the proteins are involved in mitosis, meiosis, spermiogenesis, and testicular microenvironment, respectively. Altogether, these bioinformatics-based sequencing results suggest that many protein-coding genes were expressed in a development-dependent manner in Tibetan sheep testes to contribute to the testicular cell development and their surrounding microenvironment remodeling at various stages of spermatogenesis. These findings have important implications for further understanding of the mechanisms underlying spermatogenesis in sheep and even other plateau-adapted animals.