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

Genomic patterns of selection in morphometric traits across diverse Indian cattle breeds

This study seeks a comprehensive exploration of genome-wide selective processes impacting morphometric traits across diverse cattle breeds, utilizing an array of statistical methods. Morphometric traits, encompassing both qualitative and quantitative variables, play a pivotal role in characterizing and selecting livestock breeds based on their external appearance, size, and physical attributes. While qualitative traits, such as color, horn structure, and coat type, contribute to adaptive features and breed identification, quantitative traits like body weight and conformation measurements bear a closer correlation with production characteristics. This study employs advanced genotyping technologies, including the Illumina BovineSNP50 Bead Chip and next-generation sequencing methods like Reduced Representation sequencing, to identify genomic signatures associated with these traits. We applied four intra-population methods to find evidence of selection, such as Tajima's D, CLR, iHS, and ROH. We found a total of 40 genes under the selection signature, that were associated with morphometric traits in five cattle breeds (Kankrej, Tharparkar, Nelore, Sahiwal, and Gir). Crucial genes such as ADIPDQ, DPP6, INSIG1, SLC35D2 in Kankrej, LPL, ATP6V1B2, CDC14B in Tharparkar, HPSE2, PLAG1 in Nelore, PCSK1, PRKD1 in Sahiwal, and GNAQ, HPCAL1 in Gir were identified in our study. This approach provides valuable insights into the genetic basis of variations in body weight and conformation traits, facilitating informed selection processes and offering a deeper understanding of the evolutionary and domestication processes in diverse cattle breeds.

Beef Cattle Genome Project: Advances in Genome Sequencing, Assembly, and Functional Genes Discovery

Beef is a major global source of protein, playing an essential role in the human diet. The worldwide production and consumption of beef continue to rise, reflecting a significant trend. However, despite the critical importance of beef cattle resources in agriculture, the diversity of cattle breeds faces severe challenges, with many breeds at risk of extinction. The initiation of the Beef Cattle Genome Project is crucial. By constructing a high-precision functional annotation map of their genome, it becomes possible to analyze the genetic mechanisms underlying important traits in beef cattle, laying a solid foundation for breeding more efficient and productive cattle breeds. This review details advances in genome sequencing and assembly technologies, iterative upgrades of the beef cattle reference genome, and its application in pan-genome research. Additionally, it summarizes relevant studies on the discovery of functional genes associated with key traits in beef cattle, such as growth, meat quality, reproduction, polled traits, disease resistance, and environmental adaptability. Finally, the review explores the potential of telomere-to-telomere (T2T) genome assembly, structural variations (SVs), and multi-omics techniques in future beef cattle genetic breeding. These advancements collectively offer promising avenues for enhancing beef cattle breeding and improving genetic traits.

Weighted single-step genome-wide association study for direct and maternal genetic effects associated with birth and weaning weights in sheep

Body weight is an important economic trait for sheep meat production, and its genetic improvement is considered one of the main goals in the sheep breeding program. Identifying genomic regions that are associated with growth-related traits accelerates the process of animal breeding through marker-assisted selection, which leads to increased response to selection. In this study, we conducted a weighted single-step genome-wide association study (WssGWAS) to identify potential candidate genes for direct and maternal genetic effects associated with birth weight (BW) and weaning weight (WW) in Baluchi sheep. The data used in this research included 13,408 birth and 13,170 weaning records collected at Abbas-Abad Baluchi Sheep Breeding Station, Mashhad-Iran. Genotypic data of 94 lambs genotyped by Illumina 50K SNP BeadChip for 54,241 markers were used. The proportion of variance explained by genomic windows was calculated by summing the variance of SNPs within 1 megabase (Mb). The top 10 window genomic regions explaining the highest percentages of additive and maternal genetic variances were selected as candidate window genomic regions associated with body weights. Our findings showed that for BW, the top-ranked genomic regions (1 Mb windows) explained 4.30 and 4.92% of the direct additive and maternal genetic variances, respectively. The direct additive genetic variance explained by the genomic window regions varied from 0.31 on chromosome 1 to 0.59 on chromosome 8. The highest (0.84%) and lowest (0.32%) maternal genetic variances were explained by genomic windows on chromosome 10 and 17, respectively. For WW, the top 10 genomic regions explained 6.38 and 5.76% of the direct additive and maternal genetic variances, respectively. The highest and lowest contribution of direct additive genetic variances were 1.37% and 0.42%, respectively, both explained by genomic regions on chromosome 2. For maternal effects on WW, the highest (1.38%) and lowest (0.41%) genetic variances were explained by genomic windows on chromosome 2. Further investigation of these regions identified several possible candidate genes associated with body weight. Gene ontology analysis using the DAVID database identified several functional terms, such as translation repressor activity, nucleic acid binding, dehydroascorbic acid transporter activity, growth factor activity and SH2 domain binding.

Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data

