Comparative transcriptome in large-scale human and cattle populations

A pilot study of transcriptomic preimplantation genetic testing (PGT-T): towards a new step in embryo selection?

STUDY QUESTION

Is it possible to predict an euploid chromosomal constitution and identify a transcriptomic profile compatible with extended embryonic development from RNA sequencing (RNA-Seq) data?

SUMMARY ANSWER

It has been possible to obtain a karyotype comparable to preimplantation genetic testing for aneuploidy (PGT-A), in addition to a transcriptomic signature of embryos which might be suggestive of improved implantation capacity.

WHAT IS KNOWN ALREADY

Conventional assessment of embryo competence, based on morphology and morphokinetic, lacks knowledge of molecular aspects and faces controversy in predicting ploidy status. Understanding the embryonic transcriptome is crucial, as gene expression influences development and implantation. PGT has improved pregnancy rates, but problems persist when high-quality euploid embryos do not reach term. In fact, only around 50-60% implant, of which 10% result in miscarriage. Comprehensive approaches, including RNA-Seq, offer the potential to discover molecular markers of reproductive competence, and could theoretically be combined with extended-embryo culture platforms up to Day 14 that can be utilized as a proxy to study embryo development at post-implantation stages.

STUDY DESIGN, SIZE, DURATION

This prospective pilot cohort study was conducted from March 2023 to August 2023. A total of 30 vitrified human blastocysts with previous PGT-A diagnosis on Day 5 (D5) or Day 6 (D6) of development were analysed: n = 15 euploid and n = 15 aneuploid. Finally, 21 embryo samples were included in the study; the rest (n = 9) were excluded due to poor quality pre-sequencing data (n = 7) or highly discordant data (n = 2).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Following warming and re-expansion, embryos underwent a second trophectoderm (TE) biopsy. The embryos were then cultured until day 11 to assess their development. Biopsy analysis by RNA-Seq, studied the differential expressed genes (DEG) to compare embryos which did not or did attach to the plate: unattached embryos (n = 12) versus attached embryos (n = 9). Thus, we also obtained a specific transcriptomic signature of embryos with a "theoretical" capacity for sustained implantation, based on plate attachment on day 11.

MAIN RESULTS AND THE ROLE OF CHANCE

The digital karyotype obtained by RNA-Seq showed good concordance with the earlier PGT-A data, with a sensitivity of 0.81, a specificity of 0.83, a Cohen's Kappa of 0.66, and an area under the ROC of 0.9. At the gene level, 76 statistically significant DEGs were found in the comparison unattached versus attached embryos (Padj < 0.05; FC > 1). To address the functional implications of these differences, significantly deregulated pathways according to GO and KEGG categories were identified. The mural trophectoderm (TE) of the unattached blastocysts showed 63 significantly deregulated terms, displaying upregulation in autophagy, apoptosis, protein kinase and ubiquitin-like protein ligase activity, and downregulation of ribosome, spliceosome, kinetochore, segregation, and chromosome condensation processes. The overall transcriptomic signature specific to embryos still attached to the plate on day 11 (with a theoretically higher implantation capacity) consists of 501 genes, including: EMP2, AURKB, FOLR1, NOTCH3, LRP2, FZD5, MDH1, APOD, GPX8, COLEC12, HSPA1A, CMTM7, BEX3, which are related to implantation and embryonic development (raw P-value < 0.05; shrunk LFC > 1.1). These findings indicate that it might be possible to identify euploid embryos with a greater capacity for implantation and development, after excluding those embryos that present chromosomal alterations.

LIMITATIONS, REASONS FOR CAUTION

This study included a small sample size, remarkable variability between samples, and low success rate of RNA amplification. Also, structural chromosomal abnormalities were not included, and it was not possible to diagnose mosaic embryos. TE biopsy does not assure the chromosomal status of the whole embryo. The maximum day for in vitro development was Day 11, and attachment to the plate on this day does not provide a clear indication of implantation capacity and viability, which was not tested in this study.

WIDER IMPLICATIONS OF THE FINDINGS

