Genetic variation in PLAG1 associates with early life body weight and peripubertal weight and growth in Bos taurus

Genome-wide association study of age at puberty and its (co)variances with fertility and stature in growing and lactating Holstein-Friesian dairy cattle

Reproductive performance is a key determinant of cow longevity in a pasture-based, seasonal dairy system. Unfortunately, direct fertility phenotypes such as intercalving interval or pregnancy rate tend to have low heritabilities and occur relatively late in an animal's life. In contrast, age at puberty (AGEP) is a moderately heritable, early-in-life trait that may be estimated using an animal's age at first measured elevation in blood plasma progesterone (AGEP4) concentrations. Understanding the genetic architecture of AGEP4 in addition to genetic relationships between AGEP4 and fertility traits in lactating cows is important, as is its relationship with body size in the growing animal. Thus, the objectives of this research were 3-fold. First, to estimate the genetic and phenotypic (co)variances between AGEP4 and subsequent fertility during first and second lactations. Second, to quantify the associations between AGEP4 and height, length, and BW measured when animals were approximately 11 mo old (standard deviation = 0.5). Third, to identify genomic regions that are likely to be associated with variation in AGEP4. We measured AGEP4, height, length, and BW in approximately 5,000 Holstein-Friesian or Holstein-Friesian × Jersey crossbred yearling heifers across 54 pasture-based herds managed in seasonal calving farm systems. We also obtained calving rate (CR42, success or failure to calve within the first 42 d of the seasonal calving period), breeding rate (PB21, success or failure to be presented for breeding within the first 21 d of the seasonal breeding period) and pregnancy rate (PR42, success or failure to become pregnant within the first 42 d of the seasonal breeding period) phenotypes from their first and second lactations. The animals were genotyped using the Weatherby's Versa 50K SNP array (Illumina, San Diego, CA). The estimated heritabilities of AGEP4, height, length, and BW were 0.34 (90% credibility interval [CRI]: 0.30, 0.37), 0.28 (90% CRI: 0.25, 0.31), 0.21 (90% CRI: 0.18, 0.23), and 0.33 (90% CRI: 0.30, 0.36), respectively. In contrast, the heritabilities of CR42, PB21 and PR42 were all <0.05 in both first and second lactations. The genetic correlations between AGEP4 and these fertility traits were generally moderate, ranging from 0.11 to 0.60, whereas genetic correlations between AGEP4 and yearling body-conformation traits ranged from 0.02 to 0.28. Our GWAS highlighted a genomic window on chromosome 5 that was strongly associated with variation in AGEP4. We also identified 4 regions, located on chromosomes 14, 6, 1, and 11 (in order of decreasing importance), that exhibited suggestive associations with AGEP4. Our results show that AGEP4 is a reasonable predictor of estimated breeding values for fertility traits in lactating cows. Although the GWAS provided insights into genetic mechanisms underpinning AGEP4, further work is required to test genomic predictions of fertility that use this information.

Refining the genomic profiles of North African sheep breeds through meta-analysis of worldwide genomic SNP data

Introduction

The development of reproducible tools for the rapid genotyping of thousands of genetic markers (SNPs) has promoted cross border collaboration in the study of sheep genetic diversity on a global scale.

Methods

In this study, we collected a comprehensive dataset of 239 African and Eurasian sheep breeds genotyped at 37,638 filtered SNP markers, with the aim of understanding the genetic structure of 22 North African (NA) sheep breeds within a global context.

Results and discussion

We revealed asubstantial enrichment of the gene pool between the north and south shores of the Mediterranean Sea, which corroborates the importance of the maritime route in the history of livestock. The genetic structure of North African breeds mirrors the differential composition of genetic backgrounds following the breed history. Indeed, Maghrebin sheep stocks constitute a geographically and historically coherent unit with any breed-level genetic distinctness among them due to considerable gene flow. We detected a broad east-west pattern describing the most important trend in NA fat-tailed populations, exhibited by the genetic closeness of Egyptian and Libyan fat-tailed sheep to Middle Eastern breeds rather than Maghrebin ones. A Bayesian FST scan analysis revealed a set of genes with potentially key adaptive roles in lipid metabolism (BMP2, PDGFD VEGFA, TBX15, and WARS2), coat pigmentation (SOX10, PICK1, PDGFRA, MC1R, and MTIF) and horn morphology RXFP2) in Tunisian sheep. The local ancestry method detected a Merino signature in Tunisian Noire de Thibar sheep near the SULF1gene introgressed by Merino's European breeds. This study will contribute to the general picture of worldwide sheep genetic diversity.

Multi-tissue transcriptome analysis to identify candidate genes associated with weight regulation in Hanwoo cattle

While genetic markers related to meat production traits have been identified in many other cattle breeds, research on weight in Hanwoo cattle (Korean native cattle) is still insufficient. In this study, we performed expression quantitative trait loci (eQTL) analysis and differential gene expression analysis to detect candidate genes influencing the weight characteristics of 32 castrated Hanwoo cattle across 22 tissues and, we identified variants that affect gene expression levels. In total, we identified a total of 3,298 differentially expressed genes, among which we discovered key genes such as UBD, RGS2, FASN, and SCD that have functions related to adipogenesis, body weight, obesity, and lipid metabolism. Gene-set enrichment analysis revealed that candidate genes in adipose tissue are involved in metabolic pathways linked to obesity-related traits, adipose metabolism, and lipid metabolism. Additionally, we found that decreased expression of TRIM31 contributes to weight gain which can be explained by the associated candidate cis-eQTL genotypes for TRIM31 and their effect on differential gene expression between the lower and higher weight groups. Our findings revealed candidate genes associated with the weight of Hanwoo cattle and perhaps can provide comprehensive insights into the association of weight with various tissues beyond adipose tissue and muscle, indicating the potential for expanding the focus of livestock trait research.

Pleomorphic adenoma gene1 in reproduction and implication for embryonic survival in cattle: a review

The pleomorphic adenoma gene1 (PLAG1) encodes a DNA-binding, C2H2 zinc-finger protein which acts as a transcription factor that regulates the expression of diverse genes across different organs and tissues; hence, the name pleomorphic. Rearrangements of the PLAG1 gene, and/or overexpression, are associated with benign tumors and cancers in a variety of tissues. This is best described for pleomorphic adenoma of the salivary glands in humans. The most notable expression of PLAG1 occurs during embryonic and fetal development, with lesser expression after birth. Evidence has accumulated of a role for PLAG1 protein in normal early embryonic development and placentation in mammals. PLAG1 protein influences the expression of the ike growth factor 2 (IGF2) gene and production of IGF2 protein. IGF2 is an important mitogen in ovarian follicles/oocytes, embryos, and fetuses. The PLAG1-IGF2 axis, therefore, provides one pathway whereby PLAG1 protein can influence embryonic survival and pregnancy. PLAG1 also influences over 1,000 other genes in embryos including those associated with ribosomal assembly and proteins. Brahman (Bos indicus) heifers homozygous for the PLAG1 variant, rs109815800 (G > T), show greater fertility than contemporary heifers with either one, or no copy, of the variant. Greater fertility in heifers homozygous for rs109815800 could be the result of early puberty and/or greater embryonic survival. The present review first looks at the broader roles of the PLAG1 gene and PLAG1 protein and then focuses on the emerging role of PLAG1/PLAG1 in embryonic development and pregnancy. A deeper understanding of factors which influence embryonic development is required for the next transformational increase in embryonic survival and successful pregnancy for both in vivo and in vitro derived embryos in cattle.

Genome-wide association study identifies functional genomic variants associated with young stock survival in Nordic Red Dairy Cattle