OBJECTIVE

The objective of this study was to investigate the genetic diversity, population structure and whole-genome selection signatures of Luxi cattle to reveal its genomic characteristics in terms of meat and carcass traits, skeletal muscle development, body size, and other traits.

METHODS

To further analyze the genomic characteristics of Luxi cattle, this study sequenced the whole-genome of 16 individuals from the core conservation farm in Shandong region, and collected 174 published genomes of cattle for conjoint analysis. Furthermore, three different statistics (pi, Fst, and XP-EHH) were used to detect potential positive selection signatures related to selection in Luxi cattle. Moreover, gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were performed to reveal the potential biological function of candidate genes harbored in selected regions.

RESULTS

The results showed that Luxi cattle had high genomic diversity and low inbreeding levels. Using three complementary methods (pi, Fst, and XP-EHH) to detect the signatures of selection in the Luxi cattle genome, there were 2,941, 2,221 and 1,304 potentially selected genes identified, respectively. Furthermore, there were 45 genes annotated in common overlapping genomic regions covered 0.723 Mb, including PLAG1 zinc finger (PLAG1), dedicator of cytokinesis 3 (DOCK3), ephrin A2 (EFNA2), DAZ associated protein 1 (DAZAP1), Ral GTPase activating protein catalytic subunit alpha 1 (RALGAPA1), mediator complex subunit 13 (MED13), and decaprenyl diphosphate synthase subunit 2 (PDSS2), most of which were enriched in pathways related to muscle growth and differentiation and immunity.

CONCLUSION

In this study, we provided a series of genes associated with important economic traits were found in positive selection regions, and a scientific basis for the scientific conservation and genetic improvement of Luxi cattle.

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 for primal cut lean traits in Canadian beef cattle

This study identified genomic variants and underlying candidate genes related to the whole carcass and individual primal cut lean content in Canadian commercial crossbred beef cattle. Genotyping information of 1035 crossbred beef cattle were available alongside estimated and actual carcass lean meat yield and individual primal cut lean content in all carcasses. Significant fixed effects and covariates were identified and included in the animal model. Genome-wide association analysis were implemented using the weighted single-step genomic best linear unbiased prediction (WssGBLUP). A number of candidate genes identified linked to lean tissue production were unrelated to estimated lean meat yield and were specific to the actual lean traits. Among these, 41 genes were common for actual lean traits, on specific regions of BTA4, BTA13 and BTA25 indicating potential involvement in lean mass synthesis. Therefore, the results suggested the inclusion of primal cut lean traits as a selection objective in breeding programs with consideration of further functional studies of the identified genes could help in optimizing lean yield for maximal carcass value.

