Reproductive patterns of the yak. I. Reproductive phenomena of the female yak

Serum Metabolomic Analysis of Synchronous Estrus in Yaks Based on UPLC-Q-TOF MS Technology

The yak is a unique species of livestock found in the Qinghai-Tibet Plateau and its surrounding areas. Due to factors such as late sexual maturity and a low rate of estrus, its reproductive efficiency is relatively low. The process of estrus synchronization in yaks plays a crucial role in enhancing their reproductive success and ensuring the continuation of their species. In order to clarify the characteristics of the serum metabolites of yak estrus synchronization, the yaks with inactive ovaries were compared with the estrus synchronization yaks. In this study, yaks were divided into the inactive ovaries group (IO), gonarelin-induced yak estrus group (GnRH), and chloprostenol sodium-induced yak estrus group (PGF). After the completion of the estrus synchronization treatment, blood samples were collected from the jugular veins of the non-estrus yaks in the control group and the yaks with obvious estrus characteristics in the GnRH and PGF groups. Metabolites were detected by ultra-high performance liquid chromatography-mass spectrometry, and differential metabolites were screened by multivariate statistical analysis. The results showed that a total of 70 significant differential metabolites were screened and identified in the GnRH vs. IO group, and 77 significant differential metabolites were screened and identified in the PGF vs. IO group. Compared with non-estrus yaks, 36 common significant differential metabolites were screened out after the induction of yak estrus by gonarelin (GnRH) and cloprostenol sodium (PGF), which were significantly enriched in signaling pathways such as the beta oxidation of very long chain fatty acids, bile acid biosynthesis, oxidation of branched chain fatty acids, steroidogenesis, steroid biosynthesis, and arginine and proline metabolism. This study analyzed the effects of gonadotropin releasing hormone (GnRH) and prostaglandin F (PGF) on the reproductive performance of yaks treated with estrus synchronization, which provides a theoretical basis for the optimization and application of yak estrus synchronization technology and promotes the healthy development of the yak industry.

"The Yak"-A remarkable animal living in a harsh environment: An overview of its feeding, growth, production performance, and contribution to food security

Yaks play an important role in the livelihood of the people of the Qinghai-Tibet Plateau (QTP) and contribute significantly to the economy of the different countries in the region. Yaks are commonly raised at high altitudes of ~ 3,000-5,400 m above sea level. They provide many important products, namely, milk, meat, fur, and manure, as well as social status, etc. Yaks were domesticated from wild yaks and are present in the remote mountains of the QTP region. In the summer season, when a higher quantity of pasture is available in the mountain region, yaks use their long tongues to graze the pasture and spend ~ 30-80% of their daytime grazing. The remaining time is spent walking, resting, and doing other activities. In the winter season, due to heavy snowfall in the mountains, pasture is scarce, and yaks face feeding issues due to pasture scarcity. Hence, the normal body weight of yaks is affected and growth retardation occurs, which consequently affects their production performance. In this review article, we have discussed the domestication of yaks, the feeding pattern of yaks, the difference between the normal and growth-retarded yaks, and also their microbial community and their influences. In addition, blood biochemistry, the compositions of the yaks' milk and meat, and reproduction are reported herein. Evidence suggested that yaks play an important role in the daily life of the people living on the QTP, who consume milk, meat, fur, use manure for fuel and land fertilizer purposes, and use the animals for transportation. Yaks' close association with the people's well-being and livelihood has been significant.

The Effects of Postpartum Yak Metabolism on Reproductive System Recovery

The goal of this study was to determine the metabolism of multiparous female yaks during the late perinatal period and identify its effects on reproductive recovery in order to explain the low reproduction rate of yaks. Eight multiparous female yaks were randomly selected as the sample, and serum was collected from the yaks every 7 days from the day of delivery until 28 days after the delivery (five time points). The presence of serum metabolic profiles and reproductive hormones was identified using ELISA. The key metabolites were identified using liquid chromatography-mass spectrometry, and a dynamic metabolic network representation was created using bioinformatics analysis. A total of 117 different metabolites were identified by calculating the fold change of the metabolite expression at each time point. The dynamic metabolic network was created to represent the activities of the key metabolites, metabolic indexes and reproductive hormones. The initial efficiency of the glucose metabolism in the late perinatal period was found to be low, but it increased during the final period. The initial efficiencies of the lipid and amino acid metabolisms were high but decreased during the final period. We inferred that there was a postpartum negative energy balance in female yaks and that the synthesis and secretion of estrogen were blocked due to an excessive fatty acid mobilization. As a result, the reproductive hormone synthesis and secretion were maintained at a low level in the late perinatal period, and this was the main reason for the delayed recovery of the reproductive function postpartum. However, the specific mechanism needs to be further verified.