The short-term goals following on from this pilot study is to expand the sample size with embryos of more complex abnormalities, and to perform a prospective in vitro preclinical validation. In a more distant future and with optimal results, this technique could have clinical application, thus increasing clinical outcomes by assessing both chromosomal content and transcriptomic profiling.

STUDY FUNDING/COMPETING INTEREST(S)

The Institut Valencià de Competitivitat Empresarial (IVACE) (IMIDCA/2022/39) and Generalitat Valenciana (CIACIF/2021/11) supported the present study. A.C. is an employee of JUNO Genetics. He has received honoraria for an IBSA lecture and a Merck lecture. He is also a minor shareholder of IVIRMA Global. The other authors have no conflicts of interest to declare.

TRIAL REGISTRATION NUMBER

N/A.

Do Delta-9-tetrahydrocannabinol and Cannabidiol have opposed effects on male fertility?

Cannabis sativa is a complex plant, renowned for its diverse array of bioactive compounds, the most prominent of which are delta-9-tetrahydrocannabinol (THC) and cannabidiol (CBD). These compounds exhibit markedly opposing pharmacological effects, with THC being primarily psychoactive and CBD known for its non-psychoactive properties. In recent years, there has been growing interest in the potential health implications of these compounds, particularly concerning male reproductive health. Accumulating evidence over the past decade has alluded to the potential negative effects of THC, including its association with reduced sperm quality, altered hormone levels, changes in genetic and epigenetic profiles, and potential impacts on fertility. Conversely, emerging studies suggest that CBD may exert protective and beneficial effects on male reproductive health, possibly through its anti-inflammatory and antioxidant properties. This review aims to provide a comprehensive analysis of the current scientific literature, delineating the mechanisms by which THC and CBD influence male reproductive health, highlighting the disparities in their effects, and discussing the clinical and therapeutic implications of these findings.

ZFAT (isoform-specific) and its antisense RNA 1 (ZFAT-AS1) are two allele-specific monoallelically expressed genes in cattle

In mammals, imprinted genes are characterised by a monoallelic expression, which is based on parental origin and is essential for both foetal and placental development. The ZFAT gene encodes a transcriptional factor, and its non-coding antisense RNA, ZFAT-AS1, overlaps with the ZFAT locus. Both ZFAT and ZFAT-AS1 are maternally imprinted in human placentas. In bovines, the imprinting status of the ZFAT and ZFAT-AS1 genes has yet to be reported. In this study, we analysed the allelic expression of three transcript variants (X1-X3) of the bovine ZFAT and ZFAT-AS1 genes in somatic tissues and placentas using a single nucleotide polymorphism-based method. The results showed that bovine ZFAT exhibited isoform-specific paternal expression. The ZFAT X2 variant exhibited monoallelic expression in the bovine placentas and biallelic expression in the six bovine somatic tissues (heart, liver, spleen, lung, kidney and brain). However, the ZFAT X1 and X3 variants were biallelically expressed in both bovine tissues and placentas. A 311 bp bovine ZFAT-AS1 complementary DNA (cDNA) sequence was obtained by aligning the human ZFAT-AS1 cDNA sequence with the bovine genome and conducting reverse transcription polymerase chain reaction amplification. Bovine ZFAT-AS1 have monoallelic expression in bovine placentas and somatic tissues. In addition, the DNA methylation of two regions was characterised, including the partial promoter, and exon 1 and intron 1 regions of ZFAT, and there were no differentially methylated regions.

Adipose gene expression profiles in Northern Finncattle, Mirandesa cattle, Yakutian cattle and commercial Holstein cattle