Identifying quantitative trait loci (QTL) associated with calf survival is essential for both reducing economic loss in cattle industry and understanding the genetic basis of the trait. To identify mutations and genes underlying young stock survival (YSS), we performed GWAS using de-regressed estimated breeding values of a YSS index and its component traits defined by sex and age in 3,077 Nordic Red Dairy Cattle (RDC) bulls and 2 stillbirth traits (first lactation and later lactations) in 5,141 RDC bulls. Two associated QTL regions on Bos taurus autosome (BTA) 4 and 6 were identified for the YSS index. The results of 4 YSS component traits indicate that same QTL regions were associated with bull and heifer calf mortality, but the effects were different over the growing period and suggested an additional QTL on BTA23. The GWAS on stillbirth identified 3 additional QTL regions on BTA5, 14, and 24 compared with YSS and its component traits. The conditional test of BTA6 showed at least 2 closely located QTL segregating for YSS component traits and stillbirth. We found 2 independent QTL for stillbirth on BTA23. The post-GWAS revealed LCORL, PPM1K, SSP1, MED28, and LAP3 are putative causal genes on BTA6, and a frame shift variant within LCORL, BTA6:37401770 (rs384548488) could be the putative causal variant. On BTA4, the GRB10 gene is the putative causal gene and BTA4:5296018 is the putative causal variant. In addition, NDUFA9 and FGF23 on BTA5, LYN on BTA14, and KCNK5 on BTA23 are putative causal genes for QTL for stillbirth. The gene analysis also proposed several candidate genes. Our findings shed new light on the candidate genes affecting calf survival, and the knowledge could be utilized to reduce calf mortality and thereby enhance welfare of dairy cattle.

Genome-wide association study for growth traits in Blanco Orejinegro and Romosinuano cattle

Growth traits are economically important characteristics for the genetic improvement of local cattle breeds. Genome-wide association studies (GWAS) provide valuable information to enhance the understanding on the genetics of complex traits. The aim of this study was to perform a GWAS to identify genomic regions and genes associated to birth weight, weaning weight adjusted for 240 days, 16 months, and 24 months weight in Romosinuano (ROMO) and Blanco Orejinegro (BON) cattle. A single-step genomic-BLUP was implemented using 596 BON and 569 ROMO individuals that were genotyped with an Illumina BovineSNP50 BeadChip. There were 25 regions of interest identified on different chromosomes, with few of them simultaneously associated with two or more growth traits and some were common to both breeds. The gene mapping allowed to find 173 annotations on these regions, from which 49 represent potential candidate genes with known growth-related functions in cattle and other species. Among the regions that were associated with several growth traits, that at 24 - 27 MB of BTA14, has important candidate genes such as LYPLA1, XKR4, TMEM68 and PLAG1. Another region of interest at 0.40-0.77 Mb of BTA23 was identified in both breeds, containing KHDRBS2 as a potential candidate gene influencing body weight. Future studies targeting these regions could provide more knowledge to uncover the genetic architecture underlying growth traits in BON and ROMO cattle. The genomic regions and genes identified in this study could be used to improve the prediction of genetic merit for growth traits in these creole cattle breeds.

Bovine muscle satellite cells in calves and cattle: A comparative study of cellular and genetic characteristics for cultivated meat production

This study compared the cellular and genetic characteristics of bovine skeletal muscle satellite cells (SMSCs) from Hanwoo (a Korean native cattle breed), including calves and mature cattle. SMSCs were isolated using magnetic-activated cell sorting (MACS) from tissue samples of six Hanwoo (three calves and three mature cattle) using the CD29 antibody. Calves' SMSCs exhibited significantly faster growth rates than did those from cattle (P < 0.01), with a doubling time of 2.43 days. Genetic analysis revealed higher MyoD and Pax7 expression in SMSCs from calves during proliferation than in those from mature cattle (P < 0.001). However, FASN and PLAG1 expression levels were higher in mature cattle than in calves during both proliferation and differentiation (P < 0.001). These findings highlight the need for strategies to improve bovine muscle cell growth to produce competitive cultivated meat at a competitive price.

Sequence level genome-wide associations for bull production and fertility traits in tropically adapted bulls

BACKGROUND

The genetics of male fertility is complex and not fully understood. Male subfertility can adversely affect the economics of livestock production. For example, inadvertently mating bulls with poor fertility can result in reduced annual liveweight production and suboptimal husbandry management. Fertility traits, such as scrotal circumference and semen quality are commonly used to select bulls before mating and can be targeted in genomic studies. In this study, we conducted genome-wide association analyses using sequence-level data targeting seven bull production and fertility traits measured in a multi-breed population of 6,422 tropically adapted bulls. The beef bull production and fertility traits included body weight (Weight), body condition score (CS), scrotal circumference (SC), sheath score (Sheath), percentage of normal spermatozoa (PNS), percentage of spermatozoa with mid-piece abnormalities (MP) and percentage of spermatozoa with proximal droplets (PD).

RESULTS

After quality control, 13,398,171 polymorphisms were tested for their associations with each trait in a mixed-model approach, fitting a multi-breed genomic relationship matrix. A Bonferroni genome-wide significance threshold of 5 × 10^- 8 was imposed. This effort led to identifying genetic variants and candidate genes underpinning bull fertility and production traits. Genetic variants in Bos taurus autosome (BTA) 5 were associated with SC, Sheath, PNS, PD and MP. Whereas chromosome X was significant for SC, PNS, and PD. The traits we studied are highly polygenic and had significant results across the genome (BTA 1, 2, 4, 6, 7, 8, 11, 12, 14, 16, 18, 19, 23, 28, and 29). We also highlighted potential high-impact variants and candidate genes associated with Scrotal Circumference (SC) and Sheath Score (Sheath), which warrants further investigation in future studies.

CONCLUSION

The work presented here is a step closer to identifying molecular mechanisms that underpin bull fertility and production. Our work also emphasises the importance of including the X chromosome in genomic analyses. Future research aims to investigate potential causative variants and genes in downstream analyses.

A look under the hood of genomic-estimated breed compositions for brangus cattle: What have we learned?

The Brangus cattle were developed to utilize the superior traits of Angus and Brahman cattle. Their genetic compositions are expected to be stabilized at 3/8 Brahman and 5/8 Angus. Previous studies have shown more than expected Angus lineage with Brangus cattle, and the reasons are yet to be investigated. In this study, we revisited the breed compositions for 3,605 Brangus cattle from three perspectives: genome-wise (GBC), per chromosomes (CBC), and per chromosome segments (SBC). The former (GBC) depicted an overall picture of the “mosaic” genome of the Brangus attributable to their ancestors, whereas the latter two criteria (CBC and SBC) corresponded to local ancestral contributions. The average GBC for the 3,605 Brangus cattle were 70.2% Angus and 29.8% Brahman. The K-means clustering supported the postulation of the mixture of 1/2 Ultrablack (UB) animals in Brangus. For the non-UB Brangus animals, the average GBC were estimated to be 67.4% Angus and 32.6% Brahman. The 95% confidence intervals of their overall GBC were 60.4%–73.5% Angus and 26.5%–39.6% Brahman. Possibly, genetic selection and drifting have resulted in an approximately 5% average deviation toward Angus lineage. The estimated ancestral contributions by chromosomes were heavily distributed toward Angus, with 27 chromosomes having an average Angus CBC greater than 62.5% but only two chromosomes (5 and 20) having Brahman CBC greater than 37.5%. The chromosomal regions with high Angus breed proportions were prevalent, tending to form larger blocks on most chromosomes. In contrast, chromosome segments with high Brahman breed proportion were relatively few and isolated, presenting only on seven chromosomes. Hence, genomic hitchhiking effects were strong where Angus favorable alleles resided but weak where Brahman favorable alleles were present. The functions of genes identified in the chromosomal regions with high (≥75%) Angus compositions were diverse yet may were related to growth and body development. In contrast, the genes identified in the regions with high (≥37.5%) Brahman compositions were primarily responsible for disease resistance. In conclusion, we have addressed the questions concerning the Brangus genetic make-ups. The results can help form a dynamic picture of the Brangus breed formation and the genomic reshaping.