A comparative transcriptome and proteomics study of post-partum ovarian cycle arrest in yaks (Bos grunniens)

Post-partum ovarian cycle arrest is the main factor affecting yak reproductive efficiency. There are few reports regarding the molecular regulatory mechanism of post-partum oestrus at transcriptome and proteome levels in yaks. Our previous studies focussed on the ovaries of yaks with post-partum ovarian cycle arrest and post-partum oestrus yaks. In this study, RNA sequencing transcriptomic study was combined with quantitative proteomic analyses to identify post-partum ovarian cycle-related genes and proteins. Consequently, 1,149 genes and 24 proteins were found to be up- or downregulated during post-partum oestrus. The analysis of differentially regulated genes identified three gene or protein pairs that were synchronously upregulated and no gene or protein pairs that were synchronously downregulated, suggesting that these upregulated genes may regulate the post-partum ovarian cycle. The functional classification of these differentially expressed genes and proteins indicated their connection with the oocyte meiosis, the oestrogen signalling pathway, the progesterone-mediated oocyte maturation and the gonadotrophin-releasing hormone (GnRH) signalling pathway. In this study, a total of six genes and two proteins involved in the oocyte meiosis, the oestrogen signalling pathway, the progesterone-mediated oocyte maturation and the GnRH signalling pathway were identified. The CSNK1A1, M91_09723, M91_11326, M91_21439, M91_19073, SHC2, Atf6b, M91_03062, HSPCA and calmodulin could regulate oestrus, respectively, in the post-partum so as to control the anoestrus status.

The Distribution, Expression Patterns and Functional Analysis of NR1D1 and NR4A2 in the Reproductive Axis Tissues of the Male Tianzhu White Yak

Nuclear hormone receptors NR1D1 and NR4A2 play important roles in the synthesis and metabolism of hormones that are thought to be strictly regulated by the hypothalamus-pituitary-gonad axis (HPG) tissues via gene expression. However, in the yak, the function and regulatory mechanisms of NR1D1 and NR4A2 are not clearly understood. The current study is aimed to investigate the expression patterns, distribution and functions of these two receptors in HPG tissues in male Tianzhu white yaks. Immunohistochemical staining showed NR1D1 and NR4A2 proteins were present in all yak HPG tissues with differential expression patterns and degrees of staining, particularly in Leydig cells that were strongly positive in accordance with the immunofluorescence results. qRT-PCR and Western blot results suggested that the highest expression levels of NR1D1 and NR4A2 mRNA were present in the hypothalamus, while the expression levels of NR1D1 and NR4A2 proteins were higher in the testis and epididymis than in the hypothalamus or pituitary gland. In addition, expression levels of NR1D1 and NR4A2 mRNA and protein in testicular tissues differed by age. Expression levels were significantly higher at 6 years of age. Gene ontology (GO) and pathway analysis enrichment revealed that NR1D1 may directly regulate the synthesis and metabolism of steroid hormones via interaction with different targets, while NR4A2 may indirectly regulate the synthesis and metabolism of steroid hormones. These results showed that NR1D1 and NR4A2, as important mediators, are involved in the regulation of male yak reproduction, and especially of steroid hormones and androgen metabolism. These results will be helpful for the further understanding of the regulatory mechanisms of NR1D1 and NR4A2 in yak reproduction.

Expression characteristics of immune factors in the yak (Bos grunniens) testis

