Equivalence of variance components between standard and recursive genetic models using LDL' transformations.

BACKGROUND: Recursive models are a category of structural equation models that propose a causal relationship between traits. These models are more parameterized than multiple trait models, and they require imposing restrictions on the parameter space to ensure statistical identification. Nevertheless, in certain situations, the likelihood of recursive models and multiple trait models are equivalent. Consequently, the estimates of variance components derived from the multiple trait mixed model can be converted into estimates under several recursive models through LDL' or block-LDL' transformations. RESULTS: The procedure was employed on a dataset comprising five traits (birth weight-BW, weight at 90 days-W90, weight at 210 days-W210, cold carcass weight-CCW and conformation-CON) from the Pirenaica beef cattle breed. These phenotypic records were unequally distributed among 149,029 individuals and had a high percentage of missing data. The pedigree used consisted of 343,753 individuals. A Bayesian approach involving a multiple-trait mixed model was applied using a Gibbs sampler. The variance components obtained at each iteration of the Gibbs sampler were subsequently used to estimate the variance components within three distinct recursive models. CONCLUSIONS: The LDL' or block-LDL' transformations applied to the variance component estimates achieved from a multiple trait mixed model enabled inference across multiple sets of recursive models, with the sole prerequisite of being likelihood equivalent. Furthermore, the aforementioned transformations simplify the handling of missing data when conducting inference within the realm of recursive models.

Conservation and upgradation of indigenous Cholistani cattle breed of Pakistan: a pathway to sustainable livestock.

Pakistan is endowed with many established indigenous zebu Bos indicus type (humped) cattle breeds including Sahiwal, Red Sindhi, Bhagnari and Cholistani. Amongst these indigenous cattle breeds, Sahiwal and Red Sindhi have extensively been navigated and hence these two are acclaimed as internationally recognized breeds. However, research work on Cholistani cattle breed actually initiated in 2010 and has attained a steady pace. This breed was a new entrant in Livestock Census of Pakistan since 2006. Cholistani is a hardy, tick-resistant, adaptable cattle breed being reared under pastoral nomadism of the Cholistan desert, Pakistan. The present narrative review is the first of its kind intended to sum-up all the research work conducted about this indigenous cattle breed, and to put forth research gaps for this formerly neglected cattle breed. The review discusses the research work conducted on Cholistani cattle breed under five major research subjects/domains i.e. production attributes, theriogenology-related attributes, hematochemical attributes, disease, epidemiologic and therapeutic attributes, and genetic attributes. Future horizon for research avenues has also been given. It is the dire need of time that specific breed-oriented conservation and propagation programs may be initiated in the country so that sustained livestock and enhance socioeconomic profiling of rural communities may be attained.

Nuclear sequences of mitochondrial origin in domestic yak.

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle.

BACKGROUND: Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle. RESULTS: We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits. CONCLUSIONS: This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Lifetime milk production of Holstein cattle in the humid tropics compared to Holstein-Gyr and Holstein-Brahman crosses.

Crossbred cattle are commonly used for milk production in the tropics, combining the potential benefits of pure breeds with the heterosis effects of the offspring. However, no comprehensive assessment of lifetime productivity for crossbred versus purebred cattle in low-altitude tropical environments has been carried out. The present study compares the lifetime productivity of purebred Holstein (HO, n = 17,269), Gyr (GY4, n = 435), and Brahman (BR4, n = 622) with crossbreds Gyr × Holstein (GY × HO, n = 5521) and Brahman×Holstein (BR × HO, n = 5429) cows from dairy farms located in low and medium altitude tropical regions in Costa Rica. The production traits of interest were age at first calving (AFC), days open (DO), milk production per lactation (TMP), lactation length (LLEN), age at culling (ACUL), and number of lactations (NLAC). Estimates of heterosis were also calculated. The AFC for GY × HO crosses (33-34 months) was not significantly different (p > .05) from HO (33.8 months). For BR × HO crosses, a significant (p < .05) decrease in AFC (BR3HO1 35.6 months, BR2HO2 34.5 months, and BR1H03 33.3 months) was observed as the fraction of HO breed increased. Estimates of heterosis for AFC were favourable for both crosses, of a magnitude close to 3%. The DO for F1 crosses (GY2HO2 94 days; BR2HO2 96 days) was significantly (p < .05) lower than HO (123 days). Estimates of heterosis for DO were also favourable and above 15% for both crosses. The TMP and LLEN were higher for HO (TMP = 5003 kg; LLEN = 324 days) compared with GY × HO (TMP = 4428 to 4773 kg; LLEN = 298 to 312 days) and BR × HO (TMP = 3950 to 4761 kg; LLEN = 273 to 313 days) crosses. Heterosis for TMP was favourable but low for both crosses, with a magnitude below 3.0%. The NLAC for HO (4.6 lactations) was significantly (p < .05) lower than F1 (GY2HO2, 5.8 lactations; BR2HO2, 5.4 lactations). Heterosis for NLAC was above 6.0% for both crosses. Overall, estimates of lifetime income over feed costs per cow on average were USD 2637 (30.3%) and USD 734 (8.4%) higher in F1 GY × HO and BR × HO, respectively, compared to HO. In conclusion, crossbred animals, specifically those with Gyr and Brahman genetics, extend the productive lifespan, increasing economic returns.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Deciphering genetic diversity in conserved cattle bulls to achieve sustainable development goals.