Goat Pleomorphic Adenoma Gene 1 (PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits

The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.

Assessing Genomic Diversity and Signatures of Selection in Chinese Red Steppe Cattle Using High-Density SNP Array

Chinese Red Steppe Cattle (CRS), a composite cattle breed, is well known for its milk production, high slaughter rate, carcass traits, and meat quality. Nowadays, it is widely bred in Jilin and Hebei Province and the Inner Mongolia Autonomous region. However, the population structure and the genetic basis of prominent characteristics of CRS are still unknown. In this study, we systematically describe their population structure, genetic diversity, and selection signature based on genotyping data from 61 CRS individuals with GGP Bovine 100 K chip. The results showed that CRS cattle had low inbreeding levels and had formed a unique genetic structure feature. Using two complementary methods (including comprehensive haplotype score and complex likelihood ratio), we identified 1291 and 1285 potentially selected genes, respectively. There were 141 genes annotated in common 106 overlapping genomic regions covered 5.62 Mb, including PLAG1, PRKG2, DGAT1, PARP10, TONSL, ADCK5, and BMP3, most of which were enriched in pathways related to muscle growth and differentiation, milk production, and lipid metabolism. This study will contribute to understanding the genetic mechanism behind artificial selection and give an extensive reference for subsequent breeding.

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.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

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.

Indel mutations of sheep PLAG1 gene and their associations with growth traits

Pleiomorphic adenoma gene 1 (PLAG1) is mainly expressed in embryonic development, and it is reported to take an effect on the growth performance of mice, cattle, pigs, and sheep. To explore how conservative the PLAG1 is in different sheep breeds, the effects of the two indel variants on the growth traits of the Chinese Luxi blackhead (LXBH) sheep were firstly detected. The P2-del 30 bp and P4-del 45 bp indel loci of the sheep PLAG1 gene were significantly related to 15 growth traits (P < 0.05). Genotype ID and genotype II were dominant for the P2-del 30 bp and P4-del 45 bp loci, respectively. The above findings indicated that the two indel mutations in the ovine PLAG1 gene were suggested to become the molecular markers for the selection of economic traits in sheep.

Pleomorphic adenoma gene 1 (PLAG1) promotes proliferation and inhibits apoptosis of bovine primary myoblasts through the PI3K-Akt signaling pathway

Pleomorphic adenoma gene 1 (PLAG1) is a transcription factor involved in various cellular processes in organismal growth and development. However, its role in muscle function is unclear. This work investigated the roles of PLAG1 in muscle development and explored its regulatory mechanisms. The PLAG1 was proved to promote the proliferation of bovine primary myoblasts using the cell counting kit 8 (CCK-8) assay (P < 0.001), 5-ethynyl-2'-deoxyuridine (EdU) proliferation assay (P = 0.005), quantitative real-time polymerase chain reaction (qRT-PCR) (P = 0.028), western blot, and flow cytometry (P < 0.05), and to inhibit apoptosis of bovine primary myoblasts using qRT-PCR (P = 0.038), western blot, and flow cytometry (P < 0.001). Chromatin immunoprecipitation sequencing (ChIP-seq) and western blot showed PLAG1 upregulated phosphorylated (p)-PI3K, PI3K, p-Akt, Akt, Cyclin D1, and CDK2 and inhibited the expression of p21 and p27 to enhance myoblast proliferation, and increased expression of Bcl-2, and Bcl-xL to inhibit apoptosis. Additionally, PLAG1 was identified as a target of miR-1 using dual-luciferase assay (P < 0.001), qRT-PCR (P < 0.001), and western blot. Furthermore, miR-1 might be a potential mediator of the positive feedback regulation relationship between PLAG1 and the PI3K-Akt signaling pathway.

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.

Identification of Candidate Variants Associated With Bone Weight Using Whole Genome Sequence in Beef Cattle

Bone weight is critical to affect body conformation and stature in cattle. In this study, we conducted a genome-wide association study for bone weight in Chinese Simmental beef cattle based on the imputed sequence variants. We identified 364 variants associated with bone weight, while 350 of them were not included in the Illumina BovineHD SNP array, and several candidate genes and GO terms were captured to be associated with bone weight. Remarkably, we identified four potential variants in a candidate region on BTA6 using Bayesian fine-mapping. Several important candidate genes were captured, including LAP3, MED28, NCAPG, LCORL, SLIT2, and IBSP, which have been previously reported to be associated with carcass traits, body measurements, and growth traits. Notably, we found that the transcription factors related to MED28 and LCORL showed high conservation across multiple species. Our findings provide some valuable information for understanding the genetic basis of body stature in beef cattle.

Assessing the Genetic Background and Selection Signatures of Huaxi Cattle Using High-Density SNP Array

Huaxi cattle, a specialized beef cattle breed in China, has the characteristics of fast growth, high slaughter rate, and net meat rate, good reproductive performance, strong stress resistance, and wide adaptability. In this study, we evaluated the genetic diversity, population structure, and genetic relationships of Huaxi cattle and its ancestor populations at the genome-wide level, as well as detecting the selection signatures of Huaxi cattle. Principal component analysis (PCA) and phylogenetic analysis revealed that Huaxi cattle were obviously separated from other cattle populations. The admixture analysis showed that Huaxi cattle has distinct genetic structures among all populations at K = 4. It can be concluded that Huaxi cattle has formed its own unique genetic features. Using integrated haplotype score (iHS) and composite likelihood ratio (CLR) methods, we identified 143 and 199 potentially selected genes in Huaxi cattle, respectively, among which nine selected genes (KCNK1, PDLIM5, CPXM2, CAPN14, MIR2285D, MYOF, PKDCC, FOXN3, and EHD3) related to ion binding, muscle growth and differentiation, and immunity were detected by both methods. Our study sheds light on the unique genetic feature and phylogenetic relationship of Huaxi cattle, provides a basis for the genetic mechanism analysis of important economic traits, and guides further intensive breeding improvement of Huaxi cattle.

Genomewide Association Analyses of Lactation Persistency and Milk Production Traits in Holstein Cattle Based on Imputed Whole-Genome Sequence Data

Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.

Long-range linkage disequilibrium in French beef cattle breeds

BACKGROUND

Linkage disequilibrium (LD) is a key parameter to study the history of populations and to identify and fine map quantitative trait loci (QTL) and it has been studied for many years in animal populations. The advent of new genotyping technologies has allowed whole-genome LD studies in most cattle populations. However, to date, long-range LD (LRLD) between distant variants on the genome has not been investigated in detail in cattle. Here, we present the first comprehensive study of LRLD in French beef cattle by analysing data on 672 Charolais (CHA), 462 Limousine (LIM) and 326 Blonde d'Aquitaine (BLA) individuals that were genotyped on the Illumina BovineHD Beadchip. Furthermore, whole-genome LD and haplotype block structure were analysed in these three breeds.

RESULTS

We computed linkage disequilibrium (r²) values for 5.9, 5.6 and 6.0 billion pairs of SNPs on the 29 autosomes of CHA, LIM and BLA, respectively. Mean r² values drop to less than 0.1 for distances between SNPs greater than 120 kb. However, for the first time, we detected the existence of LRLD in the three main French beef breeds. In total, 598, 266, and 795 LRLD events (r² ≥ 0.6) were detected in CHA, LIM and BLA, respectively. Each breed had predominantly population-specific LRLD interactions, although shared LRLD events occurred in a number of regions (55 LRLD events were shared between two breeds and nine between the three breeds). Examples of possible functional gene interactions and QTL co-location were observed with some of these LRLD events, which suggests epistatic selection.

CONCLUSIONS

We identified long-range linkage disequilibrium for the first time in French beef cattle populations. Epistatic selection may be the main source of the observed LRLD events, but other forces may also be involved. LRLD information should be accounted for in genome-wide association studies.

Use of gene expression and whole-genome sequence information to improve the accuracy of genomic prediction for carcass traits in Hanwoo cattle

Background

In this study, we assessed the accuracy of genomic prediction for carcass weight (CWT), marbling score (MS), eye muscle area (EMA) and back fat thickness (BFT) in Hanwoo cattle when using genomic best linear unbiased prediction (GBLUP), weighted GBLUP (wGBLUP), and a BayesR model. For these models, we investigated the potential gain from using pre-selected single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) on imputed sequence data and from gene expression information. We used data on 13,717 animals with carcass phenotypes and imputed sequence genotypes that were split in an independent GWAS discovery set of varying size and a remaining set for validation of prediction. Expression data were used from a Hanwoo gene expression experiment based on 45 animals.