The goal of this study was to characterize and evaluate the main markers of macrophages, T lymphocytes, B lymphocytes and plasmocytes in the testis of juvenile and adult yaks by quantitative real-time polymerase chain reaction and immunohistochemistry. Within the same age group, the mRNA expression of CD68 was always highest, followed by that of CD3ε, CD79α, IgG and IgA. Moreover, CD68, CD3, CD79α, IgA and IgG positive cells were all located in the testicular interstitial tissues of juvenile and adult yaks. In the same age group, the frequency of CD68 positive macrophages was higher than that of CD3 positive T lymphocytes, which was followed by that of CD79α positive B lymphocytes and IgA and IgG positive plasmocytes. No significant difference was observed between the B lymphocyte and plasmocyte frequencies in yak testes. Furthermore, CD68, CD3ε, CD79α, IgA and IgG mRNA expression levels and the frequencies of CD68, CD3, CD79α, IgA and IgG positive cells increased from juveniles to adults. Similarly, the frequencies of CD68, CD3, CD79α, IgA and IgG positive cells also increased with age. These results suggest that in the yak testis, the immune defence system against pathogens might primarily comprise macrophages and T lymphocytes in the testicular interstitial tissue. Moreover, the testicular immune environment may mature and expand to a fully functional state in adult yaks. The low frequencies of B lymphocyte and plasmocyte in yaks, differing from those in rodents and humans, might be related to the fact that yaks live in low-oxygen plateaus.

Dynamics of apoptosis-related gene expression during follicular development in yak

The potential reproduction power of domestic animals is limited by a complicated follicular atresia process. P53, caspase-9 (Casp9), Bax, Bcl-2 and Fas play a crucial role in the ovarian mitochondrion-dependent apoptosis and death receptor pathway. In accordance with this study, the expression levels of Casp9, Bax, Bcl-2 and Fas were analysed in ovaries and oviducts of yak by immunohistochemistry (IHC). P53 and the above in ovarian granulosa cells (GCs) from atretic (3-6 mm) to healthy follicles (6-8 mm) and in oviducts were examined from the luteal phase to the follicular phase during the oestrous circle by Western blot (WB) and real-time PCR (RT-PCR). Results demonstrated that typical classic apoptotic factors Casp9, Bax, Bcl-2 and Fas were expressed in the cytoplasm and zonal pellucida of oocytes, primordial follicles, primary follicles, ovarian surface epithelium, ovarian GCs, granular lutein cells, surface epithelia in oviduct uterotubal junction and oviduct ampulla during the luteal phase. RT-PCR and WB revealed that P53 and Fas significantly increased in GCs of atretic follicles. P53 and Casp9 increased in oviduct epithelium during the luteal phase, but Fas was unchanged. A contrary tendency was noted in Bcl-2 and Bax expression. Overall, P53 and Fas play an essential role in inducing GC apoptosis, and Bax, Bcl-2, Casp9 and P53 are involved in oviduct epithelial regeneration in yak.

Transcriptome analysis revealed key signaling networks regulating ovarian activities in the domestic yak

Domestic yaks are the most important livestock species on the Qinghai-Tibetan Plateau. Adult female yaks normally breed in warm season and enter anestrous in cold season. Currently, how the ovarian activity is regulated at the molecular level remains to be determined. This study was conducted to investigate follicular development and gene expression patterns of yak ovarian tissues in the warm and cold seasons. Dynamics of follicular development was evaluated based on histological analyses and global gene expression was examined by using RNA-sequencing (RNA-seq) technology. Firstly, we found that follicle development of yak cows in cold season was different from that in warm season. Interestingly, ovaries collected from yaks in cold season contained a significant higher number of antral follicles and some of these follicles showed signs of polycystic structure, indicating abnormal granulosa cell function. RNA-seq analyses of ovarian tissues from non-pregnant adult yaks in cold and warm season revealed that a list of 320 transcripts were differentially expressed, specifically, 79 were up-regulated and 241 were down-regulated in the ovaries from yaks during the cold season. Further analysis demonstrated that transcripts associated with estrogen secretion and metabolism signaling pathway were altered, including FST, CYP1A1, PIK3R1 and PIK3R2. This study showed histological features of follicle development and revealed candidate genes that may have important roles in regulating ovarian activities in the yak seasonal reproduction.

Expression and localization of follicle-stimulating hormone receptor in the yak uterus during different stages of the oestrous cycle