PLAG1 g.8795C>T Mutation Regulates Early Body Weight in Hu Sheep by Weakening miR-139 Binding

Sheep birth and weaning weights indicate their growth and survival. Thus, identifying molecular genetic markers for early body weight is important in sheep breeding. Pleomorphic adenoma gene 1 (PLAG1) is important for regulating birth weight and body length in mammals; however, its relationship with sheep body weight remains unknown. Here, the 3'-untranslated region (3'-UTR) of the Hu sheep PLAG1 gene was cloned, single nucleotide polymorphisms (SNPs) were screened, genotype-early body weight relationships were analyzed, and the possible molecular mechanism was explored. PLAG1 3'-UTR sequences with five forms of base sequences plus poly(A) tails were detected in Hu sheep and the g.8795C>T mutation was identified. Luciferase reporter assay indicated that the g.8795C>T mutation influenced PLAG1 post-transcriptional activity. miRBase prediction showed that the g.8795C>T mutation was located in the miR-139 seed sequence binding region, and miR-139 overexpression significantly decreased both PLAG1-CC and PLAG1-TT activities. Moreover, the luciferase activity of PLAG1-CC was significantly lower than that of the PLAG1-TT, but miR-139 inhibition substantially increased both PLAG1-CC and PLAG1-TT luciferase activities, suggesting that PLAG1 is the target gene of miR-139. Thus, the g.8795C>T mutation upregulates PLAG1 expression by weakening its binding with miR-139, promoting PLAG1 expression, and increasing Hu sheep birth and weaning weights.

Exploring and Identifying Candidate Genes and Genomic Regions Related to Economically Important Traits in Hanwoo Cattle

The purpose of the current review was to explore and summarize different studies concerning the detection and characterization of candidate genes and genomic regions associated with economically important traits in Hanwoo beef cattle. Hanwoo cattle, the indigenous premium beef cattle of Korea, were introduced for their marbled fat, tenderness, characteristic flavor, and juiciness. To date, there has been a strong emphasis on the genetic improvement of meat quality and yields, such as backfat thickness (BFT), marbling score (MS), carcass weight (CW), eye muscle area (EMA), and yearling weight (YW), as major selection criteria in Hanwoo breeding programs. Hence, an understanding of the genetics controlling these traits along with precise knowledge of the biological mechanisms underlying the traits would increase the ability of the industry to improve cattle to better meet consumer demands. With the development of high-throughput genotyping, genomewide association studies (GWAS) have allowed the detection of chromosomal regions and candidate genes linked to phenotypes of interest. This is an effective and useful tool for accelerating the efficiency of animal breeding and selection. The GWAS results obtained from the literature review showed that most positional genes associated with carcass and growth traits in Hanwoo are located on chromosomes 6 and 14, among which LCORL, NCAPG, PPARGC1A, ABCG2, FAM110B, FABP4, DGAT1, PLAG1, and TOX are well known. In conclusion, this review study attempted to provide comprehensive information on the identified candidate genes associated with the studied traits and genes enriched in the functional terms and pathways that could serve as a valuable resource for future research in Hanwoo breeding programs.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Association of pleomorphic adenoma gene 1 with body weight and measurement of Bali cattle (Bos javanicus)

Background and Aim:

Pleomorphic adenoma gene 1 (PLAG1) encodes a multifunctional transcription factor that controls many genes and pathways and is associated with cattle body weight and measurements. This study aimed to evaluate the association between PLAG1 polymorphisms with body weight and measurements in Bali cattle.

Materials and Methods:

A total of 87 Bali cattle, consisting of 48 bulls and 39 heifers at the Breeding Center for Bali Cattle, were used as the population in this study. Cattle were 2 years old and kept semi-intensively in the pasture. Phenotype data consisting of body weight, withers height, body length, chest girth, waist height, and chest depth were measured. Birth weight data were obtained from birth records, and weight gain, adjusted weaning weight, and yearling weight were calculated using formulas. Blood samples were taken from the jugular vein as much as 5 mL, and genomic DNA was isolated using the salting-out method. Polymerase chain reaction (PCR) was performed to amplify three target polymorphisms, namely, g.48308 C>T, g.32212 (19 bp indel), and g.45233 T>C. The presence of a 19 bp indel was determined by direct observation of the PCR product on a 2% agarose gel. Two other polymorphisms were detected by PCR-restriction fragment length polymorphism using the restriction endonuclease enzymes SacII and BclI. PLAG1 genotype and phenotype associations were analyzed using a general linear model.

Results:

The results showed that two of the target polymorphisms in PLAG1 did not vary. The DD genotype indicated by 123 bp of PCR product was the only genotype identified for g.32212 19 bp indel, and TT genotype was the only genotype found for g.45233 T>C single-nucleotide polymorphism (SNP). Conversely, g.48308 C>T SNP was found to be polymorphic. In addition, the g.48308 C>T polymorphism of PLAG1 was significantly associated with body length of Bali cattle. Cattle with the CC genotype had a greater body length than the other two genotypes.

Conclusion:

The g.48308 C>T SNP in PLAG1 was associated with Bali cattle body length characteristics. This finding could be used as a basis for selecting Bali cattle based on body length characteristics.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

Genetic variation in PLAG1 is associated with early fertility in Australian Brahman cattle

Variation in the genome region coding for PLAG1 has well-documented associations with skeletal growth and age at puberty in cattle. However, the influence of PLAG1 on other economically important traits such as cow stayability has not yet been explored. Here we investigate the effect of PLAG1 variation on early and later in life female fertility, as well as size and growth, in a well phenotyped Australian Brahman herd. Yearly pregnancy and productivity records were collected from 2,839 genotyped Brahman cows and used to generate fertility, growth, and weight phenotypes. A variant on chromosome 14 in PLAG1 (NC_037341.1:g.23338890G>T, rs109815800) was previously determined to be a putative causative mutation associated with variation in cattle stature. The imputed PLAG1 genotype at this variant was isolated for each animal and the effect of PLAG1 genotype on each trait was estimated using linear modelling. Regardless of how heifer fertility was measured, there was a significant (P < 0.05) and desirable relationship between the additive effects of PLAG1 genotype and successful heifer fertility. Heifers with two copies of the alternate allele (TT) conceived earlier and had higher pregnancy and calving rates. However, the effects of PLAG1 genotype on fertility began to diminish as cows aged and did not significantly influence stayability at later ages. While there was no effect of genotype on growth, PLAG1 had a negative effect on mature cow weight (P < 0.01), where females with two copies of the alternate allele (TT) were significantly smaller than those with either one or none. Selection emphasis on improved Brahman heifer fertility will likely increase the frequency of the T allele of rs109815800, which may also increase herd profitability and long-term sustainability through improved reproductive efficiency and reduced mature cow size.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

A genome-wide association analysis for body weight at 35 days measured on 137,343 broiler chickens

BACKGROUND

Body weight (BW) is an economically important trait in the broiler (meat-type chickens) industry. Under the assumption of polygenicity, a "large" number of genes with "small" effects is expected to control BW. To detect such effects, a large sample size is required in genome-wide association studies (GWAS). Our objective was to conduct a GWAS for BW measured at 35 days of age with a large sample size.

METHODS