Results

Using a larger number of animals in the reference set increased the accuracy of genomic prediction whereas a larger independent GWAS discovery dataset improved identification of predictive SNPs. Using pre-selected SNPs from GWAS in GBLUP improved accuracy of prediction by 0.02 for EMA and up to 0.05 for BFT, CWT, and MS, compared to a 50 k standard SNP array that gave accuracies of 0.50, 0.47, 0.58, and 0.47, respectively. Accuracy of prediction of BFT and CWT increased when BayesR was applied with the 50 k SNP array (0.02 and 0.03, respectively) and was further improved by combining the 50 k array with the top-SNPs (0.06 and 0.04, respectively). By contrast, using BayesR resulted in limited improvement for EMA and MS. wGBLUP did not improve accuracy but increased prediction bias. Based on the RNA-seq experiment, we identified informative expression quantitative trait loci, which, when used in GBLUP, improved the accuracy of prediction slightly, i.e. between 0.01 and 0.02. SNPs that were located in genes, the expression of which was associated with differences in trait phenotype, did not contribute to a higher prediction accuracy.

Conclusions

Our results show that, in Hanwoo beef cattle, when SNPs are pre-selected from GWAS on imputed sequence data, the accuracy of prediction improves only slightly whereas the contribution of SNPs that are selected based on gene expression is not significant. The benefit of statistical models to prioritize selected SNPs for estimating genomic breeding values is trait-specific and depends on the genetic architecture of each trait.