The primary objective of Sustainable Development Goal target 2.5 established by the United Nations is to ensure the preservation of genetic diversity in domesticated animals. The ICAR-National Bureau of Animal Genetic Resources in India has been actively engaged in the conservation of cattle and buffalo bull semen for long-term storage. This present study aimed to assess the genetic diversity present in the conserved cattle bull semen, which would aid in determining the most suitable strategy for future conservation management. A total of 192 bull semen belonging to 19 cattle breeds were selected to evaluate genetic diversity using 17 pairs of FAO recommended microsatellite primers. Total 267 alleles were detected across all the samples which indicates substantial amount of allelic variation is being maintained in conserved bulls. Further, all cattle bulls semen conserved showed higher observed heterozygosity than expected heterozygosity which indicates excess genetic diversity in all the populations. The FST, F IT and FIS value across the loci and population is 0.146 ± 0.009, 0.054 ± 0.038, and - 0.105 ± 0.035, respectively, which suggests lack of inbreeding in conserved cattle bull semen. This study has established genetic diversity in conserved cattle semen samples to achieve sustainable development goals. In addition, it provides compelling evidence that the current approach for conserving cattle bull semen is heading in the correct direction.

The effect of miR-23b-3p on regulating GH by targeting POU1F1 in Yanbian yellow cattle.

This study aimed to evaluate the effect of miR-23b-3p on growth hormone (GH) in pituitary cells of Yanbian yellow cattle. The mRNA and protein levels of GH and miR-23b-3p target genes were measured by real time fluorescence quantitative PCR (qPCR) and Western blot, respectively. The target relationship of miR-23b-3p was validated by double luciferase reporter gene system. The results showed that GH mRNA and protein levels in pituitary cells of Yanbian yellow cattle were significantly lower in the miR-23b-3p-mi group than in the NC group (P<0.01), while GH mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.05). The result of bioinformatics analysis and double luciferase reporter gene system validation proved that miR-23b-3p targeted 3'UTR of pituitary specific transcription factor 1 (POU1F1). POU1F1 mRNA and protein levels were lower miR-23b-3p-mi group than in the NC group (P<0.01), while POU1F1 mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.01). These results demonstrated that miR-23b-3p could regulate GH expression in pituitary cells by regulating POU1F1 gene.

Nutritional and productive performance of purebred Nellore heifers and crossed with Brangus and Braford finished in a feedlot system.

This study evaluated the nutritional and productive performance of Nellore purebred heifers and crossbred Brangus x Nellore (BGNE) and Braford x Nellore (BFNE) in a feedlot system. Thirty heifers (10 of each genetic group) with an average age of 18 months and an initial body weight of 261 kg were used. The experiment was structured and conducted according to a completely randomized design, with three treatments. Heifers received two diets (60 days each) during the experimental period. The experiment lasted 120 days with four experimental periods. Nellore heifers had a lower intake than crossbred heifers (P <0.05). There were no differences between BGNE and BFNE heifers, which had higher final body weight, average daily gain, feed efficiency, hot carcass weight and carcass length than NE heifers. Crossed heifers presented better fat cover than NE heifers. However, NE heifers had higher carcass dressing Despite presenting lower carcass yields than Nellore heifers, crossed heifers are more efficient and have higher performance and better fat cover on the carcass than purebred Nellore heifers. Crossbreeding synthetic breeds, such as Brangus and Braford breeds, with the Nellore breed is an effective way to increase the productivity and efficiency of feedlot heifers in tropical regions.

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip.

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1-5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens).