The GWAS included 137,343 broilers spanning 15 pedigree generations and 392,295 imputed single nucleotide polymorphisms (SNPs). A false discovery rate of 1% was adopted to account for multiple testing when declaring significant SNPs. A Bayesian ridge regression model was implemented, using AlphaBayes, to estimate the contribution to the total genetic variance of each region harbouring significant SNPs (1 Mb up/downstream) and the combined regions harbouring non-significant SNPs.

RESULTS

GWAS revealed 25 genomic regions harbouring 96 significant SNPs on 13 Gallus gallus autosomes (GGA1 to 4, 8, 10 to 15, 19 and 27), with the strongest associations on GGA4 at 65.67-66.31 Mb (Galgal4 assembly). The association of these regions points to several strong candidate genes including: (i) growth factors (GGA1, 4, 8, 13 and 14); (ii) leptin receptor overlapping transcript (LEPROT)/leptin receptor (LEPR) locus (GGA8), and the STAT3/STAT5B locus (GGA27), in connection with the JAK/STAT signalling pathway; (iii) T-box gene (TBX3/TBX5) on GGA15 and CHST11 (GGA1), which are both related to heart/skeleton development); and (iv) PLAG1 (GGA2). Combined together, these 25 genomic regions explained ~ 30% of the total genetic variance. The region harbouring significant SNPs that explained the largest portion of the total genetic variance (4.37%) was on GGA4 (~ 65.67-66.31 Mb).

CONCLUSIONS

To the best of our knowledge, this is the largest GWAS that has been conducted for BW in chicken to date. In spite of the identified regions, which showed a strong association with BW, the high proportion of genetic variance attributed to regions harbouring non-significant SNPs supports the hypothesis that the genetic architecture of BW35 is polygenic and complex. Our results also suggest that a large sample size will be required for future GWAS of BW35.

Detection of 15-bp Deletion Mutation within PLAG1 Gene and Its Effects on Growth Traits in Goats

Simple Summary

Goats have always served as an important domesticated livestock. PLAG1 is a major gene that affects the stature and growth of animals. Body size traits are very important for goats as it directly affects the economic characteristics of meat and cashmere production. This study showed that the 15-base pair (bp) InDel (rs637141549) can significantly affect growth traits such as body weight, height, height at hip cross, chest circumference, hip width and body index of goats through the detection of large samples (n = 1581) in four indigenous breeds. Accordingly, it is suggested that the deletion mutation can be used as a potential molecular marker that significantly affects goat growth traits. Moreover, the 15bp deletion mutation can be used as a potential molecular marker, which significantly affects the growth traits of goats and plays an important role in animal husbandry production.

Abstract

Stature and weight are important growth and development traits for animals, which also significantly affect the productivity of livestock. Polymorphic adenoma gene 1 (PLAG1) is located in the growth-related quantitative trait nucleotides (QTN), and its variation has been determined to significantly affect the body stature of bovines. This study found that novel 15-bp InDel could significantly influence important growth traits in goats. The frequencies of genotypes of the 15-bp mutation and relationship with core growth traits such as body weight, body height, height at hip cross, chest circumference, hip width and body index were explored in 1581 individuals among 4 Chinese native goat breeds. The most frequent genotypes of Shaanbei white Cashmere goat (SWCG), Inner Mongolia White Cashmere goat (IMCG) and Guanzhong Dairy goat (GZDG) were II genotypes (insertion/insertion), and the frequency of ID genotype (insertion/deletion) was found to be slightly higher than that of II genotype in Hainan Black goat (HNBG), showing that the frequency of the I allele was higher than that of the D allele. In adult goats, there were significant differences between 15-bp variation and body weight, chest circumference and body height traits in SWCG (p < 0.05). Furthermore, the locus was also found to be significantly correlated with the body index of HNBG (p = 0.044) and hip width in GZDG (p = 0.002). In regard to lambs, there were significant differences in height at the hip cross of SWCG (p = 0.036) and hip width in IMWC (p = 0.005). The corresponding results suggest that the 15-bp InDel mutation of PLAG1 is associated with the regulation of important growth characteristics of both adult and lamb of goats, which may serve as efficient molecular markers for goat breeding.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle

Simple Summary

Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest. Gene set enrichment analysis indicated that the identified candidate genes were related to biological processes and functional terms that were involved in growth and lipid metabolism. In conclusion, these results suggest that the incorporation of GWAS results and network analysis can help us to better understand the genetic bases underlying growth and carcass traits.

Abstract

In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.

Distribution and association study of PLAG1 gene between copy number variation and Chinese cattle populations

Copy number variation mainly refers to the copy number change of DNA fragments from 1 to 5 Mb. The deletion, duplication, inversion and ectopic of these fragments are collectively referred to as CNV. Numerous studies have shown that transfer factors play a vital role in regulating the growth and development of the body, for example the pleomorphic adenoma gene (PLAG). However, there is no study of CNV in PLAG1 gene. We qualified copy numbers within PLAG1 gene in 8 cattle breeds (Qinchuan, Qaidamu, Jinjiang, Guangfeng, Ji'an, Jiaxian, Pinan and Xianan cattle) by quantitative PCR, and explored their impacts on CNV of PLAG1 gene and phenotypic traits in Xianan cattle. We defined Deletion into CN = 0, Normal into CN = 1 and Duplication into CN = 2. The results showed that the individual with type of CN = 1 has a significant better effect on heart girth in JA cattle population (p < 0.01); the individual with type of CN = 1 and CN = 0 has a better effect on Rump length in JX cattle population (p < 0.05); the individual with type of CN = 0 has a better effect on cannon bone circumference in XN cattle population (p < 0.05). Association analysis showed that in JA cattle, the number of CN = 2 is great in JA cattle population, and the performance of CN = 2 in heart girth is better than CN = 1; in JX cattle, the rump length of CN = 2 is less than individual with CN = 0 and CN = 1; in XN cattle, individuals with CN = 0 have a better performance on cannon bone circumference than others. The results can provide a theoretical basis for molecular breeding of Chinese cattle, molecular mark-assist selection (MAS) of growth traits of Chinese cattle, and rapidly establish a Chinese cattle population with excellent genetic resources. Simple summaryWith the living standards rising, people's demand for beef is getting higher and higher, and there is a great significance to improve the growth performance of cattle. We measured body size data and detected copy number type of different cattle breeds (Xianan cattle, Ji'an cattle and Jiaxian cattle), and analyzed the correlation between the two object. We found that copy number variation of PLAG1 gene significantly affected some growth traits of XN cattle, JA cattle, and JX cattle. This may provide the basic material for molecular marker-assisted selection breeding of Chinese cattle breeds.

Breed- and trait-specific associations define the genetic architecture of calving performance traits in cattle

Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10⁻⁸) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10⁻⁶) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10⁻⁸). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

Genomic Breed Composition of Selection Signatures in Brangus Beef Cattle

Cattle breeding routinely uses crossbreeding between subspecies (Bos taurus taurus and Bos taurus indicus) to form composite breeds, such as Brangus. These composite breeds provide an opportunity to identify recent selection signatures formed in the new population and evaluate the genomic composition of these regions of the genome. Using high-density genotyping, we first identified runs of homozygosity (ROH) and calculated genomic inbreeding. Then, we evaluated the genomic composition of the regions identified as selected (selective sweeps) using a chromosome painting approach. The genomic inbreeding increased at approximately 1% per generation after composite breed formation, showing the need of inbreeding control even in composite breeds. Three selected regions in Brangus were also identified as Angus selection signatures. Two regions (chromosomes 14 and 21) were identified as signatures of selection in Brangus and both founder breeds. Five of the 10 homozygous regions in Brangus were predominantly Angus in origin (probability >80%), and the other five regions had a mixed origin but always with Brahman contributing less than 50%. Therefore, genetic events, such as drift, selection, and complementarity, are likely shaping the genetic composition of founder breeds in specific genomic regions. Such findings highlight a variety of opportunities to better control the selection process and explore heterosis and complementarity at the genomic level in composite breeds.

Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers

Age at first calving (AFC) plays an important role in the economic efficiency of beef cattle production. This trait can be affected by a combination of genetic and environmental factors, leading to physiological changes in response to heifers’ adaptation to a wide range of environments. Genome-wide association studies through the reaction norm model were carried out to identify genomic regions associated with AFC in Nellore heifers, raised under different environmental conditions (EC). The SNP effects for AFC were estimated in three EC levels (Low, Medium, and High, corresponding to average contemporary group effects on yearling body weight equal to 159.40, 228.6 and 297.6 kg, respectively), which unraveled shared and unique genomic regions for AFC in Low, Medium, and High EC levels, that varied according to the genetic correlation between AFC in different EC levels. The significant genomic regions harbored key genes that might play an important biological role in controlling hormone signaling and metabolism. Shared genomic regions among EC levels were identified on BTA 2 and 14, harboring candidate genes associated with energy metabolism (IGFBP2, IGFBP5, SHOX, SMARCAL1, LYN, RPS20, MOS, PLAG1, CHCD7, and SDR16C6). Gene set enrichment analyses identified important biological functions related to growth, hormone levels affecting female fertility, physiological processes involved in female pregnancy, gamete generation, ovulation cycle, and age at puberty. The genomic regions highlighted differences in the physiological processes linked to AFC in different EC levels and metabolic processes that support complex interactions between the gonadotropic axes and sexual precocity in Nellore heifers.

Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

The Montana Tropical^® Composite is a recently developed beef cattle population that is rapidly expanding in Brazil and other tropical countries. This is mainly due to its improved meat quality and adaptation to tropical climate conditions compared to Zebu and Taurine cattle breeds, respectively. This study aimed to investigate the genetic architecture of ultrasound-based carcass and meat quality traits in Montana Tropical^® Composite beef cattle. Therefore, we estimated variance components and genetic parameters and performed genome-wide association studies using the weighted single-step Genomic Best Linear Unbiased Prediction (GBLUP) approach. A pedigree dataset containing 28,480 animals was used, in which 1,436 were genotyped using a moderate-density Single Nucleotide Polymorphism panel (30K; 30,105 SNPs). A total of 9,358, 5,768, 7,996, and 1,972 phenotypic records for the traits Longissimus muscle area (LMA), backfat thickness (BFT), rump fat thickness (RFT), and for marbling score (MARB), respectively, were used for the analyses. Moderate to high heritability estimates were obtained and ranged from 0.16 ± 0.03 (RFT) to 0.33 ± 0.05 (MARB). A high genetic correlation was observed between BFT and RFT (0.97 ± 0.02), suggesting that a similar set of genes affects both traits. The most relevant genomic regions associated with LMA, BFT, RFT, and MARB were found on BTA10 (5.4–5.8 Mb), BTA27 (25.2–25.5 Mb), BTA18 (60.6–61.0 Mb), and BTA21 (14.8–15.4 Mb). Two overlapping genomic regions were identified for RFT and MARB (BTA13:47.9–48.1 Mb) and for BFT and RFT (BTA13:61.5–62.3 Mb). Candidate genes identified in this study, including PLAG1, LYN, WWOX, and PLAGL2, were previously reported to be associated with growth, stature, skeletal muscle growth, fat thickness, and fatty acid composition. Our results indicate that ultrasound-based carcass and meat quality traits in the Montana Tropical^® Composite beef cattle are heritable, and therefore, can be improved through selective breeding. In addition, various novel and already known genomic regions related to these traits were identified, which contribute to a better understanding of the underlying genetic background of LMA, BFT, RFT, and MARB in the Montana Tropical Composite population.

Generation of Retinaldehyde for Retinoic Acid Biosynthesis

The concentration of all-trans-retinoic acid, the bioactive derivative of vitamin A, is critically important for the optimal performance of numerous physiological processes. Either too little or too much of retinoic acid in developing or adult tissues is equally harmful. All-trans-retinoic acid is produced by the irreversible oxidation of all-trans-retinaldehyde. Thus, the concentration of retinaldehyde as the immediate precursor of retinoic acid has to be tightly controlled. However, the enzymes that produce all-trans-retinaldehyde for retinoic acid biosynthesis and the mechanisms responsible for the control of retinaldehyde levels have not yet been fully defined. The goal of this review is to summarize the current state of knowledge regarding the identities of physiologically relevant retinol dehydrogenases, their enzymatic properties, and tissue distribution, and to discuss potential mechanisms for the regulation of the flux from retinol to retinaldehyde.

A SNP in PLAG1 is associated with body height trait in Chinese cattle

Stature is an important quantitative trait for cattle performance, which influences herd productivity. Previous studies have reported that an SNP (AC_000171.1:g.25015640G>T, rs109815800) in Pleomorphic adenoma gene 1 (PLAG1) on chromosome 14 (CHR14) is associated with bovine stature. To validate whether rs109815800 is associated with the body height of Chinese cattle, we carried out an association analysis using 558 adult cattle samples from seven populations. Then, 1038 samples from 38 Chinese cattle breeds were used to show the geographical distribution of this variant in China. The results showed that the Q allele (G allele) increased the height of cattle. Furthermore, the frequencies of Q allele in Chinese native breeds tend to decrease from northern China to southern China, and the frequency of Q allele in two Chinese beef cattle breeds is much higher than that in another 36 Chinese local cattle breeds. Our data suggest that the prevalence of the Q allele is correlated with latitude in China.

Mice lacking the epidermal retinol dehydrogenases SDR16C5 and SDR16C6 display accelerated hair growth and enlarged meibomian glands

Retinol dehydrogenases catalyze the rate-limiting step in the biosynthesis of retinoic acid, a bioactive lipid molecule that regulates the expression of hundreds of genes by binding to nuclear transcription factors, the retinoic acid receptors. Several enzymes exhibit retinol dehydrogenase activities in vitro; however, their physiological relevance for retinoic acid biosynthesis in vivo remains unclear. Here, we present evidence that two murine epidermal retinol dehydrogenases, short-chain dehydrogenase/reductase family 16C member 5 (SDR16C5) and SDR16C6, contribute to retinoic acid biosynthesis in living cells and are also essential for the oxidation of retinol to retinaldehyde in vivo Mice with targeted knockout of the more catalytically active SDR16C6 enzyme have no obvious phenotype, possibly due to functional redundancy, because Sdr16c5 and Sdr16c6 exhibit an overlapping expression pattern during later developmental stages and in adulthood. Mice that lack both enzymes are viable and fertile but display accelerated hair growth after shaving and also enlarged meibomian glands, consistent with a nearly 80% reduction in the retinol dehydrogenase activities of skin membrane fractions from the Sdr16c5/Sdr16c6 double-knockout mice. The up-regulation of hair-follicle stem cell genes is consistent with reduced retinoic acid signaling in the skin of the double-knockout mice. These results indicate that the retinol dehydrogenase activities of murine SDR16C5 and SDR16C6 enzymes are not critical for survival but are responsible for most of the retinol dehydrogenase activity in skin, essential for the regulation of the hair-follicle cycle, and required for the maintenance of both sebaceous and meibomian glands.

Detection of InDel variations within seven candidate genes and their associations with phenotypic traits in three cattle breeds