The drastic change in global climate has led to in-depth studies of the geneticresources of native cattle adapted to challenging environments. Native cattle breeds may harbor unique genetic mechanisms that have enabled them adapt to their given environmental conditions. Adipose tissues are key factors in the regulation of metabolism and energy balance and are crucial for the molecular switches needed to adapt to rapid environmental and nutritional changes. The transcriptome landscape of four adipose tissues was used in this study to investigate the differential gene expression profiles in three local breeds, Yakutian cattle (Sakha Republic), Northern Finncattle (Finland), Mirandesa cattle (Portugal) and commercial Holstein cattle. A total of 26 animals (12 cows, 14 bulls) yielded 81 samples of perirenal adipose tissue (n = 26), metacarpal adipose tissue (n = 26), tailhead adipose tissue (n = 26) and prescapular adipose tissue (n = 3). More than 17,000 genes were expressed in our dataset. Principal component analysis of the normalized expression profiles revealed a differential expression profile of the metacarpal adipose tissue. We found that the genes upregulated in the metacarpal adipose tissue of Yakutian cattle, such as NR4A3, TEKT3, and FGGY, were associated with energy metabolism and response to cold temperatures. In Mirandesa cattle, the upregulated genes in perirenal adipose tissue were related to immune response and inflammation (AVPR2, CCN1, and IL6), while in Northern Finncattle, the upregulated genes appeared to be involved in various physiological processes, including energy metabolism (IGFBP2). According to the sex-based comparisons, the most interesting result was the upregulation of the TPRG1 gene in three tissues of Yakutian cattle females, suggesting that adaptation is related to feed efficiency. The highest number of differentially expressed genes was found between Yakutian cattle and Holstein, several of which were associated with immunity in Yakutian cattle, indicating potential differences in disease resistance and immunity between the two breeds. This study highlights the vast difference in gene expression profiles in adipose tissues among breeds from different climatic environments, most likely highlighting selective pressure and the potential significance of the uniquely important regulatory functions of metacarpal adipose tissue.

Comparative analysis of testicular fusion in Spodoptera litura (cutworm) and Bombyx mori (silkworm): Histological and transcriptomic insights

Spodoptera litura commonly known as the cutworm, is among the most destructive lepidopteran pests affecting over 120 plants species. The powerful destructive nature of this lepidopteran is attributable to its high reproductive capacity. The testicular fusion that occurs during metamorphosis from larvae to pupa in S.litura positively influences the reproductive success of the offspring. In contrast, Bombyx mori, the silkworm, retains separate testes throughout its life and does not undergo this fusion process. Microscopic examination reveals that during testicular fusion in S.litura, the peritoneal sheath becomes thinner and more translucent, whereas in B.mori, the analogous region thickens. The outer basement membrane in S.litura exhibits fractures, discontinuity, and uneven thickness accompanied by a significant presence of cellular secretions, large cell size, increased vesicles, liquid droplets, and a proliferation of rough endoplasmic reticulum and mitochondria. In contrast, the testicular peritoneal sheath of B.mori at comparable developmental stage exhibits minimal change. Comparative transcriptomic analysis of the testicular peritoneal sheath reveals a substantial difference in gene expression between the two species. The disparity in differential expressed genes (DEGs) is linked to an enrichment of numerous transcription factors, intracellular signaling pathways involving Ca2+ and GTPase, as well as intracellular protein transport and signaling pathways. Meanwhile, structural proteins including actin, chitin-binding proteins, membrane protein fractions, cell adhesion, extracellular matrix proteins are predominantly identified. Moreover, the study highlights the enrichment of endopeptidases, serine proteases, proteolytic enzymes and matrix metalloproteins, which may play a role in the degradation of the outer membrane. Five transcription factors-Slforkhead, Slproline, Slcyclic, Slsilk, and SlD-ETS were identified, and their expression pattern were confirmed by qRT-PCR. they are candidates for participating in the regulation of testicular fusion. Our findings underscore significant morphological and trancriptomic variation in the testicular peritoneal sheath of S.litura compared to the silkworm, with substantial changes at the transcriptomic level coinciding with testicular fusion. The research provides valuable clues for understanding the complex mechanisms underlying this unique phenomenon in insects.

Meiotic maturation failure in primary ovarian insufficiency: insights from a bovine model

PURPOSE