BACKGROUND: There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity. RESULTS: A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks. CONCLUSIONS: Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Genetic parameters for milk yield and reproductive traits in Honduran Holstein cattle.

The aim of this study was to estimate the (co)variance components and genetic parameters for milk yield adjusted to 305d (MY305), calving-to-conception interval (CCI), number of services per conception (NSC) and calving interval (CI) of Honduran Holstein cows, by fitting a bivariate animal model using Maximum Restricted Likelihood procedures. Model included the fixed effects of calving number, the contemporary calving group (farm-season-year of calving and the cow age as covariate). The estimated means and standard deviations for MY, CCI, NSC and CI were, 5098.60 ± 1564.32 kg, 168.27 ± 104.71 days, 2.46 ± 1.69 services, and 448.73 ± 109.16 days, respectively; and their estimated heritabilities were 0.21 ± 0.05, 0.03 ± 0.028, 0.02 ± 0.024 and 0.06 ± 0.04, respectively. The genetic correlations between MY305 and CCI, NSC and CI were positive and antagonist, with values of 0.64 ± 0.52, 0.99 ± 0.56, and 0.32 ± 0.24 respectively. Even though moderate to low heritability was estimated for MY305, systematic selection for milk yield, with a reduction in reproductive efficiency, if considered as the only selection criterion is important to be considered. By including reproductive traits and considering permanent environment effects into the breeding program, might yield a slow, but constant and permanent improvement over time.

Double-Tube Multiplex TaqMan Real-Time PCR for the Detection of Eight Animal-Derived Dairy Ingredients.

Milk and dairy products represent important sources of nutrition in our daily lives. The identification of species within dairy products holds importance for monitoring food adulteration and ensuring traceability. This study presented a method that integrated double-tube and duplex real-time polymerase chain reaction (PCR) with multiplex TaqMan probes to enable the high-throughput detection of animal-derived ingredients in milk and dairy products. The detection system utilized one pair of universal primers, two pairs of specific primers, and eight animal-derived specific probes for cow, buffalo, goat, sheep, camel, yak, horse, and donkey. These components were optimized within a double-tube and four-probe PCR multiplex system. The developed double-tube detection system could simultaneously identify the above eight targets with a detection limit of 10-0.1 pg/µL. Validation using simulated adulterated milk samples demonstrated a detection limit of 0.1%. The primary advantage of this method lies in the simplification of the multiplex quantitative real-time PCR (qPCR) system through the use of universal primers. This method provides an efficient approach for detecting ingredients in dairy products, providing powerful technical support for market supervision.

Genome-wide association study identifies genomic regions associated with key reproductive traits in Korean Hanwoo cows.

BACKGROUND: Conducting genome-wide association studies (GWAS) for reproductive traits in Hanwoo cattle, including age at first calving (AFC), calving interval (CI), gestation length (GL), and number of artificial inseminations per conception (NAIPC), is of paramount significance. These analyses provided a thorough exploration of the genetic basis of these traits, facilitating the identification of key markers for targeted trait improvement. Breeders can optimize their selection strategies, leading to more efficient and sustainable breeding programs, by incorporating genetic insights. This impact extends beyond individual traits and contributes to the overall productivity and profitability of the Hanwoo beef cattle industry. Ultimately, GWAS is essential in ensuring the long-term genetic resilience and adaptability of Hanwoo cattle populations. The primary goal of this study was to identify significant single nucleotide polymorphisms (SNPs) or quantitative trait loci (QTLs) associated with the studied reproductive traits and subsequently map the underlying genes that hold promise for trait improvement. RESULTS: A genome-wide association study of reproductive traits identified 68 significant single nucleotide polymorphisms (SNPs) distributed across 29 Bos taurus autosomes (BTA). Among them, BTA14 exhibited the highest number of identified SNPs (25), whereas BTA6, BTA7, BTA8, BTA10, BTA13, BTA17, and BTA20 exhibited 8, 5, 5, 3, 8, 2, and 12 significant SNPs, respectively. Annotation of candidate genes within a 500 kb region surrounding the significant SNPs led to the identification of ten candidate genes relevant to age at first calving. These genes were: FANCG, UNC13B, TESK1, TLN1, and CREB3 on BTA8; FAM110B, UBXN2B, SDCBP, and TOX on BTA14; and MAP3K1 on BTA20. Additionally, APBA3, TCF12, and ZFR2, located on BTA7 and BTA10, were associated with the calving interval; PAX1, SGCD, and HAND1, located on BTA7 and BTA13, were linked to gestation length; and RBM47, UBE2K, and GPX8, located on BTA6 and BTA20, were linked to the number of artificial inseminations per conception in Hanwoo cows. CONCLUSIONS: The findings of this study enhance our knowledge of the genetic factors that influence reproductive traits in Hanwoo cattle populations and provide a foundation for future breeding strategies focused on improving desirable traits in beef cattle. This research offers new evidence and insights into the genetic variants and genome regions associated with reproductive traits and contributes valuable information to guide future efforts in cattle breeding.