A Deletion Downstream of the CHCHD7 Gene Is Associated with Growth Traits in Sheep

Simple Summary

The genes CHCHD7 and PLAG1 are located on the same growth-related major quantitative trait locus of sheep. PLAG1 affects sheep growth, but no corresponding studies have been conducted on CHCHD7. However, polymorphisms in the CHCHD7 gene are associated with carcass weight and muscle formation in cattle, body height in cattle and humans, and weaning weight in Duroc pigs. In this study, the mathematical expectation method was used to analyze an 8-bp deletion mutation located downstream of the CHCHD7 gene in 2350 individuals from seven sheep breeds. The associations between wild-type and deletion genotypes and growth traits in Tan sheep were also analyzed. The 8-bp deletion locus was significantly associated with body length (p = 0.032), chest depth (p = 0.015), and chest width (p = 0.047) of Tan sheep. Additionally, wild-type genotype carriers were more numerous than those heterozygous for the deletion genotype. Thus, the genotyped 8-bp deletion downstream of the CHCHD7 gene may be associated with growth and development traits in sheep.

Abstract

In sheep, the coiled-coil-helix-coiled-coil-helix domain containing 7 (CHCHD7) gene and the pleiomorphic adenoma gene 1 (PLAG1) are on the same growth-related major quantitative trait locus, positioned head-to-head approximately 420 bp apart on chromosome 9. PLAG1 affects sheep growth, but the effects of CHCHD7 have not been determined. In this study, an 8-bp deletion downstream of CHCHD7 was analyzed in 2350 sheep from seven breeds. The associations between the deletion and growth traits of Tan sheep were also determined. Both genotypes (homozygous wild-type and heterozygous) for the 8-bp deletion were found in Tan (TS), Luxi Blackhead (LXBH), Small-Tail Han (STHS), and Lanzhou Fat-Tail (LFTS) sheep. However, there were no polymorphic sites for the mutation in Hu (HS), Sartuul (SS), and Australian White (AUW) sheep. In TS, LXBH, STHS, and LFTS sheep, the deletion genotype was less frequent than the wild-type genotype, and the allele frequencies of the deletion variant were 0.007 (TS), 0.011 (LBXH), 0.008 (STHS), and 0.010 (LFTS). The 8-bp deletion was significantly associated with body length (p = 0.032), chest depth (p = 0.015), and chest width (p = 0.047) in Tan sheep. Thus, the 8-bp deletion downstream of the CHCHD7 gene might be associated with growth and development traits of sheep.

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.

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.

Genome-wide association and pathway analysis of carcass and meat quality traits in Piemontese young bulls