The Xinjiang brown cattle, Red steppe cattle, and Yunling cattle are indigenous cultivated cattle breeds in Chinese frontier provinces, and they produce high-grade beef and milk products, however, their genetic diversity in many important genes related to excellent meat and milk production is still unknown. Our previous studies have found that several candidate genes (e.g., SREBP1c and PAX7) were associated with bovine economically important phenotypic traits, but none has been reported in the above-mentioned three cattle breeds. Since the InDel (insertion/deletion) marker becomes a useful tool applied in the animal molecular breeding, herein, we firstly found that the InDel variations of seven candidate genes in these cattle. Results showed that the genotypic and allelic distributions of these seven genes were remarkably different among these three cattle (p < 0.05 or p < 0.01). Furthermore, the InDel variations of SREBP1c and PAX7 genes were significantly associated with eight phenotypic traits in Xinjiang brown cattle (p < 0.05 or p < 0.01), respectively, suggesting that they can become the useful DNA markers.

Genome-wide association study identifies the PLAG1-OXR1 region on BTA14 for carcass meat yield in cattle

Carcass meat yield is an important carcass trait that contributes to the production efficiency and economic benefits in beef cattle. It is therefore critical to identify quantitative trait loci associated with carcass traits to enable selection. Our previous studies have identified several causal variants within the pleomorphic adenoma gene 1 ( PLAG1) and coiled-coil-helix-coiled-coil-helix domain-containing 7 ( CHCHD7) genes on BTA14 for carcass traits in Chinese Simmental. In the current study, we carried out a genome-wide association study for carcass meat yield in 472 Wagyu cattle with Bovine HD SNP array. Our results showed that 27 single nucleotide polymorphisms (SNPs) were identified for tenderloin weight (TDW), striploin weight (SPW), chuck roll weight (CRW), bicep weight (BPW), knuckle weight (KCW), and flank steak weight (FSW) in Wagyu cattle. Of these SNPs, 10 distinct SNPs were detected within the oxidation resistance 1 ( OXR1), fatty acid binding protein 5 ( FABP5), TNF receptor superfamily member 11b ( TNFRSF11B), and zinc finger CCCH-type containing 3 ( ZC3H3) genes on BTA14. Notably, three significant SNPs, BovineHD1400016738, BovineHD1400016743, and BovineHD1400016665 within OXR1, were shown strong linkage disequilibrium (r2 > 0.8) and significantly associated with CRW ( P = 1.37 × 10−8 ~ 1.94 × 10−8). Moreover, Ingenuity Pathway Analysis showed that OXR1, FABP5, and CAP1A genes were involved in a single network and FABP5 may regulate the expression of OXR1 gene via node gene, peroxisome proliferator-activated receptor gamma ( PPARG). Overall, this study suggests that OXR1 and FABP5 are candidate genes affecting carcass traits in Wagyu and the PLAG1-OXR1 region on BTA14 as a putative susceptibility locus for carcass meat yield for both Chinese Simmental and Wagyu.

A novel SNP of PLAG1 gene and its association with growth traits in Chinese cattle

Genetic polymorphism has great influences on the improvement of cattle traits. The polymorphism of the same gene family will greatly change the growth traits of cattle, such as the pleomorphic adenoma gene (PLAG) family. Many studies have shown that the PLAG family proteins are the transcription regulators of nuclear protein, which mainly regulates the expression of many important genes in the body. In cattle, single nucleotide polymorphisms (SNPs) within or near the PLAG1 gene is associated with economic traits of height in cattle. Here we investigate a SNPs of bovine PLAG1 in 646 cattle from five breeds. We found three different genotypes by using cleaved amplification polymorphism sequence-tagged sites (CAPs), includes some significant differences in body height, chest circumference and other shapes (P < 0.05), also we found that the TT genotype had no advantage in body shape. These results indicate that the selection of PLAG1 gene could be used to ensure the breeding direction for growth traits of the beef cattle.

Exploring genotype-phenotype relationships of the CRABP2 gene on growth traits in beef cattle

Cellular retinoic acid binding protein 2 (CRABP2) is essential to myoblast differentiation. However, there was no report about the function of CRABP2 gene in cattle. This study explored the association of CRABP2 gene polymorphisms with growth traits in cattle breeds by several methods, such as DNA sequencing, PCR, PCR-RFLP and forced PCR-RFLP. Two sequence variants were determined. There were 621 individuals in six cattle breeds from China for the experiment, and three breeds were used to test validation of polymorphisms and extent of linkage disequilibrium (LD). The results showed that both SNPs (SNP1, g.2458 G > T, SNP2, g.3878 G > A) were in intron1. Two SNPs were in low linkage disequilibrium. Association analysis suggested that SNP1 had the significant difference on growth traits with body height, height at hip cross and body slanting length (P < .05), while SNP2 showed a significant difference in growth traits with body height, height at hip cross and body slanting length(P < .05). The results of this investigation displayed that the CRABP2 gene is an available candidate gene and may be used for breed selection and conservation.

Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle

Multitrait meta-analyses are a strategy to produce more accurate genome-wide association studies, especially for complex phenotypes. We carried out a meta-analysis study for traits related to sexual precocity in tropical beef cattle (Nellore and Brahman) aiming to identify important genomic regions affecting these traits. The traits included in the analyses were age at first calving (AFC), early pregnancy (EP), age at first corpus luteum (AGECL), first postpartum anoestrus interval (PPAI), and scrotal circumference (SC). The traits AFC, EP, and SCN were measured in Nellore cattle, while AGECL, PPAI, and SCB were measured in Brahman cattle. Meta-analysis resulted in 108 significant single-nucleotide polymorphisms (SNPs), at an empirical threshold P-value of 1.39 × 10-5 (false discovery rate [FDR] < 0.05). Within 0.5 Mb of the significant SNP, candidate genes were annotated and analyzed for functional enrichment. Most of the closest genes to the SNP with higher significance in each chromosome have been associated with important roles in reproductive function. They are TSC22D2, KLF7, ARHGAP29, 7SK, MAP3K5, TLE3, WDR5, TAF3, TMEM68, PPP1R15B, NR2F2, GALR1, SUFU, and KCNU1. We did not observe any significant SNP in BTA5, BTA12, BTA17, BTA18, BTA19, BTA20, BTA22, BTA23, BTA25, and BTA28. Although the majority of significant SNPs are in BTA14, it was identified significant associations in multiple chromosomes (19 out of 29 autosomes), which is consistent with the postulation that reproductive traits are complex polygenic phenotypes. Five proposed association regions harbor the majority of the significant SNP (76%) and were distributed over four chromosomes (P < 1.39 × 10-5, FDR < 0.05): BTA2 (5.55%) from 95 to 96 Mb, BTA4 (5.55%) from 94.1 to 94.8 Mb, BTA14 (59.26%) from 24 to 25 Mb and 29 to 30 Mb, and BTA21 (5.55%) from 6.7 Mb to 11.4 Mb. These regions harbored key genes related to reproductive function. Moreover, these genes were enriched for functional groups associated with immune response, maternal-fetal tolerance, pregnancy maintenance, embryo development, fertility, and response to stress. Further studies including other breeds and precocity traits could confirm the importance of these regions and identify new candidate regions for sexual precocity in beef cattle.

Functional Partitioning of Genomic Variance and Genome-Wide Association Study for Carcass Traits in Korean Hanwoo Cattle Using Imputed Sequence Level SNP Data