Spatially resolved transcriptomic profiling of placental development in dairy cow.

The placenta plays a crucial role in successful mammalian reproduction. Ruminant animals possess a semi-invasive placenta characterized by a highly vascularized structure formed by maternal endometrial caruncles and fetal placental cotyledons, essential for full-term fetal development. The cow placenta harbors at least two trophoblast cell populations: uninucleate (UNC) and binucleate (BNC) cells. However, the limited capacity to elucidate the transcriptomic dynamics of the placental natural environment has resulted in a poor understanding of both the molecular and cellular interactions between trophoblast cells and niches, and the molecular mechanisms governing trophoblast differentiation and functionalization. To fill this knowledge gap, we employed Stereo-seq to map spatial gene expression patterns at near single-cell resolution in the cow placenta at 90 and 130 days of gestation, attaining high-resolution, spatially resolved gene expression profiles. Based on clustering and cell marker gene expression analyses, key transcription factors, including YBX1 and NPAS2, were shown to regulate the heterogeneity of trophoblast cell subpopulations. Cell communication and trajectory analysis provided a framework for understanding cell-cell interactions and the differentiation of trophoblasts into BNCs in the placental microenvironment. Differential analysis of cell trajectories identified a set of genes involved in regulation of trophoblast differentiation. Additionally, spatial modules and co-variant genes that help shape specific tissue structures were identified. Together, these findings provide foundational insights into important biological pathways critical to the placental development and function in cows.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak.

BACKGROUND: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY). RESULTS: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak. CONCLUSIONS: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Gene expression supports a single origin of horns and antlers in hoofed mammals.

Horns, antlers, and other bony cranial appendages of even-toed hoofed mammals (ruminant artiodactyls) challenge traditional morphological homology assessments. Cranial appendages all share a permanent bone portion with family-specific integument coverings, but homology determination depends on whether the integument covering is an essential component or a secondary elaboration of each structure. To enhance morphological homology assessments, we tested whether juvenile cattle horn bud transcriptomes share homologous gene expression patterns with deer antlers relative to pig outgroup tissues, treating the integument covering as a secondary elaboration. We uncovered differentially expressed genes that support horn and antler homology, potentially distinguish them from non-cranial-appendage bone and other tissues, and highlight the importance of phylogenetic outgroups in homology assessments. Furthermore, we found differentially expressed genes that could support a shared cranial neural crest origin for horns and antlers and expression patterns that refine our understanding of the timing of horn and antler differentiation.

Development of PCR-based genetic test for detection of TTF1 mutation causing abortion in Holstein Friesian cattle.

A stop-gain mutation (rs715966442; BTA11: 1,02,463,944 nucleotide position) in transcription termination factor, RNA polymerase I (TTF1) gene causes abortion in Holstein Friesian (HF) cattle. A PCR-restriction fragment length polymorphism (PCR-RFLP)-based genetic test has been developed and validated to screen the TTF1 mutation locus in HF cattle. The mutation locus was screened in 80 HF and HF crossbreds using the protocol, which revealed two animals as carriers of the mutant TTF1 allele. The test employed is cost-effective, rapid and precise and can be utilized as an effective tool for the screening of TTF1 mutation carriers in HF cattle population.

