Genome-wide profiling of gene expression and DNA methylation provides insight into low-altitude acclimation in Tibetan pigs

DNA Methylation of Postnatal Liver Development in Pigs

DNA methylation plays an important role in the development and tissue differentiation of eukaryotes. In this study, bisulfite sequencing (BS-seq) technology was used to analyze the DNA methylation profiles of liver tissues taken from Rongchang pigs at three postnatal feeding stages, including newborn, suckling, and adult. The DNA methylation pattern across the genomes or genic region showed little difference between the three stages. We observed 419 differentially methylated regions (DMRs) in promoters, corresponding to 323 genes between newborn and suckling stages, in addition to 288 DMRs, corresponding to 134 genes, between suckling and adult stages and 351 DMRs, corresponding to 293 genes, between newborn and adult stages. These genes with DMRs were mainly enriched in metabolic, immune-related functional processes. Correlation analysis showed that the methylation level of gene promoters was significantly negatively correlated with gene expression. Further, we found that genes related to nutritional metabolism, e.g., carbohydrate metabolism (FAHD1 and GUSB) or fatty acid metabolism (LPIN1 and ACOX2), lost DNA methylation in their promoter, with mRNA expression increased in newborn pigs compared with those in the suckling stage. A few fatty acid metabolism-related genes (SLC27A5, ACOX2) were hypomethylated and highly expressed in the newborn stage, which might satisfy the nutritional requirements of Rongchang pigs with high neonatal birth rates. In the adult stage, HMGCS2-which is related to fatty acid β-oxidation-was hypomethylated and highly expressed, which explains that the characteristics of high energy utilization in adult Rongchang pigs and their immune-related genes (CD68, STAT2) may be related to the establishment of liver immunity. This study provides a comprehensive analysis of genome-wide DNA methylation patterns in pig liver postnatal development and growth. Our findings will serve as a valuable resource in hepatic metabolic studies and the agricultural food industry.

Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing

Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.

Genome-wide comparison of DNA methylation patterns between yak and three cattle strains and their potential association with mRNA transcription

Yak has evolved specific adaptative mechanisms to high-altitude environment. Up to date, only a few studies reported the DNA methylation in yak. In the present study, genome-wide DNA methylome and transcriptome profiles in lung, mammary, and biceps brachii muscle tissues were compared between yak and three cattle breeds (Tibetan cattle, Sanjiang cattle, and Holstein cattle). The association between differentially expressed genes (DEGs) and differentially methylated regions (DMRs) was analyzed, and the biological functions of DEGs potentially driven by DMRs were explored by KEGG enrichment analysis. Finally, we found that yak-specific DMRs-driven DEGs were mainly involved in neuromodulation, respiration, lung development, blood pressure regulation, cardiovascular protection, energy metabolism, DNA repair, and immune functions. The higher levels of the key genes associated with these functions were observed in yak than in cattle, suggesting that DNA methylation might regulate these genes. Overall, the present study contributes basic data at the DNA methylation level to further understand the physiological metabolism in yak.

Plateau Adaptation Gene Analyses Reveal Transcriptomic, Proteomic, and Dual Omics Expression in the Lung Tissues of Tibetan and Yorkshire Pigs

Elevated environments such as plateaus are often classified as low oxygen environments. The hypoxic adaptation mechanisms utilized by organisms in these conditions are not well understood. To address this, the differentially expressed genes (DEGs) involved in hypoxia adaptation were assessed using two pig breeds (Tibetan pig [TP] and Yorkshire sow [YY]). Genes related to lung tissue responses to hypoxia were assessed using transcriptomic (using RNA-seq) and proteomic (using iTRAQ) analysis. A total of 1021 DEGs were screened out. In the iTRAQ omics data, a total of 22,100 peptides were obtained and 4518 proteins were found after filtering. A total of 271 differentially expressed proteins [DEPs] were screened using the conditions of p < 0.05; FC ≤ 0.833; and FC ≥ 1.2. A total of 14 DEGs at the mRNA and protein levels were identified and found to be associated with regulation of the inflammatory response; blood particles; and MAPK cascade response regulation. Among the DEGs, six were associated with hypoxia adaptation function (mitochondria and glycolysis) in pigs. The results of this study identify novel candidate genes involved in porcine hypoxia adaptation mechanisms.

Critical roles of the activation of ethylene pathway genes mediated by DNA demethylation in Arabidopsis hyperhydricity