Quantitative traits are usually controlled by numerous genomic variants with small individual effects, and variances associated with those traits are explained in a continuous manner. However, the relative contributions of genomic regions to observed genetic variations have not been well explored using sequence level single nucleotide polymorphism (SNP) information. Here, imputed sequence level SNP data (11,278,153 SNPs) of 2109 Hanwoo steers (Korean native cattle) were partitioned according to functional annotation, chromosome, and minor allele frequency (MAF). Genomic relationship matrices (GRMs) were constructed for each classified region and fitted in the model both separately and together for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BFT), and marbling score (MS) traits. A genome-wide association study (GWAS) was performed to identify significantly associated variants in genic and exon regions using a linear mixed model, and the genetic contribution of each exonic SNP was determined using a Bayesian mixture model. Considering all SNPs together, the heritability estimates for CWT, EMA, BFT, and MS were 0.57 ± 0.05, 0.46 ± 0.05, 0.45 ± 0.05, and 0.49 ± 0.05, respectively, which reflected substantial genomic contributions. Joint analysis revealed that the variance explained by each chromosome was proportional to its physical length with weak linear relationships for all traits. Moreover, genomic variances explained by functional category and MAF class differed greatly among the traits studied in joint analysis. For example, exon regions had larger contributions for BFT (0.13 ± 0.08) and MS (0.22 ± 0.08), whereas intron and intergenic regions explained most of the total genomic variances for CWT and EMA (0.22 ± 0.09–0.32 ± 0.11). Considering different functional classes of exon regions and the per SNP contribution revealed the largest proportion of genetic variance was attributable to synonymous variants. GWAS detected 206 and 27 SNPs in genic and exon regions, respectively, on BTA4, BTA6, and BTA14 that were significantly associated with CWT and EMA. These SNPs were harbored by 31 candidate genes, among which TOX, FAM184B, PPARGC1A, PRKDC, LCORL, and COL1A2 were noteworthy. BayesR analysis found that most SNPs (>93%) had very small effects and the 4.02–6.92% that had larger effects (10^-4 × σA2, 10^-3 × σA2, and 10^-2 × σA2) explained most of the total genetic variance, confirming polygenic components of the traits studied.

A PLAG1 mutation contributed to stature recovery in modern cattle

The recent evolution of cattle is marked by fluctuations in body size. Height in the Bos taurus lineage was reduced by a factor of ~1.5 from the Neolithic to the Middle Ages, and increased again only during the Early Modern Ages. Using haplotype analysis, we found evidence that the bovine PLAG1 mutation (Q) with major effects on body size, weight and reproduction is a >1,000 years old derived allele that increased rapidly in frequency in Northwestern European B. taurus between the 16^th and 18^th centuries. Towards the 19^th and 20^th centuries, Q was introgressed into non-European B. taurus and Bos indicus breeds. These data implicate a major role of Q in recent changes in body size in modern cattle, and represent one of the first examples of a genomic sweep in livestock that was driven by selection on a complex trait.

Novel Muscle-Specific Genes TCAP, TNNI1, and FHL1 in Cattle: SNVs, Linkage Disequilibrium, Combined Genotypes, Association Analysis of Growth Performance, and Carcass Quality Traits and Expression Studies

TCAP, TNNI1, and FHL1 regulate muscle growth and development. In this study, four single nucleotide variants (SNVs) were discovered in almost all of the exon and intron regions of the TCAP, TNNI1, and FHL1 genes using DNA pooled sequencing, polymerase chain reaction (PCR)-RFLP, and forced-PCR-RFLP methods in 576 cattle. Four SNVs were significantly associated with the growth performance and carcass quality traits of the cattle. In addition, the haplotype, haplotype frequency, and linkage disequilibrium coefficient of three sequence variants were also evaluated in the cattle population. Haplotype analysis demonstrated that eight haplotypes were present in the Qinchuan cattle population and no haplotypes were present in the Chinese Holstein population; haplotype 1 had the highest frequency in the Qinchuan (42.7%) population. Statistical analyses of 12 combined genotypes indicated that some were significantly associated with the growth performance and carcass quality traits of the Qinchuan cattle population. Moreover, the quantitative real-time polymerase chain reaction results demonstrated that the bovine TCAP, TNNI1, and FHL1 genes were exclusively expressed in muscle tissue. These data support the high potentials of the TCAP, TNNI1, and FHL1 as marker genes to improve the growth performance and carcass quality traits of Qinchuan cattle or other animals selection programs.

Deciphering signature of selection affecting beef quality traits in Angus cattle

Artificial selection towards a desired phenotype/trait has modified the genomes of livestock dramatically that generated breeds that greatly differ in morphology, production and environmental adaptation traits. Angus cattle are among the famous cattle breeds developed for superior beef quality. This paper aimed at exploring genomic regions under selection in Angus cattle that are associated with meat quality traits and other associated phenotypes. The whole genome of 10 Angus cattle was compared with 11 Hanwoo (A-H) and 9 Jersey (A-J) cattle breeds using a cross-population composite likelihood ratio (XP-CLR) statistical method. The top 1% of the empirical distribution was taken as significant and annotated using UMD3.1. As a result, 255 and 210 genes were revealed under selection from A-H and A-J comparisons, respectively. The WebGestalt gene ontology analysis resulted in sixteen (A-H) and five (A-J) significantly enriched KEGG pathways. Several pathways associated with meat quality traits (insulin signaling, type II diabetes mellitus pathway, focal adhesion pathway, and ECM-receptor interaction), and feeding efficiency (olfactory transduction, tight junction, and metabolic pathways) were enriched. Genes affecting beef quality traits (e.g., FABP3, FTO, DGAT2, ACS, ACAA2, CPE, TNNI1), stature and body size (e.g., PLAG1, LYN, CHCHD7, RPS20), fertility and dystocia (e.g., ESR1, RPS20, PPP2R1A, GHRL, PLAG1), feeding efficiency (e.g., PIK3CD, DNAJC28, DNAJC3, GHRL, PLAG1), coat color (e.g., MC1-R) and genetic disorders (e.g., ITGB6, PLAG1) were found to be under positive selection in Angus cattle. The study identified genes and pathways that are related to meat quality traits and other phenotypes of Angus cattle. The findings in this study, after validation using additional or independent dataset, will provide useful information for the study of Angus cattle in particular and beef cattle in general.

Exploring evidence of positive selection signatures in cattle breeds selected for different traits

Since domestication, the genome landscape of cattle has been changing due to natural and artificial selection forces resulting in several general and specialized cattle breeds of the world. Identifying genomic regions affected due to these forces in livestock gives an insight into the history of selection for economically important traits and genetic adaptation to specific environments of the populations under consideration. This study explores the genes/genomic regions under selection in relation to the phenotypes of Holstein, Hanwoo, and N'Dama cattle breeds using Tajima's D, XP-CLR, and XP-EHH population statistical methods. The whole genomes of 10 Holstein (South Korea), 11 Hanwoo (South Korea), and 10 N'Dama (West Africa-Guinea) cattle breeds re-sequenced to ~11x coverage and retained 37 million SNPs were used for the study. Selection signature analysis revealed 441, 512, and 461 genes under selection from Holstein, Hanwoo, and N'Dama cattle breeds, respectively. Among all these, seven genes including ARFGAP3, SNORA70, and other RNA genes were common between the breeds. From each of the gene lists, significant functional annotation cluster terms including milk protein and thyroid hormone signaling pathway (Holstein), histone acetyltransferase activity (Hanwoo), and renin secretion (N'Dama) were enriched. Genes that are related to the phenotypes of the respective breeds were also identified. Moreover, significant breed-specific missense variants were identified in CSN3, PAPPA2 (Holstein), C1orf116 (Hanwoo), and COMMD1 (N'Dama) genes. The genes identified from this study provide an insight into the biological mechanisms and pathways that are important in cattle breeds selected for different traits of economic significance.

Functional confirmation of PLAG1 as the candidate causative gene underlying major pleiotropic effects on body weight and milk characteristics