Morphological changes of the uterus and alterations in its secretory activity under the influence of steroid hormones been well documented. The oestrous cycle is also associated with significant changes in plasma follicle-stimulating hormone (FSH), whose effects are mediated through its receptor (FSHR). Reports showed that in many mammals, FSHR was expressed in gonadal and extragonadal tissues including cervix, female reproductive tract, and pituitary gland. Follicle-stimulating hormone (FSH) signals through endothelial FSHR directly stimulate angiogenesis and involved in the timing of birth in human, and the presence of FSHR in the placenta is essential for normal pregnancy in mice. But little is known about FSHR expression in the yak uterus. The main objective of the present study was to determine the expression and localization of FSHR in the yak uterus during different phases of the oestrous cycle. Results showed that FSHR protein was localized in the surface and glandular epithelial cells, stroma cells, myometrial smooth muscle cells and blood vessel endothelial cells. The expression of FSHR protein peaked at oestrus, significantly decreased at dioestrus (p < 0.05) and increased again at proestrus. FSHR mRNA was highly expressed at both proestrus and oestrus, and decreased at metestrus with the lowest values at dioestrus (p < 0.05). In conclusion, FSHR expression in the yak uterus changed with the stage of the oestrous cycle suggesting that FSHR plays an essential role in regulating the endometrial and myometrial functions during the oestrus cycle in the yak.

The histological characteristics, age-related thickness change of skin, and expression of the HSPs in the skin during hair cycle in yak (Bos grunniens)

OBJECTIVE

This experiment was conducted to study the histological characteristics, age-related thickness changes, and expression of HSPs in the skin of yak.

METHODS

A total of 20 yaks (10 males and 10 females) were used. Different regions of the normal skin of three different ages (newborn, half-year-old and adult) of yaks were harvested for histological study and thickness measurement. Biopsy samples were taken from the scapula regions of the skin from the same five approximately 1-year-old yaks during the hair cycle (telogen, anagen and catagen). RT-PCR, western blot and immunohistochemistry methods using the mRNA and protein levels were used to detect the expression of HSP27, HSP70 and HSP90. RT-PCR method was used to detect the mRNA expression of CGI-58 and KDF1. The IPP6.0 software was used to analyze the immunohistochemistry and measure the thickness of the skin.

RESULTS

The general histological structure of hairy yak skin was similar to other domestic mammals. The unique features included prominent cutaneous vascular plexuses, underdeveloped sweat glands, a large number of nasolabial glands in the nasolabial plate, and hair follicle groups composed of one primary follicle and several secondary follicles. The skin, epidermis and dermis thickness did vary significantly between different body regions and different ages. The thickness of the skin, epidermis and dermis increased from newborn to adult in yaks. Yak skin thickness decreased from dorsally to ventrally on the trunk. The skin on the lateral surface was thicker than the skin on the medial surface on the limbs. HSP27, HSP70 and HSP90 showed different expression patterns during the hair cycle using RT-PCR, western blot and immunohistochemistry methods. The expression of HSP27 mRNA and protein in the anagen stage was the highest, followed by the catagen stage, and the expression in the telogen stage was the lowest. The expression of HSP70 mRNA and protein in the telogen stage was the highest, followed by the anagen stage, and the expression in the catagen stage was the lowest. The expression of HSP90 mRNA and protein in the anagen stage was the highest, followed by the telogen stage, and the expression in the catagen stage was the lowest. HSPs were mainly expressed in the outer root sheath of hair follicle during the hair cycle, also expressed in epidermis, sebaceous gland and sweat gland in the skin of Yak. The expression of CGI-58 mRNA in the anagen stage was the highest, followed by the catagen stage, and the expression in the telogen stage was the lowest. The expression of KDF1 mRNA in the telogen stage was the highest, followed by the catagen stage, and the expression in the anagen stage was the lowest.

MEANING

In this study, we examined and fully described the histology of normal skin in Yak and measured the skin thickness of different ages and different regions in Yak. These data may be useful to better understand and appreciate the adaptability features of yak skin. Our investigation reports the expression patterns of HSPs in yak skin for the first time. The different expression pattern of HSPs during the hair cycle suggests they may play different roles in yak hair follicle biology.

Morphometric analysis of yak placentomes during gestation