Hyperhydricity (HH) often occurs in plant tissue culture, seriously influencing the commercial micropropagation and genetic improvement. DNA methylation has been studied for its function in plant development and stress responses. However, its potential role in HH is unknown. In this study, we report the first comparative DNA methylome analysis of normal and hyperhydric Arabidopsis thaliana (L.) Heynh. seedlings using whole-genome bisulfite sequencing (BS-seq). We found that the global methylation level decreased in hyperhydric seedlings, and most of the differentially methylated genes were CHH hypomethylated genes. Moreover, the bisulfite sequencing results showed that hyperhydric seedlings displayed CHH demethylation patterns in the promoter of the ACS1 and ETR1 genes, resulting in upregulated expression of both genes and increased ethylene accumulation. Furthermore, hyperhydric seedling displayed reduced stomatal aperture accompanied by decreased water loss and increased phosphorylation of aquaporins accompanied by increased water uptake. While silver nitrate (AgNO₃ ) prevented HH by maintained the degree of methylation in the promoter regions of ACS1 and ETR1 and downregulated the transcription of both genes. AgNO₃ also reduced the content of ethylene together with the phosphorylation of aquaporins and water uptake. Taken together, this study suggested that DNA demethylation is a key switch that activates ethylene pathway genes to enable ethylene synthesis and signal transduction, which may subsequently influence aquaporin phosphorylation and stomatal aperture, eventually causing HH; thus, DNA demethylation plays a crucial role in HH. These results provide insights into the epigenetic regulation mechanism of HH and confirm the role of ethylene and AgNO₃ in hyperhydricity control.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

DNA methylomic homogeneity and heterogeneity in muscles and testes throughout pig adulthood

DNA methylome pattern is significantly different among tissues, ages, breeds, and genders. We assessed 20 methylome and transcriptome data in longissimus dorsi (LD) or testicles from Bamaxiang (BMX) and Large White pigs (LW) by deep sequencing technology. We identified ~55.7M CpGs and 5.30M, 0.20M, 1.20M, and 0.16M differential CpGs (P<0.01) between tissues, ages, breeds, and genders, respectively. Interestingly, 7.54% of differentially methylated regions (DMRs) are co-localized with promoters, which potentially regulate gene expression. RNA-seq analysis revealed that 23.42% CpGs are significantly correlated with gene expression (mean |r|=0.58, P<0.01), most of which are enriched in tissue-specific functions. Specially, we also found that the methylation levels in promoters of 655 genes were strongly associated with their expression levels (mean |r|=0.66, P<0.01). In addition, differentially methylated CpGs (DMCpGs) between breeds in HOXC gene cluster imply important regulatory roles in myocytes hypertrophy and intermuscular fat (IMF) deposition. Dramatically, higher similarity of methylation pattern was observed within pedigree than across pedigrees, which indicates the existence of heritable methylation regions. In summary, a part of CpGs in promoter can change its methylation pattern and play a marked regulatory function in different physiological or natural environments.

Transcriptomic comparison of liver tissue between Anqing six-end-white pigs and Yorkshire pigs based on RNA sequencing

Chinese indigenous pig and Western commercial pig breeds show different patterns of lipid metabolism, fat deposition, and fatty acid composition; for these reasons, they have become vitally important models of energy metabolism and obesity in humans. To compare the mechanisms underlying lipid metabolism between Yorkshire pigs (lean type) and Anqing six-end-white pigs (obese type), the liver transcriptomes of six castrated boars with a body weight of approximately 100 kg (three Yorkshire and three Anqing) were analyzed by RNA-seq. The total number of reads produced for each liver sample ranged from 47.05 to 62.6 million. Among 362 differentially expressed genes, 142 were up-regulated and 220 were down-regulated in Anqing six-end-white pigs. Based on these data, 79 GO terms were significantly enriched. The top 10 (the 10 with lowest corrected P-value) significantly enriched GO terms were identified, including lipid metabolic process and carboxylic acid metabolic process. Pathway analysis revealed three significantly enriched KEGG pathways including PPAR signaling pathway, steroid hormone biosynthesis, and retinol metabolism. Based on protein–protein interaction networks, multiple genes responsible for lipid metabolism were identified, such as PCK1, PPARA, and CYP7A1, and these were considered promising candidate genes that could affect porcine liver lipid metabolism and fat deposition. Our results provide abundant transcriptomic information that will be useful for animal breeding and biomedical research.

In-depth mapping of the proteome of Tibetan pig tenderloin (longissimus dorsi) using offline high-pH reversed-phase fractionation and LC-MS/MS