Molecular characterization and expression profile of the ALDH1A1 gene and its functions in yak luteal cells.

The ALDH1A1 gene encodes a cytoplasmic member of the aldehyde dehydrogenase 1 family, which plays an important role in regulating animal reproductive performance, including estrus cycle and embryonic development. The aim of this study was to characterize ALDH1A1 activity in ovaries of 3-5 year-old yaks and to determine its effects on cell proliferation, apoptosis, and progesterone secretion in luteal cells (LCs). The coding sequence (CDS) of the ALDH1A1 gene was cloned by reverse transcription-PCR and immunohistochemical analysis was used to confirm localization of the ALDH1A1 protein in the ovary. To assess the activity of ALDH1A1 in regulating progesterone secretion, si-ALDH1A1 was transfected into LCs in vitro and progesterone levels in LC supernatants were measured by ELISA. The interference efficiency was assessed by real-time quantitative PCR (RT-qPCR) and immunofluorescence staining, and cell proliferation and apoptosis were evaluated by EdU and TUNEL staining, respectively. The cloned ALDH1A1 sequence contained 1462 bp, encoding 487 amino acids. Immunohistochemical analysis showed that ALDH1A1 protein expression, which was significantly higher in LCs, was mainly found in antral follicles and the corpus luteum (CL). The expression of ALDH1A1 mRNA in LCs was effectively inhibited by si-ALDH1A1transfection, and progesterone secretion was markedly decreased along with the significant down-regulation of progesterone pathway-related genes, STAR, CYP11A1, CYP19A1, CYP17A1, 3ß-HSD, and HSD17B1. Knockdown of ALDH1A1 mRNA expression decreased cell proliferation and increased apoptosis in LCs. The mRNA expression of the proliferation-related genes, PCNA, CCND1, CCNB1 and CDC25A, was significantly down-regulated, while expression of the apoptosis-promoting CASP3 gene was significantly increased. In summary, we characterized the yak ALDH1A1 gene and revealed that ALDH1A1 knockdown promoted apoptosis, repressed cell proliferation, and decreased progesterone secretion by yak LCs, potentially by regulating the mRNA expression of genes related to proliferation, apoptosis, and progesterone synthesis and secretion.

Evaluation of non-additive genetic effects on carcass and meat quality traits in Korean Hanwoo cattle using genomic models.

The traditional genetic evaluation methods generally consider additive genetic effects only and often ignore non-additive (dominance and epistasis) effects that may have contributed to genetic variation of complex traits of livestock species. The available dense single nucleotide polymorphisms (SNPs) panels offer to investigate the potential benefits of including non-additive genetic effects in the genomic evaluation models. Data from 16 971 genotyped (Illumina Bovine 50 K SNP chip) Korean Hanwoo cattle were used to estimate genetic variance components and prediction accuracy of genomic breeding values (GEBVs) for four carcass and meat quality traits: carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT) and marbling score (MS). Five different genetic models were evaluated through including additive, dominance and epistatic interactions (additive by additive, A × A; additive by dominance, A × D and dominance by dominance, D × D) successively in the models. The estimates of additive genetic variances and narrow sense heritabilities (ha2) were found similar across the evaluated models and traits except when additive interaction (A × A) was included. The dominance variance estimates relative to phenotypic variance ranged from 1.7-3.4% for CWT and MS traits, whereas, they were close to zero for EMA and BFT traits. The magnitude of A × A epistatic heritability (haa2) ranged between 14.8 and 27.7% in all traits. However, heritability estimates for A × D and D × D epistatic interactions (had2 and hdd2) were quite low compared to haa2 and were contributed only 0.0-9.7% of the total phenotypic variation. In general, broad sense heritability (hG2) estimates were almost twice (ranging between 0.54 and 0.68) the ha2 for all of the investigated traits. The inclusion of dominance effects did not improve the prediction accuracy of GEBV but improved 2.0-3.0% when epistatic effects were included in the model. More importantly, rank correlation revealed that partitioning of variance components considering dominance and epistatic effects in the model would enable to re-rank of top animals with better prediction of GEBV. The present result suggests that dominance and epistatic effects could be included in the genomic evaluation model for better estimates of variance components and more accurate prediction of GEBV for carcass and meat quality traits in Korean Hanwoo cattle.