To study the morphological characteristics of placentomes during gestation, thirty-one yaks were used to determine anatomical and histological data of placentomes by morphometric method. The results showed that there were no significant changes in the total number of yak placentomes during gestation, but remarkably fewer placentomes in the non-pregnant horn compared with the pregnant horn. The mean size of placentomes increased significantly from <61 days to 181-210 days and decreased in 211+ days. The development of placentomes in the pregnant horn were obviously different from that of the non-pregnant horn except <61 days, which were longer, wider, and thicker. The volume density of fetal villi increased during gestation, with a concomitant reduction in the volume density of caruncular endometrium. However, there was no significant change in the surface density of fetal villi. The volume density of binucleate cell increased in prophase and then decreased. No apparent changes in the volume densities of maternal and fetal pyknotic cell occurred throughout gestation. There was a significant increase in the average total volume of yak placentomes from 61-90 days to 181-210 days and a significant decrease in 211+ days. The same patterns were showed in the average total volumes of caruncular endometrium and binucleate cell as well as the mean total surface area of fetal villi. However, the mean total volumes of fetal villi, maternal and fetal pyknotic cell increased constantly with increasing gestation time. Three types of placentome shape were identified in yak. The results showed that the maternal-fetal contact area of yak was enlarged by increasing the placentome size as well as the total volume and surface area of placentome component during pregnancy, so placental function could be enhanced further to meet the needs of fetal development.

Application of sensitive enzymeimmunoassay for determination of cortisol in blood plasma of yaks (Poephagus grunniens L.)

As an alternative to radioimmunoassays, a simple, highly sensitive and quick enzymeimmunoassay (EIA) for determination of cortisol in blood plasma of yaks on microtiterplates using second antibody coating technique and cortisol-horseradish peroxidase as a label has been developed. The wells of the microtiterplate were coated with affinity-purified goat IgG (antirabbit IgG) that binds the hormone specific antibody. The EIA was carried out directly in 20 microl of heat treated plasma after 1:5 dilution with PBS. The cortisol standard curve, with doses ranged from 0.4 to 100 pg/well. The sensitivity of the assay was 20 pg/ml. Cortisol standard curve in buffer showed parallelism with serially diluted yak plasma containing high endogenous cortisol. Intra- and inter-assay coefficients of variation (CV) determined using pooled plasma was found 6.58 and 11.35%, respectively. Recovery of known concentrations of added cortisol in charcoal stripped plasma was linear (r = 0.98). For biological validation of cortisol enzymeimmunoassay, the blood samples were collected from yak cows at -48 and -24h before and 0, 12, 24, 36, 48, 60, 72, 84 and 96 h after dexamethasone administration. The plasma cortisol before dexamethasone administration was significantly (P < 0.05) higher than after dexamethasone administration. The developed EIA was further validated biologically by estimating cortisol in peri-parturient cows beginning day 10 prior to calving till day 10 post-calving; the concentrations were along with the expected lines as reported in bovine. In conclusion, the EIA developed in this study is simple, highly sensitive, valid and sufficiently reliable method for estimation of cortisol directly in bovine plasma.

Is the yak (Poephagus grunniens L.) really a seasonal breeder?

Yaks are considered to be seasonally polyestrous and breeding occurs from July to November. Here we show that some yaks in peak non-breeding season do exhibit cyclic luteal activity without exhibiting any behavioral signs around expected estrus. A total of eight non-lactating yaks were selected from the Yak Farm belonging to National Research Centre on Yak for various sets of experiments. The animals were maintained as per semi range system of management. They were allowed to graze during daytime and fed concentrate mixture @2 kg/animal/day as per standard farm practices of the center. Blood samples were collected on alternate days for 30 days by jugular venipuncture from the yaks during peak breeding season (July to November) and from the same yaks in non-breeding season (February to March). The plasma samples were analysed for progesterone and estradiol-17beta by RIA and EIA procedures, respectively. During breeding season, the mean plasma progesterone at estrus was basal (<or=0.2 ng/ml) and there after started to increase and reached a peak at days 15-16 of the cycle and then declined rapidly to basal levels at estrus. Mean plasma estradiol-17beta concentrations were high at estrus and then declined to basal level at days 2-3 of the cycle. The mean plasma estradiol-17beta concentrations stayed low for the remaining part of the estrous cycle except for a small elevation, which was exhibited between days 6 and 12 of the cycle. During non-breeding season, in 50% of the animal studied, plasma progesterone and estradiol-17beta levels were at basal level. However, in the remaining animals distinct luteal activity was recorded in terms of plasma progesterone profile. High concentrations of progesterone were maintained in one animal, which suggested existence of cystic corpus luteum. The trend of progesterone and estradiol-17beta profiles in other three animals indicated that they could be cycling although they did not exhibit estrus symptoms. It was concluded from the present study that with improved nutritional support some yaks exhibit estrous cyclicity even during the non-breeding period of the year.