Oocytes from women presenting primary ovarian insufficiency (POI) generate viable embryos at a lower rate than non-POI women, but the mechanisms responsible for the lower oocyte quality remain elusive. Due to the scarcity of human oocytes for research, animal models provide a promising way forward. We aimed at investigating the molecular events characterizing final maturation in POI oocytes in a well-defined POI-like bovine model.

METHODS

Single-cell RNA-sequencing of bovine control and POI-like, GV, and MII oocytes (n = 5 per group) was performed. DEseq2 was used to identify differentially expressed genes. Further, a Gene set enrichment analysis and a transcriptomic meta-analysis between bovine and human oocytes were performed.

RESULTS

In control cows, we found 2223 differentially expressed genes between the GV and MII stages. Specifically, the affected genes were related to RNA processing and transport, protein synthesis, organelle remodeling and reorganization, and metabolism. The meta-analysis with a set of young human oocytes at different maturation stages revealed 315 conserved genes through the GV-MII transition in cows and humans, mostly related to meiotic progression and cell cycle. Gene expression analysis between GV and MII of POI-like oocytes showed no differences in terms of differentially expressed genes, pointing towards a substantial failure to properly remodel the transcriptome in the POI model, and with the clustering analysis indicating that the cow's genetic background had a higher impact than the oocyte's maturation stage.

CONCLUSION

Overall, we have identified and characterized a valuable animal model of POI, paving the way to identifying new molecular mechanisms involved in POI.

Improving the accuracy of genomic prediction in dairy cattle using the biologically annotated neural networks framework

BACKGROUND

Biologically annotated neural networks (BANNs) are feedforward Bayesian neural network models that utilize partially connected architectures based on SNP-set annotations. As an interpretable neural network, BANNs model SNP and SNP-set effects in their input and hidden layers, respectively. Furthermore, the weights and connections of the network are regarded as random variables with prior distributions reflecting the manifestation of genetic effects at various genomic scales. However, its application in genomic prediction has yet to be explored.

RESULTS

This study extended the BANNs framework to the area of genomic selection and explored the optimal SNP-set partitioning strategies by using dairy cattle datasets. The SNP-sets were partitioned based on two strategies-gene annotations and 100 kb windows, denoted as BANN_gene and BANN_100kb, respectively. The BANNs model was compared with GBLUP, random forest (RF), BayesB and BayesCπ through five replicates of five-fold cross-validation using genotypic and phenotypic data on milk production traits, type traits, and one health trait of 6,558, 6,210 and 5,962 Chinese Holsteins, respectively. Results showed that the BANNs framework achieves higher genomic prediction accuracy compared to GBLUP, RF and Bayesian methods. Specifically, the BANN_100kb demonstrated superior accuracy and the BANN_gene exhibited generally suboptimal accuracy compared to GBLUP, RF, BayesB and BayesCπ across all traits. The average accuracy improvements of BANN_100kb over GBLUP, RF, BayesB and BayesCπ were 4.86%, 3.95%, 3.84% and 1.92%, and the accuracy of BANN_gene was improved by 3.75%, 2.86%, 2.73% and 0.85% compared to GBLUP, RF, BayesB and BayesCπ, respectively across all seven traits. Meanwhile, both BANN_100kb and BANN_gene yielded lower overall mean square error values than GBLUP, RF and Bayesian methods.

CONCLUSION

Our findings demonstrated that the BANNs framework performed better than traditional genomic prediction methods in our tested scenarios, and might serve as a promising alternative approach for genomic prediction in dairy cattle.

Integrated analysis of genomics and transcriptomics revealed the genetic basis for goaty flavor formation in goat milk

