Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production

Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production Jinwei Xin, Zhixin Chai, Chengfu Zhang, Qiang Zhang, Yong Zhu, and 6 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-20775/v4 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 20 Oct, 2020 Read the published version in BMC Genomics → Version 4 posted 4 You are reading this latest preprint version Show more versions Abstract Background Domestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), by providing milk and meat. They have evolved numerous physiological adaptations to high-altitude environment, including strong blood oxygen transportation capabilities and high metabolism. The roles of DNA methylation and gene expression in milk production and high-altitudes adaptation need further exploration. Results We performed genome-wide DNA methylome and transcriptome analyses of breast, lung, and biceps brachii muscle tissues from yaks of different ages. We identified 432,350 differentially methylated regions (DMRs) across the age groups within each tissue. The post-mature breast tissue had considerably more differentially methylated regions (155,957) than that from the three younger age groups. Hypomethylated genes with high expression levels might regulate milk production by influencing protein processing in the endoplasmic reticulum. According to weighted gene correlation network analysis, the “hub” gene ZGPAT was highly expressed in the post-mature breast tissue, indicating that it potentially regulates the transcription of 280 genes that influence protein synthesis, processing, and secretion. The tissue network analysis indicated that high expression of HIF1A regulates energy metabolism in the lung. Conclusions This study provides a basis for understanding the epigenetic mechanisms underlying milk production in yaks, and the results offer insight to breeding programs aimed at improving milk production. Epigenetics & Genomics Animal Science milk production DNA methylation transcription factor epigenetic regulation Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Background Domestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), in the Himalayas, and in the connecting Central Asian highlands. They provide milk, meat, hides, fiber, fuel, and transportation [1, 2]. Milk is an important source of high-quality protein. It contains high quantities of essential amino acids, such as lysine, which is commonly deficient in many human diets [3]. Milk proteins also impact immunomodulatory processes and gastrointestinal activities [4]. Protein content and composition influence the technological properties of milk and are therefore important for the dairy industry, particularly in Europe, where the majority of the milk produced is used to produce cheese. In recent decades, there have been extraordinary advances in our knowledge of the physiology and biochemistry of the lactating mammary gland. Previous research indicates that milk protein synthesis in the mammary gland depends on hormonal and developmental cues that modulate the transcriptional and translational regulation of genes through the activity of specific transcription factors, non-coding RNAs, and alterations of the chromatin structure in mammary epithelial cells [5]. The interplay between these factors may influence milk protein synthesis, which is crucial during the entire lactation process in high-producing dairy cattle. Despite such advancements in research, little is currently known about the physiological and cellular regulation required for milk protein synthesis and secretion in yak. We hypothesized that the genes responsible for milk production were regulated by DNA methylation and that distinct sub-modules of correlated expression variation could be identified. In this study, we performed genome-wide DNA methylome and transcriptome analyses of yak lung, breast, and biceps brachii muscle tissues at four different stages of development to identify the regulatory networks associated with milk protein synthesis, metabolism, and secretion. Results Global DNA methylation and gene expression in the breast, lungs, and biceps brachii muscle at different ages We generated the methylomes and transcriptomes of lung, breast, and biceps brachii muscle tissues from 12 female Riwoqe yaks at four different stages (n=3/stage) of development: 6 months old (MO) (young), 30 MO (pre-mature), 54 MO (mature) and 90 MO (post mature). Among these, only the post-mature yaks (90 months) were lactation, with ~3.16 kg/day milk yield. After performing sequence quality control and filtering, we obtained a single-base resolution methylome covering 85.6% (27,471,373/32,092,725) of CpG sites across the genome with an average depth of 22.5×. We first calculated pairwise Pearson's correlations of CpG sites with at least 10× coverage depth across all samples, which were well clustered by tissue types (Figure 1a). The correlation of CpG methylation levels for biological replicates was strong (median Pearson's r = 0.74), the correlation of CpG methylation levels between ages (median Pearson's r = 0.72) was relatively weaker and the correlation between tissues was weakest (median Pearson's r = 0.66) (Figure 1c). We aligned the transcriptome sequencing data for all samples to our newly assembled yak genome reference, and we subsequently obtained the transcripts. In total, we obtained 2,0504 transcripts, that were then annotated to the Gene Ontology (GO) [6], InterPro [7], Kyoto Encyclopedia of Genes and Genomes (KEGG) [8], Swiss-Prot [9], and TrEMBL [10] databases (Table S1). We also calculated pairwise Pearson's correlations of all transcripts and obtained similar results to those of DNA methylation (Figure 1b). Biological replicates showed the highest correlation coefficients, while different tissues showed the lowest correlation coefficients (Figure 1d). Differentially methylated regions among the age groups We determined differentially methylated regions (DMRs) across age groups within the breast, lung, and biceps brachii muscle tissues (Table S2-4). Within the lung and biceps brachii muscle tissues, age groups did not differ in age-related DMRs (A-DMRs), but post-mature breast tissue had considerably more DMRs (155,957) than the three younger age groups (Figure 2a). We then investigated the correlations between DMRs and their corresponding differentially expressed genes. The ratio of negatively to positively correlated gene pairs was 1.02 for promoters with DMRs and 0.95 for gene bodies with DMRs. Not every methylation was correlated with the expression of its associated gene, due to the gene regulation complexity [11]. At ~90 months, ~120 days after giving birth for the third time, yaks are in the lactation period (milk yield of ~3.16kg/day)(Li & Jiang, 2019), so it is possible that the observed methylation partially controls yak lactation. Since promoter methylation decreases gene expression [11], we selected 375 hypomethylated promoter genes (highly expressed) along with 207 hypermethylated promoter genes (lowly expressed) from post-mature yak breast tissues. The hypomethylated (highly expressed) genes were only enriched in“protein processing in endoplasmic reticulum (ER)” (9 genes, 2.964-fold enrichment, p = 0.0049). Specifically, the genes are involved in vesicle trafficking (SEC23B), oligosaccharide linking (MOGS, RPN2), folding and assembly (HSPA5), transportation (LMAN2, SEL1L), and ubiquitination and degradation (UBE2J1, UBE2J2, DERL2) [12, 13]. These genes were significantly upregulated at 90 months in breast tissues, but not in lung and biceps brachii muscle tissues (Figure 2e). Based on this data, methylation might help regulate milk production by influencing protein processing in the endoplasmic reticulum during the lactation period. We also examined A-DMRs that overlapped across age groups. Young and post-mature tissues rarely shared A-DMRs when comparing the lung and muscle tissues (for young tissues, muscle: 586 A-DMRs, breast: 2,249 A-DMRs, lung: 496 A-DMRs; for post-mature tissues, muscle: 470 A-DMRs, breast: 12,050 A-DMRs, lung: 772 A-DMRs) (Figure 2b, c). Pre-mature and mature stages also rarely shared A-DMRs across muscle and lung tissues (Figure 2d), suggesting that methylation patterns were already established at the young stage and that no extensively divergent epigenetic difference occurred across different age groups under natural high-altitude conditions. Consensus network analysis for tissues and age groups We first performed a multi-way ANOVA test for each gene across all samples (n=36) to test the null hypothesis that the gene expression level did not differ among age groups and tissues. At the threshold for significance (p<0.05), 417 age-related and 8,560 tissue-related genes were selected for further weighted gene correlation network analysis (WGCNA), which uses network topography to group genes into modules based on correlations [14]. Next, we conducted WGCNA for tissue- and age-related gene expression respectively to identify a “consensus network”–a common pattern of genes that are correlated in all conditions. We performed a consensus network, module statistic, and eigengene network analyses to identify modules, to assess relationships between modules and traits, and to study the relationships between co-expression modules [15]. The consensus networks identified for tissues and age groups had clearly delineated modules (Figure 3a, 3b), and the modules identified were significantly correlated with tissues and age groups (Figure 3c, 3d). Age network analysis indicates that ZGPAT might regulate milk production Within the age-related gene network, the largest module (“turquoise”, n=356) had negative correlation for breast tissue (r=-0.57, p=3e-04) and age (r= 0.37, p=0.03) and a positive correlation for biceps brachii muscle tissues (Table S5). The breast tissue had a stronger signal than that of age and may have overwhelmed the signal from age. Genes in this module were enriched in the GO categories of “protein polyubiquitination,” “RNA polymerase II core promoter proximal region sequence-specific DNA binding,” “ATP binding,” “transcription, DNA-templated,” and “negative regulation of transcription from RNA polymerase II promoter” (p-values of 0.000344, 0.000822, 0.000991, 0.00139, and 0.00244, respectively). The “blue” (n=48) and “grey” (n=13) modules showed only a negative correlation with age (r= -0.58, p= 2e-04; r= -0.76, p= 6e-08, respectively) and exhibited no enrichment of GO categories for genes. After applying the threshold of the absolute value of gene significance for age (|GS| >0.5) and module membership measures (|MM| >0.6) in each module, we defined 20 and 7 “hub” genes in the “turquoise” and “blue” modules (Table 1). The gene expression of the “hub” genes was well clustered by modules, which was consistent with the negative correlation with age (“turquoise” r=0.37, “blue” r=-0.58). The upregulated expression level of the “hubs” in breast tissue at 90 months old indicated that the “turquoise” module had a stronger correlation with breast tissue than with age (Figure 4a). Table 1. List of “hub” genes in the consensus network for age Yak ID Gene symbol Module color Gene significance p-value Module membership p-value BmuPB009868 AP3D1 turquoise 0.513965075 0.001344116 0.868945999 6.37E-12 BmuPB004083 TMEM30A turquoise 0.505722627 0.001652657 0.848054843 6.65E-11 BmuPB000726 DDRGK1 turquoise 0.53596122 0.000754414 0.842070133 1.22E-10 BmuPB019011 OTUD3 turquoise 0.543916749 0.000606215 0.828680521 4.37E-10 BmuPB000091 CNPPD1 turquoise 0.51195228 0.001414359 0.814537409 1.50E-09 BmuPB006815 TOR1AIP1 turquoise 0.544705745 0.000593032 0.809789221 2.21E-09 BmuPB007137 MGAT4A turquoise 0.545969144 0.000572452 0.798170686 5.51E-09 BmuPB007385 HECA turquoise 0.603273223 9.84E-05 0.770629907 3.85E-08 BmuPB019443 SYAP1 turquoise 0.500377628 0.001884511 0.733120646 3.68E-07 BmuPB009825 PHAX turquoise 0.561938984 0.00036184 0.70497442 1.59E-06 BmuPB008032 UBE2G2 turquoise 0.52385538 0.001041698 0.70208134 1.83E-06 BmuPB012517 SYF2 turquoise 0.568328068 0.000299194 0.699367268 2.08E-06 BmuPB001610 FURIN turquoise 0.567906296 0.000303008 0.697046314 2.32E-06 BmuPB000679 ZGPAT turquoise 0.504624296 0.001698141 0.683958391 4.25E-06 BmuPB007049 SKP2 turquoise -0.527608908 0.000943732 -0.669790419 7.91E-06 BmuPB000224 C2orf6 turquoise 0.531191165 0.000857925 0.652871802 1.59E-05 BmuPB007420 TAB2 turquoise 0.505905132 0.001645204 0.652844304 1.59E-05 BmuPB018569 ORAOV1 turquoise 0.51035386 0.001472421 0.637043276 2.95E-05 BmuPB010550 CCPG1 turquoise 0.644779605 2.19E-05 0.612874899 7.08E-05 BmuPB012521 TMEM57 turquoise 0.57317319 0.000258349 0.60967979 7.91E-05 BmuPB013324 MCM3 blue -0.583585649 0.000186982 0.841505398 1.29E-10 BmuPB010064 SPC24 blue -0.509962181 0.001486965 0.744497748 1.93E-07 BmuPB003102 C17orf49 blue -0.526803776 0.000964031 0.741805787 2.26E-07 BmuPB015902 SERPINH1 blue -0.607769763 8.44E-05 0.721017606 7.04E-07 BmuPB016582 SRPX2 blue -0.51838468 0.001200558 0.698151166 2.20E-06 BmuPB012996 UCK2 blue -0.588274454 0.000161067 0.677420455 5.68E-06 BmuPB015372 UMPS blue -0.503577413 0.001742514 0.648111658 1.92E-05 We used the AnimalTFDB 3.0 database [16] to examine transcription factors in these 27 “hubs” and found that ZGPAT encodes a transcription regulator protein and was significantly upregulated in breast tissue at 90 months of age (Figure 4a). Previous study reported that this protein specifically binds the 5'-GGAG[GA]A[GA]A-3' consensus sequence and represses transcription by recruiting the chromatin multi-complex NuRD to target promoters [17]. ZGPAT was highly expressed in breast tissue at 90 months, and it potentially regulated the transcription of 280 genes (weight >0.15) in the network from the “turquoise” module. In order to identify the most important cellular activities controlled by this TF regulatory network, we analyzed over-represented GO biological process and molecular function terms, as well as KEGG pathways. These potential target genes were enriched in the GO categories of “protein binding,” “ATP binding,” and “zinc ion binding,” among others, and the KEGG categories of “aminoacyl-tRNA biosynthesis,” “autophagy animal,” and “protein processing in endoplasmic reticulum” (Figure 4b). These enriched GO terms and KEGG pathways likely help regulate protein synthesis, processing, and secretion in breast tissue. For example, 6 of 7 genes from the “protein processing in endoplasmic reticulum” category were also upregulated at 90 months of age in breast tissue (Figure 4c) and involved in multiple processes in the endoplasmic reticulum, including vesicle trafficking (SEC24C), folding and assembly (SELENOS), transportation (BCAP31), and ubiquitination and degradation (BAG1, UBE2G2, and MARCH6) [12, 13]. Only DNAJC10 was downregulated at 90 months of age in breast tissue. This gene encodes an endoplasmic reticulum co-chaperone that is part of the endoplasmic reticulum-associated degradation complex involved in recognizing and degrading misfolded proteins [13]. Tissue network analysis indicates regulative role of HIF1A in lung Within the tissue-related gene module network, four modules showed a positive correlation and two showed a negative correlation with the lung; all significant module-trait relationships were negative in the muscle but positive in the breast (Figure 3d). Moreover, 99.54% of the total 8,560 tissue-related genes were related to the top 4 modules (“turquoise,” n=3833; “blue,” n=2795; “brown,” n=1052; “yellow,” n=339) (Figure 5a, Table S6), and these modules were also highly correlated with other modules; for example, “brown,” “yellow,” and “black” showed a high eigengene adjacency with each other (Figure 5b). We applied the more stringent threshold absolute value of gene significance for the age and module membership measures in the top four modules to identify “hub” genes in the “turquoise,” “blue,” “brown,” and “yellow” modules. With the threshold values of |GS|>0.7 and |MM| >0.8, 34 “hub” genes were identified in the “turquoise” module. Twenty-four hubs were then filtered from the gene significance of module-lung relationships, and 10 were filtered from the gene significance of module-breast relationships. These were further divided into 3 clusters by hierarchical clustering, which showed high expression levels in the breast (cluster 1), lung (cluster 2), and biceps brachii muscle (cluster 3) tissues, with distinct clustering patterns by tissue (Figure 5c). Table 2. List of “hub” genes in the consensus network for tissue Yak ID Gene symbol Module color Tissue Gene significance p-value Module membership p-value BmuPB014336 EEF1G turquoise lung -0.73907039 2.64E-07 0.826748126 5.20E-10 BmuPB017352 PMS1 turquoise lung -0.71599797 9.13E-07 0.850552958 5.12E-11 BmuPB000878 MTPAP turquoise lung -0.715583465 9.33E-07 0.912877613 8.73E-15 BmuPB018762 CXHXorf58 turquoise lung -0.709433358 1.27E-06 0.82262518 7.50E-10 BmuPB014450 RWDD4 turquoise lung -0.708526806 1.33E-06 0.923713618 9.94E-16 BmuPB011005 HUS1 turquoise lung -0.707137304 1.43E-06 0.875052199 2.97E-12 BmuPB005540 MTUS1 turquoise lung -0.704482832 1.62E-06 0.871626441 4.58E-12 BmuPB011453 EBF3 turquoise lung -0.703373458 1.71E-06 0.870717052 5.13E-12 BmuPB010608 LAMB3 turquoise lung 0.70783756 1.38E-06 -0.85067459 5.05E-11 BmuPB004299 C3H3orf58 turquoise lung 0.708872549 1.31E-06 -0.88522675 7.61E-13 BmuPB020871 G6PD turquoise lung 0.708930411 1.30E-06 -0.90727231 2.41E-14 BmuPB007438 ZCCHC6 turquoise lung 0.709740967 1.25E-06 -0.84738754 7.13E-11 BmuPB007592 CTDSPL turquoise lung 0.711634321 1.14E-06 -0.8813581 1.30E-12 BmuPB004894 HIF1A turquoise lung 0.713237417 1.05E-06 -0.87302984 3.84E-12 BmuPB020508 FAM122B turquoise lung 0.720139618 7.37E-07 -0.82428698 6.48E-10 BmuPB018102 CCDC82 turquoise lung 0.720592854 7.20E-07 -0.90665968 2.68E-14 BmuPB014539 CNTRL turquoise lung 0.722130497 6.64E-07 -0.87149974 4.65E-12 BmuPB003882 VAV3 turquoise lung 0.725601285 5.53E-07 -0.89496378 1.82E-13 BmuPB020153 EPS8L1 turquoise lung 0.727495147 4.99E-07 -0.82286844 7.34E-10 BmuPB010308 PGM2 turquoise lung 0.746811608 1.69E-07 -0.83795918 1.83E-10 BmuPB004008 GDAP2 turquoise lung 0.754806874 1.05E-07 -0.88661303 6.26E-13 BmuPB012568 WASF2 turquoise lung 0.757428944 8.92E-08 -0.83874427 1.69E-10 BmuPB011966 ACTG1 turquoise lung 0.772114072 3.49E-08 -0.84373776 1.03E-10 BmuPB014173 STAT6 turquoise lung 0.80896084 2.37E-09 -0.84981355 5.53E-11 BmuPB015106 MAN2A1 turquoise breast 0.705124693 1.57E-06 -0.84782348 6.81E-11 BmuPB008614 VPS26A turquoise breast 0.731553581 4.01E-07 -0.88669366 6.18E-13 BmuPB018496 FAM92A turquoise breast 0.718223061 8.14E-07 -0.85598454 2.85E-11 BmuPB005078 RASA1 turquoise breast 0.705902582 1.51E-06 -0.85869193 2.11E-11 BmuPB011476 FAM53B turquoise breast -0.728187435 4.81E-07 0.848680696 6.23E-11 BmuPB008348 EPDR1 turquoise breast -0.710771064 1.19E-06 0.832481139 3.08E-10 BmuPB003453 TROVE2 turquoise breast 0.72654448 5.26E-07 -0.80669368 2.84E-09 BmuPB021200 VWA7 turquoise breast -0.722814498 6.41E-07 0.814033145 1.56E-09 BmuPB010528 RNF111 turquoise breast 0.737234077 2.92E-07 -0.80146106 4.28E-09 BmuPB015811 FRMD3 turquoise breast -0.739223421 2.61E-07 0.809147869 2.33E-09 According to the AmalTFDB 3.0 database [16], EBF3, HIF1A, and STAT6 were annotated as transcription factors. EBF3 encodes a member of the early B-cell factor (EBF) family of DNA binding transcription factors. EBF proteins are involved in B-cell differentiation, bone development, and neurogenesis and they may also function as tumor suppressors [18]. STAT6 is a member of the STAT family of transcription factors. In response to cytokines and growth factors, STAT family members are phosphorylated by the receptor-associated kinases. They form homo- or heterodimers that translocate to the cell nucleus where they act as transcription activators [19]. HIF1A, hypoxia-inducible factor-1, functions as a master regulator of the cellular and systemic homeostatic response to hypoxia by activating the transcription of genes, involved in energy metabolism, angiogenesis, and apoptosis, as well as other genes whose protein products increase oxygen delivery or facilitate metabolic adaptation to hypoxia [20]. HIF1A was upregulated in the lung to potentially control the transcription of 2008 genes (weight >0.15), which were enriched in multiple GO biological process categories, molecular function categories and KEGG pathways. Most of the enriched GO terms and KEGG pathways were related to energy metabolism. Specifically, enriched GO terms included “mitochondrial respiratory chain complex I assembly,” “NADH dehydrogenase (ubiquinone) activity,” “ATP binding,” “mitochondrial translation,” “tricarboxylic acid cycle,” and “GTP binding” while enriched KEGG pathways included “thermogenesis,” “carbon metabolism,” and “citrate cycle TCA cycle” (Figure 5d). Mitochondria function as the primary energy producers of the cell and serves as the hub for a variety of immune pathways as the center of biosynthesis, oxidative stress response, and cellular signaling [21]. NADH dehydrogenase is a core subunit of the mitochondrial membrane respiratory chain and is believed to contribute to the minimal assembly required for catalysis [22]. The protein which binds ATP or GTP, carries three phosphate groups esterified to a sugar moiety and provides energy and phosphate sources for the cell [23, 24]. The tricarboxylic acid cycle, a series of metabolic reactions in aerobic cellular respiration, occurs in the mitochondria of animals and plants. During this cycle, acetyl-CoA, which is formed from pyruvate produced during glycolysis, is completely oxidized to CO 2 via interconversion of various carboxylic acids. This results in the reduction of NAD and FAD to NADH and FADH2, respectively, and indirectly in the synthesis of ATP by oxidative phosphorylation [25]. The “thermogenesis” pathway is essential for warm-blooded animals, because it ensures normal cellular and physiological functions under challenging environmental conditions [26]. Discussion Previous studies described transcription profiling of the mammary gland in livestock (including cattle [27], sheep [28], and goat [29]) and DNA methylation profiling of the mammary gland in cattle [30]. These studies showed temporal and spatial specificity in the methylome and transcriptome profiles of the mammary gland in different species, but they only reported the differential gene expression profile of the mammary gland; and the regulatory network is still unknown. In this study, we for the first time generated the methylomes and transcriptomes of lung, breast, and biceps brachii muscle tissues from yak at four different stages of development (6, 30, 54, and 90 months; young, pre-mature, mature, and post-mature, respectively). We found that breast tissue at 90 months showed considerably differential methylation levels compared with other month groups, but lung and biceps brachii muscle tissues did not. Enrichment analysis for upregulated genes with hypomethylated DMRs from breast tissues of 90-month-old yaks showed that DNA methylation might regulate the activation of the “protein processing in endoplasmic reticulum (ER)” pathway. Because only the 90-month-old yaks were in the lactation period, it appears that DNA methylation regulates milk production by influencing protein processing in the endoplasmic reticulum. In this study, hub genes were identified by WGCNA. The data show that the hub genes with the highest MM and GS in modules of interest should be candidates for further research. This study identified turquoise module genes associated with milk yield, and 20 genes were considered hub genes, showing the highest mRNA expression level in breast tissue at 90 months when yaks enter the lactation period. In these hub genes, ZGPAT was annotated as a transcription factor that potentially regulated the transcription of 280 genes in the “turquoise” module network, which was enriched in the KEGG categories of “aminoacyl-tRNA biosynthesis,” “autophagy animal,” and “protein processing in endoplasmic reticulum”. This result suggests that ZGPAT helps with regulating protein synthesis, processing, and secretion in breast tissue. Moreover, the 7 genes potentially regulated by ZGPAT in “protein processing in endoplasmic reticulum” were totally different from the aforementioned 9 genes regulated by hypomethylation, illustrating that DNA methylation and transcription factor possibly co-regulate milk production. In addition, the tissue network analysis confirms the importance of HIF1A in regulating energy metabolism, which is necessary for adaptations to low temperature and hypoxia in high altitudes. Conclusions The results of this comprehensive study provide a solid basis for understanding the roles of DNA methylation and transcriptional network underlying milk protein synthesis and high-altitude adaptation in yaks. This information advances our understanding of the regulatory network in the mammary gland at different developmental stages and may help inform breeding programs aimed at improving milk production. Methods Animals and samples In total, twelve female yaks (belonging to an indigenous yak breed that lives at altitudes of 3800-4000 meters above sea level in Riwoqe, Tibet, China) were collected between June and December of 2016 from private farms. They were grouped with 3 replicates into age categories of 6, 30, 54, and 90 months. At the time of slaughter, their mean live weights were 44.93 kg (6 months old), 153.06 kg (30 months old), 188.3 kg (54 months old), and 243.56 kg (90 months old). Only the 90 month-old yaks were lactating, producing ~3.16 kg/day milk yield (~120 days after giving birth for the third time). The 54-month-old yaks were in a dry period (~540 days after giving birth for the first time) [31]. The blood relation between the last three generations, which were housed simultaneously and fed the same diets, was unknown. The yaks were not fed the night before they were slaughtered and were humanely sacrificed with the following procedures: (1) showered with clean water close to body temperature (35-38°C), (2) electrically stunned (120V dc, 12s) prior to exsanguination, (3) sacrificed while in the coma by bloodletting from carotid artery and jugular vein, (4) and dissected rapidly to obtain breast, lung, and biceps brachii muscle tissue samples, which were immediately frozen in liquid nitrogen, and stored at −80°C until RNA and DNA extraction. Whole genome bisulfite sequencing The QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) was used to isolate the high-quality DNA from each sample. According to the manufacturer’s instructions, 1 μg of genomic DNA was fragmented by sonication to a mean size of approximately 250 bp and subsequently used for whole genome bisulfite sequencing (WGBS) library construction using an Acegen Bisulfite-Seq Library Prep Kit (Acegen, Shenzhen, GD, China). Briefly, fragmented DNA was end-repaired, 5'-phosphorylated, 3'-dA-tailed, and then ligated to methylated adapters. The methylated adapter-ligated DNAs were purified using 1× Agencourt AMPure XP magnetic beads (Beckman Coulter, Brea, CA, USA) and subjected to bisulfite conversion with a ZYMO EZ DNA Methylation-Gold Kit (Zymo research, Irvine, CA, USA). The converted DNAs were then amplified using 25 μl HiFi HotStart U+ RM and 8-bp index primers with a final concentration of 1 μM each. The constructed WGBS libraries were then analyzed with an Agilent 2100 Bioanalyzer (Agilent Technologies, SantaClara, CA, USA), quantified with a Qubit fluorometer with Quant-iT dsDNA HS Assay Kit (Invitrogen, Carlsbad, CA,USA), and finally sequenced on an Illumina Hiseq X ten sequencer (PE150 mode) (Illumina, San Diego, CA, USA) Methylation calculation and identification of DMRs Low quality reads that contained more than 5 ‘N’s or had a low-quality value for over 50% of the sequence (Phred score < 5) were filtered. The sequencing reads of the samples were aligned to the yak reference genome [32] using BSMAP (Version 2.74) [33]. The methylated CpG (mCG) sites were identified following a previously described algorithm [34]. The methylation levels for each sample were calculated using in-house Perl scripts. Differentially methylated regions (DMRs) were identified using metilene (Version 0.2-6) within a 500 bp sliding window at 250 bp steps with at least 10 CpGs covered by over 10× sequence reads, applying the thresholds of differential methylation β >=15%, FDR for two-dimensional Kolmogorov-Smirnov-Test p-value <0.05. [35]. The enrichment analyses were conducted using WebGestalt (WEB-based Gene SeT AnaLysis Toolkit) [36]. Total RNA extraction, library preparation, and sequencing The TRIzol reagent (Invitrogen, Carlsbad, CA, USA) was used to isolate the total RNA of each sample. The purity, concentration, and integrity of RNA were checked using the NanoPhotometer spectrophotometer (IMPLEN, Westlake Village, CA, USA), the Qubit RNA Assay Kit in Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), and the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 System (Agilent Technologies, SantaClara, CA, USA), respectively. For each sample, 3 µg high-quality RNA was used as input material for RNA-seq library preparation. First, ribosomal RNA was removed using the Epicentre Ribo-Zero rRNA Removal Kit (Epicentre, Madison, WI, USA). Next, the rRNA-depleted RNA was used to create sequencing libraries using the NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). Finally, the library products were purified using 1× Agencourt AMPure XP magnetic beads (Beckman Coulter, Brea, CA, USA) and the Agilent Bioanalyzer 2100 System (Agilent Technologies, SantaClara, CA, USA) was employed to assess the library quality. Clustering of the index-coded samples was completed on a cBot Cluster Generation System using the TruSeq PE Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA), and then the libraries were sequenced on the Illumina HiSeq X Ten Platform to generate 150 bp paired-end reads. Quality analysis, transcriptome assembly, and abundance estimation Clean reads were obtained by removing reads containing the adapter or poly-N and by removing low quality reads (over 10% of the sequence with a quality value < 30) from the raw data using in-house Perl scripts. All downstream analyses were based on the good-quality clean reads. Paired-end clean reads were mapped to the yak reference genome [32] with STAR (available at https://github.com/alexdobin/STAR/releases). The mapped reads of each sample were assembled using StringTie [37]. Next, all sample transcriptomes were merged to reconstruct a comprehensive transcriptome using Perl scripts. After the final transcriptome was generated, StringTie and edgeR were used to estimate the expression levels of all transcripts [38]. StringTie was used to assess the expression level of mRNAs by calculating fragments per kilobase of transcript per million fragments mapped (FPKM). Differentially expressed mRNAs were identified using the DESeq2 package, with the criteria of fold-change log2 > 1 or log2 < -1 and with the statistical significance set to FDR < 0.05. Weighted gene correlation network analysis A WGCNA network [14] was generated for both age-related and tissue-related genes. Consensus networks and module statistics followed the overall approach described by Langfelder et al. (2008). The network was derived based on a signed Spearman correlation using a b of 10 as a weight function. The topological overlap metric (TOM) [15] was derived from the resulting adjacency matrix and used to cluster the modules using the blockwiseModules function (blockwise Consensus Modules, for the consensus modules) and the dynamic tree cut algorithm [15] with a height of 0.25 and a deep split level of 2, a reassign threshold of 0.2, and a minimum module size of 30 (100 for the consensus network). Each eigenmodules–the first principal component of the module and a signature of gene expression–was then correlated with the dose, and each module that was correlated with the dose-response curve with a p-value < 0.01 (p-value < 0.05 for the consensus network) was considered statistically significant. List Of Abbreviations QTP: Qinghai-Tibetan Plateau DMRs: differentially methylated regions WGBS: whole genome bisulfite sequencing FPKM: fragments per kilobase of transcript per million fragments WGCNA: Weighted gene correlation network analysis TOM: topological overlap metric GO: Gene Ontology KEGG: Kyoto Encyclopedia of Genes and Genomes A-DMRs: age-related DMRs ER: endoplasmic reticulum GS: gene significance MM: module membership measures Declarations Ethics approval and consent to participate All protocols for collection of the semen samples of yaks were reviewed and approved by the Ethics Committee at Institute of Animal Science and Veterinary, Tibet Academy of Agricultural and Animal Husbandry Sciences (Permit Number: 2015-216). Consent for publication Not applicable. Availability of data and material The DNA methylation data and RNA transcriptome data in this study are available in SRA under the accession numbers PRJNA530286 and PRJNA512958, respectively. Competing interests The authors declare that they have no competing interests. Funding This work was supported by a program of Provincial Department of Finance of the Tibet Autonomous Region (No: XZNKY-2019-C-052), Program National Beef Cattle and Yak Industrial Technology System (No: CARS-37), Basic Research Programs of Sichuan Province (No: 2019YJ0256) and the Open Project Program of State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement (NO:XZNKY-2019-C-007K10). The funding bodies had no role in the study design; collection, analysis, and interpretation of data; or in writing the manuscript. Authors' contributions JX, QJ and JZ planned and coordinated the study and wrote the manuscript. CY, XC and HJ collected the samples. ZC and CZ performed the library construction and sequencing and the quality control analysis. QZ, YZ and HC performed downstream analysis of the data and assisted in the generation of additional files for the manuscript. All authors read and approved the final manuscript. Acknowledgements The authors thank the animal husbandry station of Chang Du and the agricultural bureau of Riwoqe for all their support. References Wiener G, Han J-L, Long R-J: The Yak , 2nd edn. Bangkok, Thailand: Regional Office for Asia and the Pacific, Food and Agriculture Organization of the United Nations; 2003. Qiu Q, Zhang G, Ma T, Qian W, Wang J, Ye Z, Cao C, Hu Q, Kim J, Larkin DM: The yak genome and adaptation to life at high altitude . Nature Genetics 2012, 44 (8):946-949. Nichols K, Doelman J, Kim JJM, Carson M, Metcalf JA, Cant JP: Exogenous essential amino acids stimulate an adaptive unfolded protein response in the mammary glands of lactating cows . J Dairy Sci 2017, 100 (7):5909-5921. Kanwar JR, Kanwar RK, Sun X, Punj V, Matta H, Morley SM, Parratt A, Puri M, Sehgal R: Molecular and biotechnological advances in milk proteins in relation to human health. Curr Protein Pept Sci 2009, 10 (4):308-38. Osorio JS, Lohakare J, Bionaz M: Biosynthesis of milk fat, protein, and lactose: roles of transcriptional and posttranscriptional regulation . Physiol Genomics 2016, 48 (4):231-56. Gene Ontology Consortium: Expansion of the Gene Ontology knowledgebase and resources . Nucleic Acids Res 2017, 45 (D1):D331-D338. Finn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, Chang HY, Dosztányi Z, El-Gebali S, Fraser M, Gough J, Haft D, Holliday GL, Huang H, Huang X, Letunic I, Lopez R, Lu S, Marchler-Bauer A, Mi H, Mistry J, Natale DA, Necci M, Nuka G, Orengo CA, Park Y, Pesseat S, Piovesan D, Potter SC, Rawlings ND, Redaschi N, Richardson L, Rivoire C, Sangrador-Vegas A, Sigrist C, Sillitoe I, Smithers B, Squizzato S, Sutton G, Thanki N, Thomas PD, Tosatto SC, Wu CH, Xenarios I, Yeh LS, Young SY, Mitchell AL: InterPro in 2017-beyond protein family and domain annotations . Nucleic Acids Res 2017, 45 (D1):D190-D199. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K: KEGG: new perspectives on genomes, pathways, diseases and drugs . Nucleic Acids Res 2017, 45 (D1):D353-D361. Boutet E, Lieberherr D, Tognolli M, Schneider M, Bansal P, Bridge AJ, Poux S, Bougueleret L, Xenarios I: UniProtKB/Swiss-Prot, the Manually Annotated Section of the UniProt KnowledgeBase: How to Use the Entry View . Methods Mol Biol 2016, 1374 :23-54. Junker V, Contrino S, Fleischmann W, Hermjakob H, Lang F, Magrane M, Martin MJ, Mitaritonna N, O'Donovan C, Apweiler R: The role SWISS-PROT and TrEMBL play in the genome research environment. J Biotechnol 2000, 78 (3):221-34. Burgess, D.J: Gene expression: Principles of gene regulation across tissues. Nature Reviews Genetics, 2017, 18 (12):701-701. Preston GM, Brodsky JL: The evolving role of ubiquitin modification in endoplasmic reticulum-associated degradation . Biochem J 2017, 474 (4):445-469. McCaffrey K, Braakman I: Protein quality control at the endoplasmic reticulum . Essays Biochem 2016, 60 (2):227-235 Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis . BMC Bioinformatics 2008, 9 :559. Liu X, Hu AX, Zhao JL, Chen FL: Identification of Key Gene Modules in Human Osteosarcoma by Co-Expression Analysis Weighted Gene Co-Expression Network Analysis (WGCNA) . J Cell Biochem 2017, 118 (11):3953-3959. Hu H, Miao YR, Jia LH, Yu QY, Zhang Q, Guo AY: AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors . Nucleic Acids Res 2019, 7 (D1):D33-D38. Sun LD, Xiao FL, Li Y, Zhou WM, Tang HY, Tang XF, Zhang H, Schaarschmidt H, Zuo XB, Foelster-Holst R, He SM, Shi M, Liu Q, Lv YM, Chen XL, Zhu KJ, Guo YF, Hu DY, Li M, Li M, Zhang YH, Zhang X, Tang JP, Guo BR, Wang H, Liu Y, Zou XY, Zhou FS, Liu XY, Chen G, Ma L, Zhang SM, Jiang AP, Zheng XD, Gao XH, Li P, Tu CX, Yin XY, Han XP, Ren YQ, Song SP, Lu ZY, Zhang XL, Cui Y, Chang J, Gao M, Luo XY, Wang PG, Dai X, Su W, Li H, Shen CP, Liu SX, Feng XB, Yang CJ, Lin GS, Wang ZX, Huang JQ, Fan X, Wang Y, Bao YX, Yang S, Liu JJ, Franke A, Weidinger S, Yao ZR, Zhang XJ: Genome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population . Nat Genet 2011, 43 (7):690-4. Liao D: Emerging roles of the EBF family of transcription factors in tumor suppression . Mol Cancer Res 2009, 7 (12):1893-901. Jargosch M, Kröger S, Gralinska E, Klotz U, Fang Z, Chen W, Leser U, Selbig J, Groth D, Baumgrass R: Data integration for identification of important transcription factors of STAT6-mediated cell fate decisions . Genet Mol Res 2016, 15 (2). Moniz S, Biddlestone J, Rocha S: Grow2: the HIF system, energy homeostasis and the cell cycle . Histol Histopathol 2014, 29 (5):589-600. van der Bliek AM, Sedensky MM, Morgan PG: Cell Biology of the Mitochondrion . Genetics 2017, 207 (3):843-871. Ganesan V, Sivanesan D, Yoon S: Correlation between the Structure and Catalytic Activity of [Cp*Rh(Substituted Bipyridine)] Complexes for NADH Regeneration . Inorg Chem 2017, 56 (3):1366-1374. Li T, Liu J, Smith WW: Synphilin-1 binds ATP and regulates intracellular energy status . PLoS One 2014, 9 (12):e115233. Goody RS: The significance of the free energy of hydrolysis of GTP for signal-transducing and regulatory GTPases . Biophys Chem 2003, 100 (1-3):535-44. Nunes-Nesi A, Araújo WL, Obata T, Fernie AR: Regulation of the mitochondrial tricarboxylic acid cycle . Curr Opin Plant Biol 2013, 16 (3):335-43. Bal NC, Singh S, Reis FCG, Maurya SK, Pani S, Rowland LA, Periasamy M: Both brown adipose tissue and skeletal muscle thermogenesis processes are activated during mild to severe cold adaptation in mice . J Biol Chem 2017, 292 (40):16616-16625. Xiaogang Cui, Yali Hou, Shaohua Yang, Yan Xie, Shengli Zhang, Yuan Zhang, Qin Zhang, Xuemei Lu, George E Liu, Dongxiao Sun: Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing . BMC Genomics . 2014, 15 :226. Suárez-Vega A, Gutiérrez-Gil B, Klopp C, Tosser-Klopp G, Arranz JJ: Comprehensive RNA-Seq profiling to evaluate lactating sheep mammary gland transcriptome . Sci Data . 2016, 3 :160051. Mobuchon L, Marthey S, Boussaha M, Le Guillou S, Leroux C, Le Provost F: Annotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches . BMC Genomics . 2015, 16 (1):285. Yang Zhou, Erin E Connor, Derek M Bickhart, Congjun Li, Ransom L Baldwin, Steven G Schroeder, Benjamin D Rosen, Liguo Yang, Curtis P Van Tassell, George E Liu: Comparative whole genome DNA methylation profiling of cattle sperm and somatic tissues reveals striking hypomethylated patterns in sperm . Gigascience. 2018, 7 (5):giy039. Zhixiong Li , Mingfeng Jiang : Metabolomic profiles in yak mammary gland tissue during the lactation cycle . PLoS One 2019, 14 (7):e0219220. Qiumei Ji, Jinwei Xin, Zhixin Chai, Chengfu Zhang, Dawa Yangla, Luo Sang, Qiang Zhang, Pingcuo Zhandui, Minsheng Peng, Yong Zhu, Hanwen Cao, Hui Wang, Jianlin Han, Jincheng Zhong: A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak . Molecular Ecology Resources , DOI:10.1111/1755-0998.13236. In press. Xi Y, Li W: BSMAP: whole genome bisulfite sequence MAPping program . BMC Bioinformatics 2009, 10 (1):232-232. Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon GC, Tontifilippini J, Nery JR, Lee LK, Ye Z, Ngo Q: Human DNA methylomes at base resolution show widespread epigenomic differences . Nature 2009, 462 (7271):315-322. Juhling F, Kretzmer H, Bernhart SH, Otto C, Stadler PF, Hoffmann S: metilene: Fast and sensitive calling of differentially methylated regions from bisulfite sequencing data . Genome Research 2016, 26 (2):256-262. Wang, J., Duncan, D., Shi, Z., Zhang, B: WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update. Nucleic Acids Res 2013 , 41: W77-83. Pertea M, Pertea G, Antonescu C, Chang TC, Mendell JT, Salzberg SL: StringTie enables improved reconstruction of a transcriptome from RNA-seq reads . Nature Biotechnology 2015, 33 (3):290-295. Robinson MD, Mccarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data . Bioinformatics 2010, 26(1):139-140. Supplementary Files ARRIVEchecklist.docx TableS1.xlsx TableS2.xlsx TableS3.xlsx TableS4.xlsx TableS5.xlsx TableS6.xlsx Cite Share Download PDF Status: Published Journal Publication published 20 Oct, 2020 Read the published version in BMC Genomics → Version 4 posted Editorial decision: Minor revision 28 Sep, 2020 Editor assigned by journal 24 Sep, 2020 Submission checks completed at journal 23 Sep, 2020 Editor invited by journal 23 Sep, 2020 You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-20775","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":2806478,"identity":"353e54cb-3c54-4724-a37a-8f35aadf811a","order_by":0,"name":"Jinwei Xin","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jinwei","middleName":"","lastName":"Xin","suffix":""},{"id":2806479,"identity":"bb54add4-18b3-4ee9-9456-d08863dc1f89","order_by":1,"name":"Zhixin Chai","email":"","orcid":"","institution":"Southwest Minzu University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Zhixin","middleName":"","lastName":"Chai","suffix":""},{"id":2806480,"identity":"0ab75da0-d13a-452c-ae64-73d5dede9c59","order_by":2,"name":"Chengfu Zhang","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Chengfu","middleName":"","lastName":"Zhang","suffix":""},{"id":2806481,"identity":"c0033185-7afc-48b2-82d8-1ad620a1d190","order_by":3,"name":"Qiang Zhang","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Qiang","middleName":"","lastName":"Zhang","suffix":""},{"id":2806482,"identity":"a480f564-a393-4dec-9645-1bb57a1f2f41","order_by":4,"name":"Yong Zhu","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yong","middleName":"","lastName":"Zhu","suffix":""},{"id":2806483,"identity":"10d02837-2625-4fef-bc84-8bcf3ebffe55","order_by":5,"name":"Hanwen Cao","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Hanwen","middleName":"","lastName":"Cao","suffix":""},{"id":2806484,"identity":"e4df57b9-8bed-4d0e-9bb8-a7d79fee0ca1","order_by":6,"name":"Yangji Cidan","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yangji","middleName":"","lastName":"Cidan","suffix":""},{"id":2806485,"identity":"32db82af-357a-437e-b879-0a02b9320234","order_by":7,"name":"Xiaoying Chen","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xiaoying","middleName":"","lastName":"Chen","suffix":""},{"id":2806486,"identity":"a5c6a78b-0262-4553-b348-dbb3a98a6ad3","order_by":8,"name":"Hui Jiang","email":"","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Hui","middleName":"","lastName":"Jiang","suffix":""},{"id":2806487,"identity":"2af07741-1851-4b07-b7a9-02823f958d52","order_by":9,"name":"Jincheng Zhong","email":"","orcid":"","institution":"Southwest minzu university","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jincheng","middleName":"","lastName":"Zhong","suffix":""},{"id":2806488,"identity":"32a7067f-ee3d-4914-841d-9cae061e7226","order_by":10,"name":"Qiumei Ji","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAyUlEQVRIiWNgGAWjYBACxmYGhgNAmodN/vGBAx9+kKKFnyEt8eDMHlKsk2zIMT7MwUaESuZ25o2HeXNsZAwOnPlwmIGHQZ5f7AAhh7EVHObdlsZjcLB3w+ECCwbDmbMTCGnhMQBqOQwiNxyewcOQYHCbOC3/eQyO8Tw4zMNGvJYDPJI9PAzEamErODh3WzIPvwSbATCQJQj7xbD/8OYPb7fZ2bNJMD/+8OGHjTy/NCEtDQwGyHwJ/MpBQJ4BVcsoGAWjYBSMAkwAALx/RJBk78/OAAAAAElFTkSuQmCC","orcid":"","institution":"State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Qiumei","middleName":"","lastName":"Ji","suffix":""}],"badges":[],"createdAt":"2020-04-01 14:44:40","currentVersionCode":4,"declarations":"","doi":"10.21203/rs.3.rs-20775/v4","doiUrl":"https://doi.org/10.21203/rs.3.rs-20775/v4","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12864-020-07151-3","type":"published","date":"2020-10-20T12:00:00+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":2689199,"identity":"7637adca-6c38-43f4-bd0a-1184769e0e22","added_by":"auto","created_at":"2020-09-29 21:32:37","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":109882,"visible":true,"origin":"","legend":"Global DNA methylation and gene expression among samples.\nPearson’s correlation analysis based on the methylation of CpG sites (a) and gene expression (b) among samples. Boxplot of Pearson’s correlation coefficients between replicates, ages, or tissues for methylation (c) and gene expression (d). M6, M30, M54 and M90 represent different 6, 30, 54 and 90 months old, respectively. B, L, and M represent breast, lung, and biceps brachii muscle, respectively. 1, 2, and 3 represent different replicates.\n","description":"","filename":"fig1.png","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/fig1.png"},{"id":2689200,"identity":"f1b697c0-fe46-4350-ba8f-d31ee00afd03","added_by":"auto","created_at":"2020-09-29 21:32:38","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":189604,"visible":true,"origin":"","legend":"Overview of age-associated DMRs. (a) Basic statistics for A-DMRs within each tissue. Overlap of A-DMRs associated with the 6-month group (b), 90-month group (c), and 30- and 54- months groups (d) in the muscle, breast, and lung respectively.","description":"","filename":"2.PNG","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/2.PNG"},{"id":2689203,"identity":"c243f0e8-e4ad-4fd7-8b79-fd269ad2d6b1","added_by":"auto","created_at":"2020-09-29 21:32:39","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":56785,"visible":true,"origin":"","legend":"Modules of consensus networks and correlation with traits. Consensus networks from the age (a) or tissue (b) curves. Gene expression similarity was determined using a pair-wise weighted correlation metric and clustered according to a topological overlap metric into modules; assigned modules are colored on the bottom, and gray genes were not assigned to any module. Consensus network modules for age (c) and tissue (d) correlated with traits using the eigenmodule (the first principal component of the module). The correlation coefficients and the p-value in parenthesis are provided underneath; color-coding refers to the correlation coefficient (legend at right).","description":"","filename":"fig3.png","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/fig3.png"},{"id":2689205,"identity":"c4f666ae-0c59-48cd-8027-edb4b95693f1","added_by":"auto","created_at":"2020-09-29 21:32:39","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":110010,"visible":true,"origin":"","legend":"“Hub” genes and potential target genes of ZGPAT in the age network. (a) Expression level of 27 “hub” genes. (b) Enrichment analysis of ZGPAT’s potential target genes. (c) Expression level of 7 genes in “protein processing in endoplasmic reticulum”, which was enriched from potential target genes of ZGPAT.","description":"","filename":"fig4.png","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/fig4.png"},{"id":2689207,"identity":"4091758a-dfee-49ab-9ebf-dceca05376d5","added_by":"auto","created_at":"2020-09-29 21:32:40","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":282769,"visible":true,"origin":"","legend":"Modules and “hub” genes in the tissue network. \n(a) WGCNA modules of the tissue-related genes, (b) correlations between modules showed by the eigenmodule adjacency heatmap, (c) expression level of “hub” genes in the tissue network, (d) enrichment analysis of potential target genes of HIF1A, and the number of enriched genes and enrichment fold are indicated on the right.\n","description":"","filename":"fig5.png","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/fig5.png"},{"id":13597652,"identity":"4ae79507-087b-4909-b437-b46dfbeffd17","added_by":"auto","created_at":"2021-09-17 05:33:06","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1728412,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/3180a08b-9c2a-4c31-9cef-32711dee58ea.pdf"},{"id":2689198,"identity":"da9cbc3b-cc81-44df-a1cb-0266c371c4e1","added_by":"auto","created_at":"2020-09-29 21:32:37","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":26153,"visible":true,"origin":"","legend":"","description":"","filename":"ARRIVEchecklist.docx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/ARRIVEchecklist.docx"},{"id":2689201,"identity":"19544a4a-9879-4613-bec4-d6ec7fbd55ba","added_by":"auto","created_at":"2020-09-29 21:32:38","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":8930391,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS1.xlsx"},{"id":2689202,"identity":"7afbd393-c931-4375-b659-c4d4ee1b6174","added_by":"auto","created_at":"2020-09-29 21:32:39","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":31358266,"visible":true,"origin":"","legend":"","description":"","filename":"TableS2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS2.xlsx"},{"id":2689204,"identity":"1902e2e7-5b43-40b3-848f-491180be5df4","added_by":"auto","created_at":"2020-09-29 21:32:39","extension":"xlsx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":13472105,"visible":true,"origin":"","legend":"","description":"","filename":"TableS3.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS3.xlsx"},{"id":2689206,"identity":"f5fc1fac-e6ca-4fc6-9a8d-3bfb088062d4","added_by":"auto","created_at":"2020-09-29 21:32:40","extension":"xlsx","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":22699694,"visible":true,"origin":"","legend":"","description":"","filename":"TableS4.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS4.xlsx"},{"id":2689208,"identity":"fe36af57-6232-4431-986d-868308b8337c","added_by":"auto","created_at":"2020-09-29 21:32:40","extension":"xlsx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":118319,"visible":true,"origin":"","legend":"","description":"","filename":"TableS5.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS5.xlsx"},{"id":2689209,"identity":"6337dbee-c0b3-41a2-86c6-55d2f888879d","added_by":"auto","created_at":"2020-09-29 21:32:40","extension":"xlsx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":3768055,"visible":true,"origin":"","legend":"","description":"","filename":"TableS6.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-20775/v4/TableS6.xlsx"}],"financialInterests":"","formattedTitle":"Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production","fulltext":[{"header":"Background","content":"\u003cp\u003eDomestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), in the Himalayas, and in the connecting Central Asian highlands. They provide milk, meat, hides, fiber, fuel, and transportation [1, 2]. Milk is an important source of high-quality protein. It contains high quantities of essential amino acids, such as lysine, which is commonly deficient in many human diets [3]. Milk proteins also impact immunomodulatory processes and gastrointestinal activities [4]. Protein content and composition influence the technological properties of milk and are\u0026nbsp;therefore important for the dairy industry, particularly in Europe, where the majority of the milk produced is used to produce cheese. In recent decades, there have been extraordinary advances in our knowledge of the physiology and biochemistry of the lactating mammary gland. Previous research indicates that milk protein synthesis in the mammary gland depends on hormonal and developmental cues that modulate the transcriptional and translational regulation of genes through the activity of specific transcription factors, non-coding RNAs, and alterations of the chromatin structure in mammary epithelial cells [5]. The interplay between these factors may influence milk protein synthesis, which is crucial during the entire lactation process in high-producing dairy cattle.\u0026nbsp;Despite such advancements in research, little is currently known about the physiological and cellular regulation required for milk protein synthesis and secretion in yak. We hypothesized that the genes responsible for milk production were regulated by DNA methylation and that distinct sub-modules of correlated expression variation could be identified. In this study, we performed genome-wide DNA methylome and transcriptome analyses of yak lung, breast, and biceps brachii muscle tissues at four different stages of development to identify the regulatory networks associated with milk protein synthesis, metabolism, and secretion.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eGlobal DNA methylation and gene expression in the breast, lungs, and biceps brachii muscle at different ages \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe generated the methylomes and transcriptomes of lung, breast, and biceps brachii muscle tissues from 12 female Riwoqe yaks at four different stages (n=3/stage) of development: 6 months old (MO) (young), 30 MO (pre-mature), 54 MO (mature) and 90 MO (post mature). Among these, only the post-mature yaks (90 months) were lactation, with ~3.16 kg/day milk yield. After performing sequence quality control and filtering, we obtained a single-base resolution methylome covering 85.6% (27,471,373/32,092,725) of CpG sites across the genome with an average depth of 22.5\u0026times;. We first calculated pairwise Pearson's correlations of CpG sites with at least 10\u0026times; coverage depth across all samples, which were well clustered by tissue types (Figure 1a). The correlation of CpG methylation levels for biological replicates was strong (median Pearson's r = 0.74), the correlation of CpG methylation levels between ages (median Pearson's r = 0.72) was relatively weaker and the correlation between tissues was weakest (median Pearson's r = 0.66) (Figure 1c). We aligned the \u0026nbsp;transcriptome sequencing data for all samples to our newly assembled yak genome reference, and we subsequently obtained the transcripts. In total, we obtained 2,0504 transcripts, that were then annotated to the Gene\u0026nbsp;Ontology (GO) [6], InterPro [7], Kyoto Encyclopedia of Genes and Genomes (KEGG) [8], Swiss-Prot [9], and TrEMBL [10] databases (Table S1). We also calculated pairwise Pearson's correlations of all transcripts and obtained similar results to those of DNA methylation (Figure 1b). Biological replicates showed the highest correlation coefficients, while different tissues showed the lowest correlation coefficients (Figure 1d).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDifferentially methylated regions among the age groups \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe determined differentially methylated regions (DMRs) across age groups within the breast, lung, and biceps brachii muscle tissues (Table S2-4). Within the lung and biceps brachii muscle tissues, age groups did not differ in age-related DMRs (A-DMRs), but post-mature breast tissue had considerably more DMRs (155,957) than the three younger age groups (Figure 2a). We then investigated the correlations between DMRs and their corresponding differentially expressed genes. The ratio of negatively to positively correlated gene pairs was 1.02 for promoters with DMRs and 0.95 for gene bodies with DMRs. Not every methylation was correlated with the expression of its associated gene, due to the gene regulation complexity [11].\u003c/p\u003e\n\u003cp\u003eAt ~90 months, ~120 days after giving birth for the third time, yaks are in the lactation period (milk yield of ~3.16kg/day)(Li \u0026amp; Jiang, 2019), so it is possible that the observed methylation partially controls yak lactation. Since promoter methylation decreases gene expression [11], we selected 375 hypomethylated promoter genes (highly expressed) along with 207 hypermethylated promoter genes (lowly expressed) from post-mature yak breast tissues. The hypomethylated (highly expressed) genes were only enriched in\u0026ldquo;protein processing in endoplasmic reticulum (ER)\u0026rdquo; (9 genes, 2.964-fold enrichment, p = 0.0049). Specifically, the genes are involved in vesicle trafficking (SEC23B), oligosaccharide linking (MOGS, RPN2), folding and assembly (HSPA5), transportation (LMAN2, SEL1L), and ubiquitination and degradation (UBE2J1, UBE2J2, DERL2) [12, 13]. These genes were significantly upregulated at 90 months in breast tissues, but not in lung and biceps brachii muscle tissues (Figure 2e). Based on this data, methylation might help regulate milk production by influencing protein processing in the endoplasmic reticulum during the lactation period.\u003c/p\u003e\n\u003cp\u003eWe also examined A-DMRs that overlapped across age groups. Young and post-mature tissues rarely shared A-DMRs when comparing the lung and muscle tissues (for young tissues, muscle: 586 A-DMRs, breast: 2,249 A-DMRs, lung: 496 A-DMRs; for post-mature tissues, muscle: 470 A-DMRs, breast: 12,050 A-DMRs, lung: 772 A-DMRs) (Figure 2b, c). Pre-mature and mature stages also rarely shared A-DMRs across muscle and lung tissues (Figure 2d), suggesting that methylation patterns were already established at the young stage and that no extensively divergent epigenetic difference occurred across different age groups under natural high-altitude conditions.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsensus network analysis for tissues and age groups\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe first performed a multi-way ANOVA test for each gene across all samples (n=36) to test the null hypothesis that the gene expression level did not differ among age groups and tissues. At the threshold\u0026nbsp;for significance (p\u0026lt;0.05), 417 age-related and 8,560 tissue-related genes were selected for further weighted gene correlation network analysis (WGCNA), which uses network topography to group genes into modules based on correlations [14]. Next, we conducted WGCNA for tissue- and age-related gene expression respectively to identify a \u0026ldquo;consensus network\u0026rdquo;\u0026ndash;a common pattern of genes that are correlated in all conditions. We performed a consensus network, module statistic, and eigengene network analyses to identify modules, to assess relationships between modules and traits, and to study the relationships between co-expression modules [15]. The consensus networks identified for tissues and age groups had clearly delineated modules (Figure 3a, 3b), and the modules identified were significantly correlated with tissues and age groups (Figure 3c, 3d).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAge network analysis indicates that ZGPAT might regulate milk production\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWithin the age-related gene network, the largest module (\u0026ldquo;turquoise\u0026rdquo;, n=356) had negative correlation for breast tissue (r=-0.57, p=3e-04) and age (r= 0.37, p=0.03) and a positive correlation for biceps brachii muscle tissues (Table S5). The breast tissue had a stronger signal than that of age and may have overwhelmed the signal from age. Genes in this module were enriched in the GO categories of \u0026ldquo;protein polyubiquitination,\u0026rdquo; \u0026ldquo;RNA polymerase II core promoter proximal region sequence-specific DNA binding,\u0026rdquo; \u0026ldquo;ATP binding,\u0026rdquo; \u0026ldquo;transcription, DNA-templated,\u0026rdquo; and \u0026ldquo;negative regulation of transcription from RNA polymerase II promoter\u0026rdquo; (p-values of 0.000344, 0.000822, 0.000991, 0.00139, and 0.00244, respectively). The \u0026ldquo;blue\u0026rdquo; (n=48) and \u0026ldquo;grey\u0026rdquo; (n=13) modules showed only a negative correlation with age (r= -0.58, p= 2e-04; r= -0.76, p= 6e-08, respectively) and exhibited no enrichment of GO categories for genes. After applying the threshold of the absolute\u0026nbsp;value of gene significance for age (|GS| \u0026gt;0.5) and module membership measures (|MM| \u0026gt;0.6) in each module, we defined 20 and 7 \u0026ldquo;hub\u0026rdquo; genes in the \u0026ldquo;turquoise\u0026rdquo; and \u0026ldquo;blue\u0026rdquo; modules (Table 1). The gene expression of the \u0026ldquo;hub\u0026rdquo; genes was well clustered by modules, which was consistent with the negative correlation with age (\u0026ldquo;turquoise\u0026rdquo; r=0.37, \u0026ldquo;blue\u0026rdquo; r=-0.58). The upregulated expression level of the \u0026ldquo;hubs\u0026rdquo; in breast tissue at 90 months old indicated that the \u0026ldquo;turquoise\u0026rdquo; module had a stronger correlation with breast tissue than with age (Figure 4a).\u003c/p\u003e\n\u003cp\u003eTable 1. List of \u0026ldquo;hub\u0026rdquo; genes in the consensus network for age\u003c/p\u003e\n\u003ctable border=\"1\" width=\"643\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd width=\"99\"\u003e\n\u003cp\u003eYak ID\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"87\"\u003e\n\u003cp\u003eGene symbol\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"91\"\u003e\n\u003cp\u003eModule color\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"107\"\u003e\n\u003cp\u003eGene significance\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003ep-value\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"116\"\u003e\n\u003cp\u003eModule membership\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003ep-value\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB009868\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eAP3D1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.513965075\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001344116\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.868945999\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e6.37E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB004083\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eTMEM30A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.505722627\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001652657\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.848054843\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e6.65E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB000726\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eDDRGK1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.53596122\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000754414\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.842070133\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.22E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB019011\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eOTUD3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.543916749\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000606215\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.828680521\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e4.37E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB000091\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eCNPPD1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.51195228\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001414359\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.814537409\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.50E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB006815\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eTOR1AIP1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.544705745\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000593032\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.809789221\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.21E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB007137\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eMGAT4A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.545969144\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000572452\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.798170686\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e5.51E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB007385\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eHECA\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.603273223\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e9.84E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.770629907\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e3.85E-08\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB019443\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSYAP1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.500377628\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001884511\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.733120646\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e3.68E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB009825\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003ePHAX\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.561938984\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.00036184\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.70497442\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.59E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB008032\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eUBE2G2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.52385538\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001041698\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.70208134\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.83E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB012517\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSYF2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.568328068\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000299194\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.699367268\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.08E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB001610\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eFURIN\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.567906296\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000303008\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.697046314\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.32E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB000679\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eZGPAT\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.504624296\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001698141\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.683958391\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e4.25E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB007049\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSKP2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.527608908\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000943732\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.669790419\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e7.91E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB000224\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eC2orf6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.531191165\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000857925\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.652871802\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.59E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB007420\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eTAB2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.505905132\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001645204\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.652844304\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.59E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB018569\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eORAOV1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.51035386\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001472421\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.637043276\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.95E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB010550\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eCCPG1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.644779605\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.19E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.612874899\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e7.08E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB012521\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eTMEM57\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.57317319\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000258349\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.60967979\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e7.91E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB013324\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eMCM3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.583585649\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000186982\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.841505398\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.29E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB010064\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSPC24\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.509962181\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001486965\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.744497748\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.93E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB003102\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eC17orf49\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.526803776\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000964031\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.741805787\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.26E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB015902\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSERPINH1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.607769763\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e8.44E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.721017606\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e7.04E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB016582\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eSRPX2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.51838468\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001200558\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.698151166\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e2.20E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB012996\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eUCK2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.588274454\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.000161067\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.677420455\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e5.68E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd\u003e\n\u003cp\u003eBmuPB015372\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eUMPS\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003eblue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e-0.503577413\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.001742514\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e0.648111658\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd\u003e\n\u003cp\u003e1.92E-05\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe used the AnimalTFDB 3.0 database [16] to examine transcription factors in these 27 \u0026ldquo;hubs\u0026rdquo; and found that ZGPAT encodes a transcription regulator protein and was significantly upregulated in breast tissue at 90 months of age (Figure 4a). Previous study reported that this protein specifically binds the 5'-GGAG[GA]A[GA]A-3' consensus sequence and represses transcription by recruiting the chromatin multi-complex NuRD to target promoters [17]. ZGPAT was highly expressed in breast tissue at 90 months, and it potentially regulated the transcription of 280 genes (weight \u0026gt;0.15) in the network from the \u0026ldquo;turquoise\u0026rdquo; module. In order to identify the most important cellular activities controlled by this TF regulatory network, we analyzed over-represented GO biological process and molecular function terms, as well as KEGG pathways. These potential target genes were enriched in the GO categories of \u0026ldquo;protein binding,\u0026rdquo; \u0026ldquo;ATP binding,\u0026rdquo; and \u0026ldquo;zinc ion binding,\u0026rdquo; among others, and the KEGG categories of \u0026ldquo;aminoacyl-tRNA biosynthesis,\u0026rdquo; \u0026ldquo;autophagy animal,\u0026rdquo; and \u0026ldquo;protein processing in endoplasmic reticulum\u0026rdquo; (Figure 4b). These enriched GO terms and KEGG pathways likely help regulate protein\u0026nbsp;synthesis, processing, and secretion in breast tissue. For example, 6 of 7 genes from the \u0026ldquo;protein processing in endoplasmic reticulum\u0026rdquo; category were also upregulated at 90 months of age in breast tissue (Figure 4c) and involved in multiple processes in the endoplasmic reticulum, including vesicle trafficking (SEC24C), folding and assembly (SELENOS), transportation (BCAP31), and ubiquitination and degradation (BAG1, UBE2G2, and MARCH6) [12, 13]. Only DNAJC10 was downregulated at 90 months of age in breast tissue. This gene encodes an endoplasmic reticulum co-chaperone that is part of the endoplasmic reticulum-associated degradation complex involved in recognizing and degrading misfolded proteins [13].\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTissue network analysis indicates regulative role of HIF1A in lung\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWithin the tissue-related gene module network, four modules showed a positive correlation and two showed a negative correlation with the lung; all significant module-trait relationships were negative in the muscle but positive in the breast (Figure 3d). Moreover, 99.54% of the total 8,560 tissue-related genes were related to the top 4 modules (\u0026ldquo;turquoise,\u0026rdquo; n=3833; \u0026ldquo;blue,\u0026rdquo; n=2795; \u0026ldquo;brown,\u0026rdquo; n=1052; \u0026ldquo;yellow,\u0026rdquo; n=339) (Figure 5a, Table S6), and these modules were also highly correlated with other modules; for example, \u0026ldquo;brown,\u0026rdquo; \u0026ldquo;yellow,\u0026rdquo; and \u0026ldquo;black\u0026rdquo; showed a high eigengene adjacency with each other (Figure 5b).\u003c/p\u003e\n\u003cp\u003eWe applied the more stringent threshold absolute\u0026nbsp;value of gene significance for the age and module membership measures in the top four modules to identify \u0026ldquo;hub\u0026rdquo; genes in the \u0026ldquo;turquoise,\u0026rdquo; \u0026ldquo;blue,\u0026rdquo; \u0026ldquo;brown,\u0026rdquo; and \u0026ldquo;yellow\u0026rdquo; modules. With the threshold values of |GS|\u0026gt;0.7 and |MM| \u0026gt;0.8, 34 \u0026ldquo;hub\u0026rdquo; genes were identified in the \u0026ldquo;turquoise\u0026rdquo; module. Twenty-four hubs were then filtered from the gene significance of module-lung relationships, and 10 were filtered from the gene significance of module-breast relationships. These were further divided into 3 clusters by hierarchical clustering, which showed high expression levels in the breast (cluster 1), lung (cluster 2), and biceps brachii muscle (cluster 3) tissues, with distinct clustering patterns by tissue (Figure 5c).\u003c/p\u003e\n\u003cp\u003eTable 2. List of \u0026ldquo;hub\u0026rdquo; genes in the consensus network for tissue\u003c/p\u003e\n\u003ctable border=\"1\" width=\"684\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eYak ID\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eGene symbol\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eModule color\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003eTissue\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003eGene significance\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003ep-value\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003eModule membership\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003ep-value\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB014336\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eEEF1G\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.73907039\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e2.64E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.826748126\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e5.20E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB017352\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003ePMS1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.71599797\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e9.13E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.850552958\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e5.12E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB000878\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eMTPAP\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.715583465\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e9.33E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.912877613\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e8.73E-15\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB018762\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eCXHXorf58\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.709433358\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.27E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.82262518\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e7.50E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB014450\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eRWDD4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.708526806\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.33E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.923713618\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e9.94E-16\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB011005\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eHUS1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.707137304\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.43E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.875052199\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.97E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB005540\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eMTUS1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.704482832\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.62E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.871626441\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e4.58E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB011453\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eEBF3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.703373458\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.71E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.870717052\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e5.13E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB010608\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eLAMB3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.70783756\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.38E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.85067459\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e5.05E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB004299\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eC3H3orf58\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.708872549\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.31E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.88522675\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e7.61E-13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB020871\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eG6PD\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.708930411\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.30E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.90727231\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.41E-14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB007438\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eZCCHC6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.709740967\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.25E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.84738754\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e7.13E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB007592\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eCTDSPL\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.711634321\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.14E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.8813581\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.30E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB004894\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eHIF1A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.713237417\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.05E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.87302984\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e3.84E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB020508\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eFAM122B\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.720139618\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e7.37E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.82428698\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e6.48E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB018102\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eCCDC82\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.720592854\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e7.20E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.90665968\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.68E-14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB014539\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eCNTRL\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.722130497\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e6.64E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.87149974\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e4.65E-12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB003882\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eVAV3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.725601285\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e5.53E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.89496378\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.82E-13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB020153\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eEPS8L1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.727495147\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e4.99E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.82286844\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e7.34E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB010308\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003ePGM2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.746811608\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.69E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.83795918\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.83E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB004008\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eGDAP2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.754806874\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.05E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.88661303\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e6.26E-13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB012568\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eWASF2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.757428944\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e8.92E-08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.83874427\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.69E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB011966\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eACTG1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.772114072\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e3.49E-08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.84373776\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.03E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB014173\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eSTAT6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003elung\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.80896084\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e2.37E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.84981355\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e5.53E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB015106\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eMAN2A1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.705124693\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.57E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.84782348\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e6.81E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB008614\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eVPS26A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.731553581\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e4.01E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.88669366\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e6.18E-13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB018496\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eFAM92A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.718223061\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e8.14E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.85598454\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.85E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB005078\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eRASA1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.705902582\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.51E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.85869193\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.11E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB011476\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eFAM53B\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.728187435\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e4.81E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.848680696\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e6.23E-11\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB008348\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eEPDR1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.710771064\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e1.19E-06\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.832481139\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e3.08E-10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB003453\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eTROVE2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.72654448\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e5.26E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.80669368\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.84E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB021200\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eVWA7\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.722814498\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e6.41E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.814033145\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e1.56E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB010528\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eRNF111\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e0.737234077\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e2.92E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e-0.80146106\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e4.28E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"97\"\u003e\n\u003cp\u003eBmuPB015811\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"88\"\u003e\n\u003cp\u003eFRMD3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eturquoise\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"60\"\u003e\n\u003cp\u003ebreast\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"95\"\u003e\n\u003cp\u003e-0.739223421\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e2.61E-07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"84\"\u003e\n\u003cp\u003e0.809147869\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"76\"\u003e\n\u003cp\u003e2.33E-09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAccording to the AmalTFDB 3.0 database [16], EBF3, HIF1A, and STAT6 were annotated as transcription factors. EBF3 encodes a member of the early B-cell factor (EBF) family of DNA binding transcription factors. EBF proteins are involved in B-cell differentiation, bone development, and neurogenesis and they may also function as tumor suppressors [18]. STAT6\u0026nbsp;is a member of the STAT family of transcription factors. In response to cytokines and growth factors, STAT family members are phosphorylated by the receptor-associated kinases. They form homo- or heterodimers that translocate to the cell nucleus where they act as transcription activators [19]. HIF1A, hypoxia-inducible factor-1, functions as a master regulator of the cellular and systemic homeostatic response to hypoxia by activating the transcription of genes, involved in energy metabolism, angiogenesis, and apoptosis, as well as other genes whose protein products increase oxygen delivery or facilitate metabolic adaptation to hypoxia [20]. HIF1A was upregulated in the lung to potentially control the transcription of 2008 genes (weight \u0026gt;0.15), which were enriched in multiple GO biological process categories, molecular function categories and KEGG pathways. Most of the enriched GO terms and KEGG pathways were related to energy metabolism. Specifically, enriched GO terms included \u0026ldquo;mitochondrial respiratory chain complex I assembly,\u0026rdquo; \u0026ldquo;NADH dehydrogenase (ubiquinone) activity,\u0026rdquo; \u0026ldquo;ATP binding,\u0026rdquo; \u0026ldquo;mitochondrial translation,\u0026rdquo; \u0026ldquo;tricarboxylic acid cycle,\u0026rdquo; and \u0026ldquo;GTP binding\u0026rdquo; while enriched KEGG pathways included \u0026ldquo;thermogenesis,\u0026rdquo; \u0026ldquo;carbon metabolism,\u0026rdquo; and \u0026ldquo;citrate cycle TCA cycle\u0026rdquo; (Figure 5d).\u0026nbsp;Mitochondria\u0026nbsp;function as the primary energy producers of the cell and serves as the hub for a variety of immune pathways as the center of biosynthesis, oxidative stress response, and cellular signaling [21].\u0026nbsp;NADH dehydrogenase is a core subunit of the mitochondrial membrane respiratory chain and is believed to contribute to the minimal assembly required for catalysis [22].\u0026nbsp;\u0026nbsp;The protein which binds ATP or GTP, carries three phosphate groups esterified to a sugar moiety and provides energy and phosphate sources for the cell [23, 24]. The tricarboxylic acid cycle, a series of metabolic reactions in aerobic cellular respiration, occurs in the mitochondria of animals and plants. During this cycle, acetyl-CoA, which is formed from pyruvate produced during glycolysis, is completely oxidized to CO\u003csub\u003e2\u003c/sub\u003e via interconversion of various carboxylic acids. This results in the reduction of NAD and FAD to NADH and FADH2, respectively, and indirectly in the synthesis of ATP by oxidative phosphorylation [25]. The \u0026ldquo;thermogenesis\u0026rdquo; pathway is essential for warm-blooded animals, because it ensures normal cellular and physiological functions under challenging environmental conditions [26].\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003ePrevious studies described transcription profiling of the mammary gland in livestock (including cattle [27], sheep [28], and goat [29]) and DNA methylation profiling of the mammary gland in cattle [30]. These studies showed temporal and spatial specificity in the methylome and transcriptome profiles of the mammary gland in different species, but they only reported the\u0026nbsp;differential gene expression profile of the mammary gland; and the regulatory network is still unknown. In this study, we for the first time generated the methylomes and transcriptomes of lung, breast, and biceps brachii muscle tissues from yak at four different stages of development (6, 30, 54, and 90 months; young, pre-mature, mature, and post-mature, respectively). We found that breast tissue at 90 months showed considerably differential methylation levels compared with other month groups, but lung and biceps brachii muscle tissues did not. Enrichment analysis for upregulated genes with hypomethylated DMRs from breast tissues of 90-month-old yaks showed that DNA methylation might regulate the activation of the \u0026ldquo;protein processing in endoplasmic reticulum (ER)\u0026rdquo; pathway. Because only the 90-month-old yaks were in the lactation period, it appears that DNA methylation regulates milk production by influencing protein processing in the endoplasmic reticulum.\u003c/p\u003e\n\u003cp\u003eIn this study, hub genes were identified by WGCNA. The data show that the hub genes with the highest MM and GS in modules of interest should be candidates for further research. This study identified turquoise module genes associated with milk yield, and 20 genes were considered hub genes, showing the highest mRNA expression level in breast tissue at 90 months when yaks enter the lactation period. In these hub genes, ZGPAT was annotated as a transcription factor that potentially regulated the transcription of 280 genes in the \u0026ldquo;turquoise\u0026rdquo; module network, which was enriched in the KEGG categories of \u0026ldquo;aminoacyl-tRNA biosynthesis,\u0026rdquo; \u0026ldquo;autophagy animal,\u0026rdquo; and \u0026ldquo;protein processing in endoplasmic reticulum\u0026rdquo;. This result suggests that ZGPAT helps with regulating\u0026nbsp; protein\u0026nbsp;synthesis, processing, and secretion in breast tissue. Moreover, the 7 genes potentially regulated by ZGPAT in \u0026ldquo;protein processing in endoplasmic reticulum\u0026rdquo; were totally different from the aforementioned 9 genes regulated by hypomethylation, illustrating that DNA methylation and transcription factor possibly co-regulate milk production. In addition, the tissue network analysis confirms the importance of HIF1A in regulating energy metabolism, which is necessary for adaptations to low temperature and hypoxia in high altitudes.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eThe results of this comprehensive study provide a solid basis for understanding the roles of DNA methylation and transcriptional network underlying milk protein synthesis and high-altitude adaptation in yaks. This information advances our understanding of the regulatory network in the mammary gland at different developmental stages and may help inform breeding programs aimed at improving milk production.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eAnimals and samples \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn total, twelve female yaks (belonging to an indigenous yak breed that lives at altitudes of 3800-4000 meters above sea level in Riwoqe, Tibet, China) were collected between\u0026nbsp;June and December of 2016 from private farms.\u0026nbsp;They were grouped with 3 replicates into age categories of 6, 30, 54, and 90 months. At the time of slaughter, their mean live weights were 44.93 kg (6 months old), 153.06 kg (30 months old), 188.3 kg (54 months old), and 243.56 kg (90 months old). Only the 90 month-old yaks were lactating, producing ~3.16 kg/day milk yield (~120 days after giving birth for the third time). The 54-month-old yaks were in a dry period (~540 days after giving birth for the first time) [31]. The blood relation between the last three generations, which were housed simultaneously and fed the same diets, was unknown. \u0026nbsp;The yaks were not fed the night before they were slaughtered and were humanely sacrificed with the following procedures: (1) showered with clean water close to body temperature (35-38\u0026deg;C), (2) electrically stunned (120V dc, 12s) prior to exsanguination, (3) sacrificed while in the coma by bloodletting from carotid artery and jugular vein, (4) and dissected rapidly to obtain breast, lung, and biceps brachii muscle tissue samples, which were immediately frozen in liquid nitrogen, and stored at \u0026minus;80\u0026deg;C until RNA and DNA extraction.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWhole genome bisulfite sequencing \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) was used to isolate the high-quality DNA from each sample. According to the manufacturer\u0026rsquo;s instructions, 1 \u0026mu;g of genomic DNA was fragmented by sonication to a mean size of approximately 250 bp and subsequently used for whole genome bisulfite sequencing (WGBS) library construction using an Acegen Bisulfite-Seq Library Prep Kit (Acegen, Shenzhen, GD, China). Briefly, fragmented DNA was end-repaired, 5'-phosphorylated, 3'-dA-tailed, and then ligated to methylated adapters. The methylated adapter-ligated DNAs were purified using 1\u0026times; Agencourt AMPure XP magnetic beads (Beckman\u0026nbsp;Coulter, Brea,\u0026nbsp;CA, USA) and subjected to bisulfite conversion with a ZYMO EZ DNA Methylation-Gold Kit (Zymo\u0026nbsp;research, Irvine, CA, USA). The converted DNAs were then amplified using 25 \u0026mu;l HiFi HotStart U+ RM and 8-bp index primers with a final concentration of 1 \u0026mu;M each. The constructed WGBS libraries were then analyzed with an Agilent 2100 Bioanalyzer (Agilent Technologies, SantaClara, CA, USA), quantified with a Qubit fluorometer with Quant-iT dsDNA HS Assay Kit (Invitrogen,\u0026nbsp;Carlsbad, CA,USA), and finally sequenced on an Illumina Hiseq X ten sequencer (PE150 mode) (Illumina, San Diego, CA, USA)\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMethylation calculation and identification of DMRs \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLow quality reads that contained more than 5 \u0026lsquo;N\u0026rsquo;s or had a low-quality value for over 50% of the sequence (Phred score \u0026lt; 5) were filtered. The sequencing reads of the samples were aligned to the yak reference genome [32] using BSMAP (Version 2.74) [33]. The methylated CpG (mCG) sites were identified following a previously described algorithm [34]. The methylation levels for each sample were calculated using in-house Perl scripts. Differentially methylated regions (DMRs) were identified using metilene (Version 0.2-6) within a 500 bp sliding window at 250 bp steps with at least 10 CpGs covered by over 10\u0026times; sequence reads, applying the thresholds of differential methylation \u0026beta; \u0026gt;=15%, FDR for two-dimensional Kolmogorov-Smirnov-Test p-value \u0026lt;0.05. [35]. The enrichment analyses were conducted using WebGestalt (WEB-based Gene SeT AnaLysis Toolkit)\u0026nbsp;[36].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTotal RNA extraction, library preparation, and sequencing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe TRIzol reagent (Invitrogen, Carlsbad, CA, USA) was used to isolate the total RNA of each sample. The purity, concentration, and integrity of RNA were checked using the NanoPhotometer spectrophotometer (IMPLEN, Westlake Village, CA, USA), the Qubit RNA Assay Kit in Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), and the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 System (Agilent Technologies, SantaClara, CA, USA), respectively. For each sample, \u0026nbsp;3 \u0026micro;g high-quality RNA was used as input material for RNA-seq library preparation. First, ribosomal RNA was removed using the Epicentre Ribo-Zero rRNA Removal Kit (Epicentre, Madison, WI, USA). Next, the rRNA-depleted RNA was used to create sequencing libraries using the NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). Finally, the library products were purified using 1\u0026times; Agencourt AMPure XP magnetic beads (Beckman\u0026nbsp;Coulter, Brea,\u0026nbsp;CA, USA) and the Agilent Bioanalyzer 2100 System (Agilent Technologies, SantaClara, CA, USA) was employed to assess the library quality. Clustering of the index-coded samples was completed on a cBot Cluster Generation System using the TruSeq PE Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA), and then the libraries were sequenced on the Illumina HiSeq X Ten Platform to generate 150\u0026thinsp;bp paired-end reads.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eQuality analysis, transcriptome assembly, and abundance estimation\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eClean reads were obtained by removing reads containing the adapter or poly-N and by removing low quality reads (over 10% of the sequence with a quality value \u0026lt; 30) from the raw data using in-house Perl scripts. All downstream analyses were based on the good-quality clean reads. Paired-end clean reads were mapped to the yak reference genome [32]\u0026nbsp; with STAR (available at https://github.com/alexdobin/STAR/releases). The mapped reads of each sample were assembled using StringTie [37]. Next, all sample transcriptomes were merged to reconstruct a comprehensive transcriptome using Perl scripts. After the final transcriptome was generated, StringTie and edgeR were used to estimate the expression levels of all transcripts [38]. StringTie was used to assess the expression level of mRNAs by calculating fragments per kilobase of transcript per million fragments mapped (FPKM). Differentially expressed mRNAs were identified using the DESeq2 package, with the criteria of fold-change log2 \u0026gt; 1 or log2 \u0026lt; -1 and with the statistical significance set to FDR \u0026lt; 0.05.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWeighted gene correlation network analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA WGCNA network [14] was generated for both age-related and tissue-related genes. Consensus networks and module statistics followed the overall approach described by Langfelder et al. (2008). The network was derived based on a signed Spearman correlation using a b of 10 as a weight function. The topological overlap metric (TOM) [15] was derived from the resulting adjacency matrix and used to cluster the modules using the blockwiseModules function (blockwise Consensus Modules, for the consensus modules) and the dynamic tree cut algorithm [15] with a height of 0.25 and a deep split level of 2, a reassign threshold of 0.2, and a minimum module size of 30 (100 for the consensus network). Each eigenmodules\u0026ndash;the first principal component of the module and a signature of gene expression\u0026ndash;was then correlated with the dose, and each module that was correlated with the dose-response curve with a p-value \u0026lt; 0.01 (p-value \u0026lt; 0.05 for the consensus network) was considered statistically significant.\u003c/p\u003e"},{"header":"List Of Abbreviations","content":"\u003cp\u003eQTP: Qinghai-Tibetan Plateau\u003c/p\u003e\n\u003cp\u003eDMRs: differentially methylated regions\u003c/p\u003e\n\u003cp\u003eWGBS: whole genome bisulfite sequencing\u003c/p\u003e\n\u003cp\u003eFPKM: fragments per kilobase of transcript per million fragments\u003c/p\u003e\n\u003cp\u003eWGCNA: Weighted gene correlation network analysis\u003c/p\u003e\n\u003cp\u003eTOM: topological overlap metric\u003c/p\u003e\n\u003cp\u003eGO: Gene\u0026nbsp;Ontology\u003c/p\u003e\n\u003cp\u003eKEGG: Kyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n\u003cp\u003eA-DMRs: age-related DMRs\u003c/p\u003e\n\u003cp\u003eER: endoplasmic reticulum\u003c/p\u003e\n\u003cp\u003eGS: gene significance\u003c/p\u003e\n\u003cp\u003eMM: module membership measures\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll protocols for collection of the semen samples of yaks were reviewed and approved by the Ethics Committee at Institute of Animal Science and Veterinary, Tibet Academy of Agricultural and Animal Husbandry Sciences (Permit Number: 2015-216).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and material\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe DNA methylation data and RNA transcriptome data in this study are available in SRA under the accession numbers PRJNA530286 and PRJNA512958, respectively.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by a program of Provincial Department of Finance of the Tibet Autonomous Region (No: XZNKY-2019-C-052), Program National Beef Cattle and Yak Industrial Technology System (No: CARS-37), Basic Research Programs of Sichuan Province (No: 2019YJ0256) and the Open Project Program of State Key Laboratory of Hulless Barley and Yak Germplasm Resources and Genetic Improvement (NO:XZNKY-2019-C-007K10). The funding bodies had no role in the study design; collection, analysis, and interpretation of data; or in writing the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors' contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eJX, QJ and JZ planned and coordinated the study and wrote the manuscript. CY, XC and HJ collected the samples. ZC and CZ performed the library construction and sequencing and the quality control analysis. QZ, YZ and HC performed downstream analysis of the data and assisted in the generation of additional files for the manuscript. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank\u0026nbsp; the animal husbandry station of Chang Du and the agricultural bureau of Riwoqe\u0026nbsp;for all their support.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eWiener G, Han J-L, Long R-J: \u003cstrong\u003eThe Yak\u003c/strong\u003e, 2nd edn. Bangkok, Thailand: Regional Office for Asia and the Pacific, Food and Agriculture Organization of the United Nations; 2003.\u003c/li\u003e\n\u003cli\u003eQiu Q, Zhang G, Ma T, Qian W, Wang J, Ye Z, Cao C, Hu Q, Kim J, Larkin DM: \u003cstrong\u003eThe yak genome and adaptation to life at high altitude\u003c/strong\u003e. \u003cem\u003eNature Genetics \u003c/em\u003e2012, \u003cstrong\u003e44\u003c/strong\u003e(8):946-949.\u003c/li\u003e\n\u003cli\u003eNichols K, Doelman J, Kim JJM, Carson M, Metcalf JA, Cant JP: \u003cstrong\u003eExogenous\u003c/strong\u003e \u003cstrong\u003eessential amino acids\u003c/strong\u003e\u003cstrong\u003estimulate an adaptive unfolded protein response in the mammary glands of lactating cows\u003c/strong\u003e. \u003cem\u003eJ Dairy Sci\u003c/em\u003e 2017, \u003cstrong\u003e100\u003c/strong\u003e(7):5909-5921.\u003c/li\u003e\n\u003cli\u003eKanwar JR, Kanwar RK, Sun X, Punj V, Matta H, Morley SM, Parratt A, Puri M, Sehgal R: \u003cstrong\u003eMolecular and biotechnological advances in\u003c/strong\u003e\u003cstrong\u003emilk\u003c/strong\u003e\u0026nbsp;\u003cstrong\u003eproteins in relation to\u003c/strong\u003e\u0026nbsp;\u003cstrong\u003ehuman\u003c/strong\u003e\u0026nbsp;\u003cstrong\u003ehealth. \u003c/strong\u003e\u003cem\u003eCurr Protein Pept Sci\u003c/em\u003e 2009, \u003cstrong\u003e10\u003c/strong\u003e(4):308-38.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eOsorio JS, Lohakare J, Bionaz M: \u003cstrong\u003eBiosynthesis of\u003c/strong\u003e\u003cstrong\u003emilk\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003efat,\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eprotein, and lactose: roles of transcriptional and posttranscriptional regulation\u003c/strong\u003e. \u003cem\u003ePhysiol Genomics\u003c/em\u003e 2016,\u003cstrong\u003e 48\u003c/strong\u003e(4):231-56.\u003c/li\u003e\n\u003cli\u003eGene Ontology\u0026nbsp;Consortium: \u003cstrong\u003eExpansion of the\u003c/strong\u003e\u003cstrong\u003eGene Ontology\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eknowledgebase and resources\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 2017, \u003cstrong\u003e45\u003c/strong\u003e(D1):D331-D338.\u003c/li\u003e\n\u003cli\u003eFinn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, Chang HY, Doszt\u0026aacute;nyi Z, El-Gebali S, Fraser M, Gough J, Haft D, Holliday GL, Huang H, Huang X, Letunic I, Lopez R, Lu S, Marchler-Bauer A, Mi H, Mistry J, Natale DA, Necci M, Nuka G, Orengo CA, Park Y, Pesseat S, Piovesan D, Potter SC, Rawlings ND, Redaschi N, Richardson L, Rivoire C, Sangrador-Vegas A, Sigrist C, Sillitoe I, Smithers B, Squizzato S, Sutton G, Thanki N, Thomas PD, Tosatto SC, Wu CH, Xenarios I, Yeh LS, Young SY, Mitchell AL: \u003cstrong\u003eInterPro\u0026nbsp;in 2017-beyond protein family and domain annotations\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 2017, \u003cstrong\u003e45\u003c/strong\u003e(D1):D190-D199.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eKanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K: \u003cstrong\u003eKEGG: new perspectives on genomes, pathways, diseases and drugs\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 2017, \u003cstrong\u003e45\u003c/strong\u003e(D1):D353-D361.\u003c/li\u003e\n\u003cli\u003eBoutet E, Lieberherr D, Tognolli M, Schneider M, Bansal P, Bridge AJ, Poux S, Bougueleret L, Xenarios I: \u003cstrong\u003eUniProtKB/Swiss-Prot, the Manually Annotated Section of the UniProt KnowledgeBase: How to Use the Entry View\u003c/strong\u003e. \u003cem\u003eMethods Mol Biol \u003c/em\u003e2016, \u003cstrong\u003e1374\u003c/strong\u003e:23-54.\u003c/li\u003e\n\u003cli\u003eJunker V, Contrino S, Fleischmann W, Hermjakob H, Lang F, Magrane M, Martin MJ, Mitaritonna N, O'Donovan C, Apweiler R: \u003cstrong\u003eThe role SWISS-PROT and TrEMBL play in the genome research environment.\u003c/strong\u003e \u003cem\u003eJ Biotechnol\u003c/em\u003e 2000, \u003cstrong\u003e78\u003c/strong\u003e(3):221-34.\u003c/li\u003e\n\u003cli\u003eBurgess, D.J: \u003cstrong\u003eGene expression: Principles of gene regulation across tissues. \u003c/strong\u003e\u003cem\u003eNature Reviews Genetics, \u003c/em\u003e2017, \u003cstrong\u003e18\u003c/strong\u003e(12):701-701.\u003c/li\u003e\n\u003cli\u003ePreston GM, Brodsky JL: \u003cstrong\u003eThe evolving role of ubiquitin modification in\u003c/strong\u003e \u003cstrong\u003eendoplasmic reticulum-associated degradation\u003c/strong\u003e. \u003cem\u003eBiochem J\u003c/em\u003e 2017, \u003cstrong\u003e474\u003c/strong\u003e(4):445-469.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eMcCaffrey K, Braakman I: \u003ca href=\"https://www.ncbi.nlm.nih.gov/pubmed/27744338\"\u003e\u003cstrong\u003eProtein\u0026nbsp;quality control at the\u0026nbsp;endoplasmic reticulum\u003c/strong\u003e.\u003c/a\u003e \u003cem\u003eEssays Biochem\u003c/em\u003e 2016, \u003cstrong\u003e60\u003c/strong\u003e(2):227-235\u003c/li\u003e\n\u003cli\u003eLangfelder P, Horvath S: WGCNA: \u003cstrong\u003ean R package for weighted correlation network analysis\u003c/strong\u003e.\u003cem\u003e BMC Bioinformatics\u003c/em\u003e 2008, \u003cstrong\u003e9\u003c/strong\u003e:559.\u003c/li\u003e\n\u003cli\u003eLiu X, Hu AX, Zhao JL, Chen FL:\u003cstrong\u003e Identification of Key Gene Modules in Human Osteosarcoma by Co-Expression Analysis Weighted Gene Co-Expression Network Analysis (WGCNA)\u003c/strong\u003e.\u003cem\u003e J Cell Biochem\u003c/em\u003e 2017, \u003cstrong\u003e118\u003c/strong\u003e(11):3953-3959.\u003c/li\u003e\n\u003cli\u003eHu H, Miao YR, Jia LH, Yu QY, Zhang Q, Guo AY: \u003cstrong\u003eAnimalTFDB\u0026nbsp;3.0: a comprehensive resource for annotation and prediction of animal transcription factors\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 2019, \u003cstrong\u003e7\u003c/strong\u003e(D1):D33-D38.\u003c/li\u003e\n\u003cli\u003eSun LD, Xiao FL, Li Y, Zhou WM, Tang HY, Tang XF, Zhang H, Schaarschmidt H, Zuo XB, Foelster-Holst R, He SM, Shi M, Liu Q, Lv YM, Chen XL, Zhu KJ, Guo YF, Hu DY, Li M, Li M, Zhang YH, Zhang X, Tang JP, Guo BR, Wang H, Liu Y, Zou XY, Zhou FS, Liu XY, Chen G, Ma L, Zhang SM, Jiang AP, Zheng XD, Gao XH, Li P, Tu CX, Yin XY, Han XP, Ren YQ, Song SP, Lu ZY, Zhang XL, Cui Y, Chang J, Gao M, Luo XY, Wang PG, Dai X, Su W, Li H, Shen CP, Liu SX, Feng XB, Yang CJ, Lin GS, Wang ZX, Huang JQ, Fan X, Wang Y, Bao YX, Yang S, Liu JJ, Franke A, Weidinger S, Yao ZR, Zhang XJ: \u003cstrong\u003eGenome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population\u003c/strong\u003e. \u003cem\u003eNat Genet\u003c/em\u003e 2011, \u003cstrong\u003e43\u003c/strong\u003e(7):690-4.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eLiao D: \u003cstrong\u003eEmerging roles of the\u0026nbsp;EBF\u0026nbsp;family of transcription factors in tumor\u0026nbsp;suppression\u003c/strong\u003e. \u003cem\u003eMol Cancer Res\u003c/em\u003e 2009, \u003cstrong\u003e7\u003c/strong\u003e(12):1893-901.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eJargosch M, Kr\u0026ouml;ger S, Gralinska E, Klotz U, Fang Z, Chen W, Leser U, Selbig J, Groth D, Baumgrass R: \u003cstrong\u003eData integration for identification of important transcription factors of\u0026nbsp;STAT6-mediated cell fate decisions\u003c/strong\u003e. \u003cem\u003eGenet Mol Res\u003c/em\u003e 2016, \u003cstrong\u003e15\u003c/strong\u003e(2).\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eMoniz S, Biddlestone J, Rocha S: \u003cstrong\u003eGrow2: the\u0026nbsp;HIF\u0026nbsp;system,\u0026nbsp;energy\u0026nbsp;homeostasis and the cell cycle\u003c/strong\u003e. \u003cem\u003eHistol Histopathol\u003c/em\u003e 2014, \u003cstrong\u003e29\u003c/strong\u003e(5):589-600.\u003c/li\u003e\n\u003cli\u003evan der Bliek AM, Sedensky MM, Morgan PG: \u003cstrong\u003eCell Biology of the\u003c/strong\u003e\u003cstrong\u003eMitochondrion\u003c/strong\u003e. \u003cem\u003eGenetics\u003c/em\u003e 2017, \u003cstrong\u003e207\u003c/strong\u003e(3):843-871.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eGanesan V, Sivanesan D, Yoon S:\u003cstrong\u003e Correlation between the Structure and Catalytic Activity of [Cp*Rh(Substituted Bipyridine)] Complexes for NADH Regeneration\u003c/strong\u003e. \u003cem\u003eInorg Chem\u003c/em\u003e 2017, \u003cstrong\u003e56\u003c/strong\u003e(3):1366-1374.\u0026nbsp;\u003c/li\u003e\n\u003cli\u003eLi T, Liu J, Smith WW: \u003cstrong\u003eSynphilin-1 binds\u0026nbsp;ATP\u0026nbsp;and regulates intracellular energy status\u003c/strong\u003e.\u003cem\u003e PLoS One\u003c/em\u003e 2014, \u003cstrong\u003e9\u003c/strong\u003e(12):e115233.\u003c/li\u003e\n\u003cli\u003eGoody RS: \u003cstrong\u003eThe significance of the free\u0026nbsp;energy\u0026nbsp;of hydrolysis of\u0026nbsp;GTP\u0026nbsp;for signal-transducing and regulatory GTPases\u003c/strong\u003e. \u003cem\u003eBiophys Chem\u003c/em\u003e 2003, \u003cstrong\u003e100\u003c/strong\u003e(1-3):535-44.\u003c/li\u003e\n\u003cli\u003eNunes-Nesi A, Ara\u0026uacute;jo WL, Obata T, Fernie AR:\u003cstrong\u003e Regulation of the mitochondrial tricarboxylic acid cycle\u003c/strong\u003e. \u003cem\u003eCurr Opin Plant Biol\u003c/em\u003e 2013, \u003cstrong\u003e16\u003c/strong\u003e(3):335-43.\u003c/li\u003e\n\u003cli\u003eBal NC, Singh S, Reis FCG, Maurya SK, Pani S, Rowland LA, Periasamy M: \u003cstrong\u003eBoth brown adipose tissue and skeletal muscle\u003c/strong\u003e\u003cstrong\u003ethermogenesis\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eprocesses are activated during mild to severe cold adaptation in mice\u003c/strong\u003e. \u003cem\u003eJ Biol Chem\u003c/em\u003e 2017, \u003cstrong\u003e292\u003c/strong\u003e(40):16616-16625.\u003c/li\u003e\n\u003cli\u003eXiaogang Cui,\u0026nbsp;Yali Hou,\u0026nbsp;Shaohua Yang,\u0026nbsp;Yan Xie,\u0026nbsp;Shengli Zhang,\u0026nbsp;Yuan Zhang,\u0026nbsp;Qin Zhang,\u0026nbsp;Xuemei Lu,\u0026nbsp;George E Liu,\u0026nbsp;Dongxiao Sun: \u003cstrong\u003eTranscriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing\u003c/strong\u003e.\u003cem\u003eBMC Genomics\u003c/em\u003e. 2014, \u003cstrong\u003e15\u003c/strong\u003e:226.\u003c/li\u003e\n\u003cli\u003eSu\u0026aacute;rez-Vega A, Guti\u0026eacute;rrez-Gil B, Klopp C, Tosser-Klopp G, Arranz JJ: \u003cstrong\u003eComprehensive RNA-Seq profiling to evaluate lactating sheep mammary gland transcriptome\u003c/strong\u003e.\u003cem\u003eSci Data\u003c/em\u003e. 2016, \u003cstrong\u003e3\u003c/strong\u003e:160051.\u003c/li\u003e\n\u003cli\u003eMobuchon L, Marthey S, Boussaha M, Le Guillou S, Leroux C, Le Provost F: \u003cstrong\u003eAnnotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches\u003c/strong\u003e.\u0026nbsp;\u003cem\u003eBMC Genomics\u003c/em\u003e. 2015, \u003cstrong\u003e16\u003c/strong\u003e(1):285.\u003c/li\u003e\n\u003cli\u003eYang Zhou,\u0026nbsp;Erin E Connor,\u0026nbsp;Derek M Bickhart, Congjun Li,\u0026nbsp;Ransom L Baldwin, Steven G Schroeder, Benjamin D Rosen, Liguo Yang,\u0026nbsp;Curtis P Van Tassell,\u0026nbsp;George E Liu: \u003cstrong\u003eComparative whole genome DNA methylation profiling of cattle sperm and somatic tissues reveals striking hypomethylated patterns in sperm\u003c/strong\u003e.\u0026nbsp;\u003cem\u003eGigascience.\u003c/em\u003e 2018, \u003cstrong\u003e7\u003c/strong\u003e(5):giy039.\u003c/li\u003e\n\u003cli\u003e\u003ca href=\"https://pubmed.ncbi.nlm.nih.gov/?term=Li+Z\u0026amp;cauthor_id=31276563\"\u003eZhixiong Li\u003c/a\u003e,\u0026nbsp;\u003ca href=\"https://pubmed.ncbi.nlm.nih.gov/?term=Jiang+M\u0026amp;cauthor_id=31276563\"\u003eMingfeng Jiang\u003c/a\u003e: \u003cstrong\u003eMetabolomic profiles in yak mammary gland tissue during the lactation cycle\u003c/strong\u003e.\u003cem\u003e PLoS One \u003c/em\u003e2019, \u003cstrong\u003e14\u003c/strong\u003e(7):e0219220.\u003c/li\u003e\n\u003cli\u003eQiumei Ji, Jinwei Xin, Zhixin Chai, Chengfu Zhang, Dawa Yangla, Luo Sang, Qiang Zhang, Pingcuo Zhandui, Minsheng Peng, Yong Zhu, Hanwen Cao, Hui Wang, Jianlin Han, Jincheng Zhong:\u003cstrong\u003e A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak\u003c/strong\u003e. \u003cem\u003eMolecular Ecology Resources\u003c/em\u003e, DOI:10.1111/1755-0998.13236. In press.\u003c/li\u003e\n\u003cli\u003eXi Y, Li W: \u003cstrong\u003eBSMAP: whole genome bisulfite sequence MAPping program\u003c/strong\u003e. \u003cem\u003eBMC Bioinformatics \u003c/em\u003e2009, \u003cstrong\u003e10\u003c/strong\u003e(1):232-232.\u003c/li\u003e\n\u003cli\u003eLister R, Pelizzola M, Dowen RH, Hawkins RD, Hon GC, Tontifilippini J, Nery JR, Lee LK, Ye Z, Ngo Q: \u003cstrong\u003eHuman DNA methylomes at base resolution show widespread epigenomic differences\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2009, \u003cstrong\u003e462\u003c/strong\u003e(7271):315-322.\u003c/li\u003e\n\u003cli\u003eJuhling F, Kretzmer H, Bernhart SH, Otto C, Stadler PF, Hoffmann S: \u003cstrong\u003emetilene: Fast and sensitive calling of differentially methylated regions from bisulfite sequencing data\u003c/strong\u003e. \u003cem\u003eGenome Research \u003c/em\u003e2016, \u003cstrong\u003e26\u003c/strong\u003e(2):256-262.\u003c/li\u003e\n\u003cli\u003eWang, J., Duncan, D., Shi, Z., Zhang, B: \u003cstrong\u003eWEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update. \u003c/strong\u003e\u003cem\u003eNucleic Acids Res \u003c/em\u003e\u003cstrong\u003e2013\u003c/strong\u003e, 41: W77-83.\u003c/li\u003e\n\u003cli\u003ePertea M, Pertea G, Antonescu C, Chang TC, Mendell JT, Salzberg SL: \u003cstrong\u003eStringTie enables improved reconstruction of a transcriptome from RNA-seq reads\u003c/strong\u003e. \u003cem\u003eNature Biotechnology \u003c/em\u003e2015, \u003cstrong\u003e33\u003c/strong\u003e(3):290-295.\u003c/li\u003e\n\u003cli\u003eRobinson MD, Mccarthy DJ, Smyth GK: \u003cstrong\u003eedgeR: a Bioconductor package for differential expression analysis of digital gene expression data\u003c/strong\u003e. \u003cem\u003eBioinformatics \u003c/em\u003e2010, 26(1):139-140.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"milk production, DNA methylation, transcription factor, epigenetic regulation","lastPublishedDoi":"10.21203/rs.3.rs-20775/v4","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-20775/v4","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eDomestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), by providing milk and meat. They have evolved numerous physiological adaptations to high-altitude environment, including strong blood oxygen transportation capabilities and high metabolism. The roles of DNA methylation and gene expression in milk production and high-altitudes adaptation need further exploration. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults \u003c/strong\u003e\u003c/p\u003e\u003cp\u003eWe performed genome-wide DNA methylome and transcriptome analyses of breast, lung, and biceps brachii muscle tissues from yaks of different ages. We identified 432,350 differentially methylated regions (DMRs) across the age groups within each tissue. The post-mature breast tissue had considerably more differentially methylated regions (155,957) than that from the three younger age groups. Hypomethylated genes with high expression levels might regulate milk production by influencing protein processing in the endoplasmic reticulum. According to weighted gene correlation network analysis, the “hub” gene ZGPAT was highly expressed in the post-mature breast tissue, indicating that it potentially regulates the transcription of 280 genes that influence\u0026nbsp;protein\u0026nbsp;synthesis, processing, and secretion. The tissue network analysis indicated that high expression of HIF1A regulates energy metabolism in the lung. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eThis study provides a basis for understanding the epigenetic mechanisms underlying milk production in yaks, and the results offer insight to breeding programs aimed at improving milk production.\u003c/p\u003e","manuscriptTitle":"Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production","msid":"","msnumber":"","nonDraftVersions":[{"code":"","date":"2020-10-20 21:32:35","doi":"","editorialEvents":[{"type":"checksComplete","content":"","date":"2020-10-10T12:00:00+00:00","index":"","fulltext":""},{"type":"decision","content":"Minor revision","date":"2020-10-09T12:00:00+00:00","index":"","fulltext":""},{"type":"notPreprinted","content":""}],"status":"timeline","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":4,"date":"2020-09-29 21:32:35","doi":"10.21203/rs.3.rs-20775/v4","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2020-09-28T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-09-24T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-09-23T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-09-23T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":3,"date":"2020-09-11 23:04:15","doi":"10.21203/rs.3.rs-20775/v3","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2020-09-14T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-09-08T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-09-07T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-09-07T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":2,"date":"2020-08-11 19:20:15","doi":"10.21203/rs.3.rs-20775/v2","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2020-08-26T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-08-24T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2020-08-21T12:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-08-10T12:00:00+00:00","index":2,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-08-08T12:00:00+00:00","index":1,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2020-08-07T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-08-06T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-08-05T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-08-05T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":1,"date":"2020-04-03 14:55:36","doi":"10.21203/rs.3.rs-20775/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2020-06-15T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-05-25T12:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-05-06T12:00:00+00:00","index":2,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-05-06T12:00:00+00:00","index":3,"fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-05-05T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-04-08T12:00:00+00:00","index":1,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2020-04-06T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-03-09T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-03-08T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-03-08T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"cafcf329-debe-407f-a330-f6862bffd076","owner":[],"postedDate":"September 29th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":497689,"name":"Epigenetics \u0026 Genomics"},{"id":497690,"name":"Animal Science"}],"tags":[],"updatedAt":"2020-10-25T15:02:47+00:00","versionOfRecord":{"articleIdentity":"rs-20775","link":"https://doi.org/10.1186/s12864-020-07151-3","journal":{"identity":"bmc-genomics","isVorOnly":false,"title":"BMC Genomics"},"publishedOn":"2020-10-20 12:00:00","publishedOnDateReadable":"October 20th, 2020"},"versionCreatedAt":"2020-09-29 21:32:35","video":"","vorDoi":"10.1186/s12864-020-07151-3","vorDoiUrl":"https://doi.org/10.1186/s12864-020-07151-3","workflowStages":[]},"version":"v4","identity":"rs-20775","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-20775","identity":"rs-20775","version":["v4"]},"buildId":"wLkW0s4AflPzk-lpfg-fK","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}