In recent years, Tibetan pig breeding and meat processing have developed rapidly. However, the basic physiological and biochemical characteristics of Tibetan pork have not been systematically explored. The present study conducted a high-throughput analysis of the tenderloin (longissimus dorsi) proteome of the Tibetan pigs and performed a functional annotation and bioinformatics analysis of the identified proteins. Based on offline two-dimensional liquid chromatography fractionation and MS/MS identification, a total of 1,723 proteins were identified in the tenderloin of Tibetan pigs. Gene ontology analysis and pathway enrichment analysis revealed that the proteins involved in respiration (oxidative phosphorylation, glycolysis/gluconeogenesis, citric acid cycle, and pyruvate metabolism) and protein synthesis and metabolism (proteasome, amino acid biosynthesis, endoplasmic reticulum protein processing, and ribosomes) were significantly enriched, indicating that the energy production and protein metabolism are the most important physiological processes in Tibetan pig tenderloin. Practical applications The in-depth mapping of the tenderloin (longissimus dorsi) proteome of the Tibetan pigs gives a panoramic perspective at the protein molecular level and provides important information on the mechanisms of postmortem muscle physiology and meat quality formation. Furthermore, the development of Tibetan pork storage and processing technologies would also benefit from the characterization of the biochemical properties of Tibetan pork.

Whole-genome analysis identifying candidate genes of altitude adaptive ecological thresholds in yak populations

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry on the Qinghai-Tibet Plateau. Long-term domestication and natural selection have led to a wide distribution of yak, forming many ecological populations to adapt to the local ecological environment. High altitude is closely related to oxygen density, and it is an important environmental ecological factor for biological survival and livestock production. The aim of the present study was to perform a preliminary analysis to identify the candidate genes of altitude distribution adapted ecological thresholds in yak using next-generation sequence technology. A total of 15,762,829 SNPs were obtained from 29 yaks with high- and low-altitude distribution by genome-wide sequencing. According to the results of the selective sweep analysis with FST and ZHp, 21 candidate genes were identified. 14 genes (serine/threonine protein kinase TNNI3K, TEN1, DYM, ITPR1, ZC4H2, KNTC1, ADGRB3, CLYBL, TANGO6, ASCC3, KLHL3, PDE4D, DEPDC1B and AGBL4) were grouped into 32 Gene Ontology terms, and four genes (RPS6KA6, ITPR1, GNAO1 and PDE4D) annotated in 35 pathways, including seven environmental information processing and one environmental adaptation. Therefore, the novel candidate genes found in the current study do not only support new theories about high-altitude adaptation, but also further explain the molecular mechanisms of altitude adaptation threshold in yaks.

Molecular identification and epidemiological comparison of Cryptosporidium spp. among different pig breeds in Tibet and Henan, China

Background

Cryptosporidium spp. are important zoonotic pathogens infecting a wide range of vertebrate hosts, and causing moderate to severe diarrhea in humans. Cryptosporidium infections are frequently reported in humans and animals worldwide, but little research has been done on local pig breeds such as Tibetan pigs and Yunan Black pigs and imported pig breeds such as Landrace pigs in China. Therefore, a total of 1089 pig fecal samples from four intensive farms in four areas of China, including Tibetan pigs from Gongbujiangda County (n = 180) and Mainling County (n = 434), Tibet, Yunan Black pigs from Sanmenxia, Henan Province (n = 246), and Landrace pigs from Kaifeng, Henan Province (n = 229), and were screened for the presence of Cryptosporidium with microscopy and nested PCR amplification of the small subunit rRNA gene.

Results

The total infection rate of Cryptosporidium in 1089 fecal samples of three different pig breeds was 2.11% (23/1089), and the infection rates of Tibetan pigs, Yunan Black pigs, and Landrace pigs were 0.49% (3/614), 0.41% (1/246), and 8.30% (19/229), respectively. The prevalence of Cryptosporidium infection was significantly higher in weaned piglets (1–2 months) (4.36%, 21/482) than in younger and older age groups (p < 0.01). Sequence analysis of positive samples revealed that there was no mixed infection in our study population, which included 12 cases of C. suis mono-infections (52.17%, 12/23) and 11 cases of C. scrofarum mono-infections (47.83%, 11/23). C. suis was identified in one pre-weaned piglet (< 1 month) and 11 weaned piglets (1–2 months), while C. scrofarum was only detected in 10 weaned piglets (1–2 months) and one finished pig (> 2 months).

Conclusions

This is the first report on the identification of Cryptosporidium spp. in Tibetan pigs, and our findings also elucidate the occurrence and distribution of Cryptosporidium in three different pig breeds in Tibet and Henan, China. More molecular epidemiological studies are required to better clarify the prevalence and public health significance of Cryptosporidium in different pigs.

Association Between the Methylation Statuses at CpG Sites in the Promoter Region of the SLCO1B3, RNA Expression and Color Change in Blue Eggshells in Lushi Chickens