Goat milk exhibits a robust and distinctive "goaty" flavor. However, the underlying genetic basis of goaty flavor remains elusive and requires further elucidation at the genomic level. Through comparative genomics analysis, we identified divergent signatures of certain proteins in goat, sheep, and cow. MMUT has undergone a goat-specific mutation in the B12 binding domain. We observed the goat FASN exhibits nonsynonymous mutations in the acyltransferase domain. Structural variations in these key proteins may enhance the capacity for synthesizing goaty flavor compounds in goat. Integrated omics analysis revealed the catabolism of branched-chain amino acids contributed to the goat milk flavor. Furthermore, we uncovered a regulatory mechanism in which the transcription factor ZNF281 suppresses the expression of the ECHDC1 gene may play a pivotal role in the accumulation of flavor substances in goat milk. These findings provide insights into the genetic basis underlying the formation of goaty flavor in goat milk. STATEMENT OF SIGNIFICANCE: Branched-chain fatty acids (BCFAs) play a crucial role in generating the distinctive "goaty" flavor of goat milk. Whether there is an underlying genetic basis associated with goaty flavor is unknown. To begin deciphering mechanisms of goat milk flavor development, we collected transcriptomic data from mammary tissue of goat, sheep, cow, and buffalo at peak lactation for cross-species transcriptome analysis and downloaded nine publicly available genomes for comparative genomic analysis. Our data indicate that the catabolic pathway of branched-chain amino acids (BCAAs) is under positive selection in the goat genome, and most genes involved in this pathway exhibit significantly higher expression levels in goat mammary tissue compared to other species, which contributes to the development of flavor in goat milk. Furthermore, we have elucidated the regulatory mechanism by which the transcription factor ZNF281 suppresses ECHDC1 gene expression, thereby exerting an important influence on the accumulation of flavor compounds in goat milk. These findings provide insights into the genetic mechanisms underlying flavor formation in goat milk and suggest further research to manipulate the flavor of animal products.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

A time-resolved multi-omics atlas of transcriptional regulation in response to high-altitude hypoxia across whole-body tissues

High-altitude hypoxia acclimatization requires whole-body physiological regulation in highland immigrants, but the underlying genetic mechanism has not been clarified. Here we use sheep as an animal model for low-to-high altitude translocation. We generate multi-omics data including whole-genome sequences, time-resolved bulk RNA-Seq, ATAC-Seq and single-cell RNA-Seq from multiple tissues as well as phenotypic data from 20 bio-indicators. We characterize transcriptional changes of all genes in each tissue, and examine multi-tissue temporal dynamics and transcriptional interactions among genes. Particularly, we identify critical functional genes regulating the short response to hypoxia in each tissue (e.g., PARG in the cerebellum and HMOX1 in the colon). We further identify TAD-constrained cis-regulatory elements, which suppress the transcriptional activity of most genes under hypoxia. Phenotypic and transcriptional evidence indicate that antenatal hypoxia could improve hypoxia tolerance in offspring. Furthermore, we provide time-series expression data of candidate genes associated with human mountain sickness (e.g., BMPR2) and high-altitude adaptation (e.g., HIF1A). Our study provides valuable resources and insights for future hypoxia-related studies in mammals.

Comparative transcriptome analysis between rhesus macaques ( Macaca mulatta) and crab-eating macaques ( M. fascicularis)

Understanding gene expression variations between species is pivotal for deciphering the evolutionary diversity in phenotypes. Rhesus macaques ( Macaca mulatta, MMU) and crab-eating macaques ( M. fascicularis, MFA) serve as crucial nonhuman primate biomedical models with different phenotypes. To date, however, large-scale comparative transcriptome research between these two species has not yet been fully explored. Here, we conducted systematic comparisons utilizing newly sequenced RNA-seq data from 84 samples (41 MFA samples and 43 MMU samples) encompassing 14 common tissues. Our findings revealed a small fraction of genes (3.7%) with differential expression between the two species, as well as 36.5% of genes with tissue-specific expression in both macaques. Comparison of gene expression between macaques and humans indicated that 22.6% of orthologous genes displayed differential expression in at least two tissues. Moreover, 19.41% of genes that overlapped with macaque-specific structural variants showed differential expression between humans and macaques. Of these, the FAM220A gene exhibited elevated expression in humans compared to macaques due to lineage-specific duplication. In summary, this study presents a large-scale transcriptomic comparison between MMU and MFA and between macaques and humans. The discovery of gene expression variations not only enhances the biomedical utility of macaque models but also contributes to the wider field of primate genomics.

Analysis of nonsynonymous SNPs in candidate genes that influence bovine temperament and evaluation of their effect in Brahman cattle