A major pleiotropic quantitative trait locus (QTL) located at ~25 Mbp on bovine chromosome 14 affects a myriad of growth and developmental traits in Bos taurus and indicus breeds. These QTL have been attributed to two functional variants in the bidirectional promoter of PLAG1 and CHCHD7. Although PLAG1 is a good candidate for mediating these effects, its role remains uncertain given that these variants are also associated with expression of five additional genes at the broader locus. In the current study, we conducted expression QTL (eQTL) mapping of this region using a large, high depth mammary RNAseq dataset representing 375 lactating cows. Here we show that of the seven previously implicated genes, only PLAG1 and LYN are differentially expressed by QTL genotype, and only PLAG1 bears the same association signature of the growth and body weight QTLs. For the first time, we also report significant association of PLAG1 genotype with milk production traits, including milk fat, volume, and protein yield. Collectively, these data strongly suggest PLAG1 as the causative gene underlying this diverse range of traits, and demonstrate new effects for the locus on lactation phenotypes.

Genome-Wide Association Study Reveals the PLAG1 Gene for Knuckle, Biceps and Shank Weight in Simmental Beef Cattle

Carcass traits of beef cattle have been genetically improved to increase yield of high quality meat. Genome-wide association study (GWAS) is a powerful method to identify genetic variants associated with carcass traits. For the 770K genotyped SNPs from 1141 Chinese Simmental cattle, we used the compressed mixed linear model (CMLM) to perform a genome-wide association study for knuckle, biceps and shank of beef carcass traits. Seventeen significantly associated SNPs were found, which are located on BTA6, BTA14 and BTA15. Interestingly, one pleiotropic quantitative trait nucleotide (QTN), named BovineHD1400007259 (p < 10⁻⁸) within the well-known gene region PLAG1-CHCHD7 on BTA14, was found to govern variation of the knuckle, biceps and shank traits. The QTN accounted for 8.6% of phenotypic variance for biceps. In addition, 16 more SNPs distributed on BTA14 were detected as being associated with the carcass traits.

Characterization of human short chain dehydrogenase/reductase SDR16C family members related to retinol dehydrogenase 10

All-trans-retinoic acid (RA) is a bioactive derivative of vitamin A that serves as an activating ligand for nuclear transcription factors, retinoic acid receptors. RA biosynthesis is initiated by the enzymes that oxidize retinol to retinaldehyde. It is well established that retinol dehydrogenase 10 (RDH10, SDR16C4), which belongs to the 16C family of the short chain dehydrogenase/reductase (SDR) superfamily of proteins, is the major enzyme responsible for the oxidation of retinol to retinaldehyde for RA biosynthesis during embryogenesis. However, several lines of evidence point towards the existence of additional retinol dehydrogenases that contribute to RA biosynthesis in vivo. In close proximity to RDH10 gene on human chromosome 8 are located two genes that are phylogenetically related to RDH10. The predicted protein products of these genes, retinol dehydrogenase epidermal 2 (RDHE2, SDR16C5) and retinol dehydrogenase epidermal 2-similar (RDHE2S, SDR16C6), share 59% and 56% sequence similarity with RDH10, respectively. Previously, we showed that the single ortholog of the human RDHE2 and RDHE2S in frogs, Xenopus laevis rdhe2, oxidizes retinol to retinaldehyde and is essential for frog embryonic development. In this study, we explored the potential of each of the two human proteins to contribute to RA biosynthesis. The results of this study demonstrate that human RDHE2 exhibits a relatively low but reproducible activity when expressed in either HepG2 or HEK293 cells. Expression of the native RDHE2 is downregulated in the presence of elevated levels of RA. On the other hand, the protein encoded by the human RDHE2S gene is unstable when expressed in HEK293 cells. RDHE2S protein produced in Sf9 cells is stable but has no detectable catalytic activity towards retinol. We conclude that the human RDHE2S does not contribute to RA biosynthesis, whereas the low-activity RA-sensitive human RDHE2 may have a role in adjusting the cellular levels of RA in accord with specific physiological conditions.

Pleiotropic Genes Affecting Carcass Traits in Bos indicus (Nellore) Cattle Are Modulators of Growth

Two complementary methods, namely Multi-Trait Meta-Analysis and Versatile Gene-Based Test for Genome-wide Association Studies (VEGAS), were used to identify putative pleiotropic genes affecting carcass traits in Bos indicus (Nellore) cattle. The genotypic data comprised over 777,000 single-nucleotide polymorphism markers scored in 995 bulls, and the phenotypic data included deregressed breeding values (dEBV) for weight measurements at birth, weaning and yearling, as well visual scores taken at weaning and yearling for carcass finishing precocity, conformation and muscling. Both analyses pointed to the pleomorphic adenoma gene 1 (PLAG1) as a major pleiotropic gene. VEGAS analysis revealed 224 additional candidates. From these, 57 participated, together with PLAG1, in a network involved in the modulation of the function and expression of IGF1 (insulin like growth factor 1), IGF2 (insulin like growth factor 2), GH1 (growth hormone 1), IGF1R (insulin like growth factor 1 receptor) and GHR (growth hormone receptor), suggesting that those pleiotropic genes operate as satellite regulators of the growth pathway.

Genome-Wide Association Study of Meat Quality Traits in Nellore Cattle

The objective of this study was to identify genomic regions that are associated with meat quality traits in the Nellore breed. Nellore steers were finished in feedlots and slaughtered at a commercial slaughterhouse. This analysis included 1,822 phenotypic records of tenderness and 1,873 marbling records. After quality control, 1,630 animals genotyped for tenderness, 1,633 animals genotyped for marbling, and 369,722 SNPs remained. The results are reported as the proportion of variance explained by windows of 150 adjacent SNPs. Only windows with largest effects were considered. The genomic regions were located on chromosomes 5, 15, 16 and 25 for marbling and on chromosomes 5, 7, 10, 14 and 21 for tenderness. These windows explained 3,89% and 3,80% of the additive genetic variance for marbling and tenderness, respectively. The genes associated with the traits are related to growth, muscle development and lipid metabolism. The study of these genes in Nellore cattle is the first step in the identification of causal mutations that will contribute to the genetic evaluation of the breed.

Polymorphisms and genes associated with puberty in heifers

Puberty onset is a multifactorial process influenced by genetic determinants and environmental conditions, especially nutritional status. Genes, genetic variations, and regulatory networks compose the molecular basis of achieving puberty. In this article, we reviewed the discovery of multiple polymorphisms and genes associated with heifer puberty phenotypes and discuss the opportunities to use this evolving knowledge of genetic determinants for breeding early pubertal Bos indicus-influenced cattle. The discovery of polymorphisms and genes was mainly achieved through candidate gene studies, quantitative trait loci analyses, genome-wide association studies, and recently, global gene expression studies (transcriptome). These studies are recapitulated and summarized in the current review.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

Emerging role of PLAG1 as a regulator of growth and reproduction

Pleomorphic adenoma gene 1 (PLAG1) belongs to the PLAG family of zinc finger transcription factors along with PLAG-like 1 and PLAG-like 2. The PLAG1 gene is best known as an oncogene associated with certain types of cancer, most notably pleomorphic adenomas of the salivary gland. While the mechanisms of PLAG1-induced tumorigenesis are reasonably well understood, the role of PLAG1 in normal physiology is less clear. It is known that PLAG1 is involved in cell proliferation by directly regulating a wide array of target genes, including a number of growth factors such as insulin-like growth factor 2. This is likely to be a central mode of action for PLAG1 both in embryonic development and in cancer. The phenotype of Plag1 knockout mice suggests an important role for PLAG1 also in postnatal growth and reproduction, as PLAG1 deficiency causes growth retardation and reduced fertility. A role for PLAG1 in growth and reproduction is further corroborated by genome-wide association studies in humans and domestic animals in which polymorphisms in the PLAG1 genomic region are associated with body growth and reproductive traits. Here we review the current evidence for PLAG1 as a regulator of growth and fertility and discuss possible endocrine mechanisms involved.