Timing of ovulation in relation to onset of estrus and LH peak in yak (Poephagus grunniens L.)

The objective of the study was to determine the timing of ovulation in relation to onset of estrus and the preovulatory LH peak in yaks. For this purpose, a sensitive LH enzymeimmunoassay previously established in buffaloes was successfully validated for measuring the hormone in yak plasma. Plasma LH and progesterone were estimated from blood samples collected from eight non-lactating cycling yaks at 2 h intervals after estrus onset until 6 h after ovulation (ovulation was confirmed by palpation of ovaries per rectum). The mean+/-S.E.M. preovulatory plasma LH peak was 10.11+/-0.35 ng/ml with the values ranging from 8.75 to 11.51 ng/ml in individual yaks. The mean+/-S.E.M. duration of the LH surge was 7.25+/-0.55 h with a range of 6-10 h. Onset of LH surge (mean+/-S.E.M.) occurred 3.0+/-0.65 h after the onset of estrus. Mean plasma progesterone stayed low (<0.25 ng/ml) during the entire duration of sampling. Ovulation occurred 30.5+/-0.82 h (range, 28-34 h) after the onset of estrus and 20.25+/-1.03 h after the end of LH surge. The occurrence of the LH peaks within a narrow time frame of 4-8h post estrus onset in yaks could have contributed to the animals ovulating within a narrow time interval.

Profiles of plasma progesterone before and at the onset of puberty in yak heifers

Sixty yak heifers were closely observed to determine their ages at sexual maturity and 15 heifers were used to analyze the plasma progesterone by radioimmunoassay before and at the attainment of puberty. Yak heifers began cyclic activity at a mean (+/-S.D.) age of 33 +/- 6.7 months, but with wild variation. Fifteen heifers used to determine plasma progesterone were divided into three age groups: group I (10-14 months, n = 5), group II (20-24 months, n = 5) and group III (30-36 months, n = 5). Group I yak heifers were found to have two progesterone profiles: inactive ovary (IO) profile and low progesterone short cycle (LPSC) profile; group II had three profiles: IO profile, LPSC profile and LPNC (low progesterone normal cycle) profile; group III had three profiles: LPSC profile, LPNC profile and a normal oestrous cycle profile. It would be concluded from the present study that one or two, or even more brief rise(s) in circulating progesterone presented in yak heifers. These rises, however, were not followed by a normal luteal phase except two yak heifers that came to heat. This characteristics found in yak was obviously different with that reported in other animals.

Reproductive patterns of the yak. II. Progesterone and oestradiol-17 beta levels in plasma and milk just before the breeding season; also during normal and short oestrous cycles

Concentrations of oestradiol-17 beta and progesterone in plasma and milk of the yak before the first oestrus at the beginning of the breeding season and during normal and short oestrus cycles were measured by radioimmunoassay (RIA). Both oestradiol-17 beta and progesterone were maintained at stable, low levels before the first oestrus. However, oestradiol-17 beta showed a peak, similar to that on the day of oestrus (P > 0.05), about 6 days before the first oestrus and a small progesterone peak was seen 2 days before it. There were three peaks of oestradiol-17 beta in plasma and milk on the day of oestrus [24.57 (SD 11.35) pg ml-1 and 97.84 (SD 26.02) pg ml-1 respectively], and on days 5 and 14 of the cycle. Progesterone levels in plasma and milk were low during oestrus but peaks were seen on day 15 [5.23 (SD 1.02) ng ml-1 and 20.44 (SD 3.64) ng ml-1 respectively]. The patterns of oestradiol-17 beta and progesterone during the short cycle were similar to those during the normal cycle (P > 0.05), but the values were lower. The levels of oestradiol-17 beta and progesterone in milk were about 3 or 4 times higher than those in plasma before the first oestrus, and about 4 or 5 times higher during the normal cycle.