BACKGROUND

Temperament is an important production trait in cattle and multiple strategies had been developed to generate molecular markers to assist animal selection. As nonsynonymous single nucleotide polymorphisms are markers with the potential to affect gene functions, they could be useful to predict phenotypic effects. Genetic selection of less stress-responsive, temperamental animals is desirable from an economic and welfare point of view.

METHODS AND RESULTS

Two nonsynonymous single nucleotide polymorphisms identified in HTR1B and SLC18A2 candidate genes for temperament were analyzed in silico to determine their effects on protein structure. Those nsSNPs allowing changes in proteins were selected for a temperament association analysis in a Brahman population. Transversion effects on protein structure were evaluated in silico for each amino acid change model, revealing structural changes in the proteins of the HTR1B and SLC18A2 genes. The selected nsSNPs were genotyped in a Brahman population (n = 138), and their genotypic effects on three temperament traits were analyzed: exit velocity, pen score, and temperament score. Only the SNP rs209984404-HTR1B (C/A) showed a significant association (P = 0.0144) with pen score. The heterozygous genotype showed a pen score value 1.17 points lower than that of the homozygous CC genotype.

CONCLUSION

The results showed that in silico analysis could direct the selection of nsSNPs with the potential to change the protein. Non-synonymous single nucleotide polymorphisms causing structural changes and reduced protein stability were identified. Only rs209984404-HTR1B shows that the allele affecting protein stability was associated with the genotype linked to docility in cattle.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Learning functional conservation between human and pig to decipher evolutionary mechanisms underlying gene expression and complex traits

Assessment of genomic conservation between humans and pigs at the functional level can improve the potential of pigs as a human biomedical model. To address this, we developed a deep learning-based approach to learn the genomic conservation at the functional level (DeepGCF) between species by integrating 386 and 374 functional profiles from humans and pigs, respectively. DeepGCF demonstrated better prediction performance compared with the previous method. In addition, the resulting DeepGCF score captures the functional conservation between humans and pigs by examining chromatin states, sequence ontologies, and regulatory variants. We identified a core set of genomic regions as functionally conserved that plays key roles in gene regulation and is enriched for the heritability of complex traits and diseases in humans. Our results highlight the importance of cross-species functional comparison in illustrating the genetic and evolutionary basis of complex phenotypes.

Harnessing male germline epigenomics for the genetic improvement in cattle

Sperm is essential for successful artificial insemination in dairy cattle, and its quality can be influenced by both epigenetic modification and epigenetic inheritance. The bovine germline differentiation is characterized by epigenetic reprogramming, while intergenerational and transgenerational epigenetic inheritance can influence the offspring's development through the transmission of epigenetic features to the offspring via the germline. Therefore, the selection of bulls with superior sperm quality for the production and fertility traits requires a better understanding of the epigenetic mechanism and more accurate identifications of epigenetic biomarkers. We have comprehensively reviewed the current progress in the studies of bovine sperm epigenome in terms of both resources and biological discovery in order to provide perspectives on how to harness this valuable information for genetic improvement in the cattle breeding industry.

A Multi-Tissue Gene Expression Atlas of Water Buffalo (Bubalus bubalis) Reveals Transcriptome Conservation between Buffalo and Cattle

We generated 73 transcriptomic data of water buffalo, which were integrated with publicly available data in this species, yielding a large dataset of 355 samples representing 20 major tissue categories. We established a multi-tissue gene expression atlas of water buffalo. Furthermore, by comparing them with 4866 cattle transcriptomic data from the cattle genotype-tissue expression atlas (CattleGTEx), we found that the transcriptomes of the two species exhibited conservation in their overall gene expression patterns, tissue-specific gene expression and house-keeping gene expression. We further identified conserved and divergent expression genes between the two species, with the largest number of differentially expressed genes found in the skin, which may be related to structural and functional differences in the skin of the two species. This work provides a source of functional annotation of the buffalo genome and lays the foundations for future genetic and evolutionary studies in water buffalo.