A key concern in beef production is how to improve carcass and meat quality traits. Identifying the genomic regions and biological pathways that contribute to explaining variability in these traits is of great importance for selection purposes. In this study, genome wide-association (GWAS) and pathway-based analyses of carcass traits (age at slaughter (AS), carcass weight (CW), carcass daily gain (CDG), conformation score and rib-eye muscle area) and meat quality traits (pH, Warner-Bratzler shear force, purge loss, cooking loss and colour parameters (lightness, redness, yellowness, chroma, hue)) were conducted using genotype data from the 'GeneSeek Genomic Profiler Bovine LD' array in a cohort of 1166 double-muscled Piemontese beef cattle. The genome wide-association analysis was based on the GRAMMAR-GC approach and identified 37 significant single nucleotide polymorphisms (SNPs), which were associated with 12 traits (P<5 × 10-5). In particular, 14 SNPs associated with CW, CDG and AS were located at 38.57 to 38.94 Mb on Bos taurus autosome 6 and mapped within four genes, that is, Leucine Aminopeptidase 3, Family with Sequence Similarity 184 Member B, Non-SMC Condensin I Complex Subunit G and Ligand-Dependent Nuclear Receptor Corepressor-Like. Strong pairwise linkage disequilibrium was found in this region. For meat quality traits, most associations were 1 SNP per trait, except for a signal on BTA25 (at ~11.96 Mb), which was significant for four of the five meat colour parameters assessed. Gene-set enrichment analyses yielded significant results for six traits (right-sided hypergeometric test, false discovery rate <0.05). In particular, several pathways related to transmembrane transport (i.e., oxygen, calcium, ion and cation) were overrepresented for meat colour parameters. The results obtained provide useful information for genomic selection for beef production and quality in the Piemontese breed.

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.

New world goat populations are a genetically diverse reservoir for future use

Western hemisphere goats have European, African and Central Asian origins, and some local or rare breeds are reported to be adapted to their environments and economically important. By-in-large these genetic resources have not been quantified. Using 50 K SNP genotypes of 244 animals from 12 goat populations in United States, Costa Rica, Brazil and Argentina, we evaluated the genetic diversity, population structure and selective sweeps documenting goat migration to the "New World". Our findings suggest the concept of breed, particularly among "locally adapted" breeds, is not a meaningful way to characterize goat populations. The USA Spanish goats were found to be an important genetic reservoir, sharing genomic composition with the wild ancestor and with specialized breeds (e.g. Angora, Lamancha and Saanen). Results suggest goats in the Americas have substantial genetic diversity to use in selection and promote environmental adaptation or product driven specialization. These findings highlight the importance of maintaining goat conservation programs and suggest an awaiting reservoir of genetic diversity for breeding and research while simultaneously discarding concerns about breed designations.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

An intersection network based on combining SNP coassociation and RNA coexpression networks for feed utilization traits in Japanese Black cattle

Genome-wide association studies (GWAS) of quantitative traits have detected numerous genetic associations, but they encounter difficulties in pinpointing prominent candidate genes and inferring gene networks. The present study used a systems genetics approach integrating GWAS results with external RNA-expression data to detect candidate gene networks in feed utilization and growth traits of Japanese Black cattle, which are matters of concern. A SNP coassociation network was derived from significant correlations between SNPs with effects estimated by GWAS across 7 phenotypic traits. The resulting network genes contained significant numbers of annotations related to the traits. Using bovine transcriptome data from a public database, an RNA coexpression network was inferred based on the similarity of expression patterns across different tissues. An intersection network was then generated by superimposing the SNP and RNA networks and extracting shared interactions. This intersection network contained 4 tissue-specific modules: nervous system, reproductive system, muscular system, and glands. To characterize the structure (topographical properties) of the 3 networks, their scale-free properties were evaluated, which revealed that the intersection network was the most scale-free. In the subnetwork containing the most connected transcription factors (URI1, ROCK2, and ETV6), most genes were widely expressed across tissues, and genes previously shown to be involved in the traits were found. Results indicated that the current approach might be used to construct a gene network that better reflects biological information, providing encouragement for the genetic dissection of economically important quantitative traits.