The formation mechanism underlying the blue eggshell characteristic has been discovered in birds, and SLCO1B3 is the key gene that regulates the blue eggshell color. Insertion of an endogenous retrovirus, EAV-HP, in the SLCO1B3 5′ flanking region promotes SLCO1B3 expression in the chicken shell gland, and this expression causes bile salts to enter the shell gland, where biliverdin is secreted into the eggshell, forming a blue shell. However, at different laying stages of the same group of chickens, the color of the eggshell can vary widely, and the molecular mechanism underlying the eggshell color change remains unknown. Therefore, to reveal the molecular mechanism of the blue eggshell color variations, we analyzed the change in the eggshell color during the laying period. The results indicated that the eggshell color in Lushi chickens can be divided into three stages: 20–25 weeks for dark blue, 26–45 weeks for medium blue, and 46–60 weeks for light blue. We further investigated the expression and methylation levels of the SLCO1B3 gene at eight different weeks, finding that the relative expression of SLCO1B3 was significantly higher at 25 and 30 weeks than at other laying weeks. Furthermore, the overall methylation rate of the SLCO1B3 gene in Lushi chickens increased gradually with increasing weeks of egg production, as shown by bisulfite sequencing PCR. Pearson correlation analysis showed that methylation of the promoter region of SLCO1B3 was significantly negatively correlated with both SLCO1B3 expression in the shell gland tissue and eggshell color. In addition, we predicted that CpG5 and CpG8 may be key sites for regulating SLCO1B3 gene transcription. Our findings show that as the level of methylation increases, methylation of the CpG5 and CpG8 sites hinders the binding of transcription factors to the promoter, reducing SLCO1B3 expression during the late period and resulting in a lighter eggshell color.

Genome-wide DNA methylation profiles in Tibetan and Yorkshire pigs under high-altitude hypoxia

Background

Tibetan pigs, which inhabit the Tibetan Plateau, exhibit distinct phenotypic and physiological characteristics from those of lowland pigs and have adapted well to the extreme conditions at high altitude. However, the genetic and epigenetic mechanisms of hypoxic adaptation in animals remain unclear.

Methods

Whole-genome DNA methylation data were generated for heart tissues of Tibetan pigs grown in the highland (TH, n = 4) and lowland (TL, n = 4), as well as Yorkshire pigs grown in the highland (YH, n = 4) and lowland (YL, n = 4), using methylated DNA immunoprecipitation sequencing.

Results

We obtained 480 million reads and detected 280679, 287224, 259066, and 332078 methylation enrichment peaks in TH, YH, TL, and YL, respectively. Pairwise TH vs. YH, TL vs. YL, TH vs. TL, and YH vs. YL comparisons revealed 6829, 11997, 2828, and 1286 differentially methylated regions (DMRs), respectively. These DMRs contained 384, 619, 192, and 92 differentially methylated genes (DMGs), respectively. DMGs that were enriched in the hypoxia-inducible factor 1 signaling pathway and pathways involved in cancer and hypoxia-related processes were considered to be important candidate genes for high-altitude adaptation in Tibetan pigs.

Conclusions

This study elucidates the molecular and epigenetic mechanisms involved in hypoxic adaptation in pigs and may help further understand human hypoxia-related diseases.

Analysis of mitochondrial DNA sequence and copy number variation across five high-altitude species and their low-altitude relatives

High-altitude inhospitable environments impose a formidable life challenge for the local animals. Training and exposure to high-altitude environments produce both distinct physiological and phenotypic characteristics. The mitochondrion, an organelle crucial for the energy production, plays an important role in hypoxia adaptation. In this study, we investigated the mitochondrial DNA (mtDNA) polymorphism and copy number variation between the population pairs from distinct altitudes across the multi-species. Higher mitochondrial DNA control region's genetic diversity is conspicuous in high-altitude animals versus low-altitude relatives. We also found an accordant decrease of mtDNA copy number in most of the tissues from high-altitude animals. Compared to mammals, chickens have significantly distinct mitogenomic characteristics, and more significant changes in the skeletal muscle mtDNA copy number between high- and low-altitude individuals. Our study catches a snapshot of the biological similarities and differences in the mitochondrial high-altitude acclimation across the species.

Identification of a novel antisense long non-coding RNA PLA2G16-AS that regulates the expression of PLA2G16 in pigs

Natural antisense transcripts (NATs) are widely present in mammalian genomes and act as pivotal regulator molecules to control gene expression. However, studies on the NATs of pigs are relatively rare. Here, we identified a novel antisense transcript, designated PLA2G16-AS, transcribed from the phospholipase A2 group XVI locus (PLA2G16) in the porcine genome, which is a well-known regulatory molecule of fat deposition. PLA2G16-AS and PLA2G16 were dominantly expressed in porcine adipose tissue, and were differentially expressed between Tibetan pigs and Rongchang pigs. In addition, PLA2G16-AS has a weak sequence conservation among different vertebrates. PLA2G16-AS was also shown to form an RNA-RNA duplex with PLA2G16, and to regulate PLA2G16 expression at the mRNA level. Moreover, the overexpression of PLA2G16-AS increased the stability of PLA2G16 mRNA in porcine cells. We envision that our findings of a NAT for a regulatory gene associated with lipolysis might further our understanding of the molecular regulation of fat deposition.