Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

Genomic regions and biological pathways associated with sex-limited reproductive traits in bovine species

Many animal species exhibit sex-limited traits, where certain phenotypes are exclusively expressed in one sex. Yet, the genomic regions that contribute to these sex-limited traits in males and females remain a subject of debate. Reproductive traits are ideal phenotypes to study sexual differences since they are mostly expressed in a sex-limited way. Therefore, this study aims to use local correlation analyses to identify genomic regions and biological pathways significantly associated with male and female sex-limited traits in two distinct cattle breeds (Brahman [BB] and Tropical Composite [TC]). We used the Correlation Scan method to perform local correlation analysis on 42 trait pairs consisting of six female and seven male reproductive traits recorded on ~1,000 animals for each sex in each breed. To pinpoint a specific region associated with these sex-limited reproductive traits, we investigated the genomic region(s) consistently identified as significant across the 42 trait pairs in each breed. The genes found in the identified regions were subjected to Quantitative Trait Loci (QTL) colocalization, QTL enrichment analyses, and functional analyses to gain biological insight into sexual differences. We found that the genomic regions associated with the sex-limited reproductive phenotypes are widely distributed across all the chromosomes. However, no single region across the genome was associated with all the 42 reproductive trait pairs in the two breeds. Nevertheless, we found a region on the X-chromosome to be most significant for 80% to 90% (BB: 33 and TC: 38) of the total 42 trait pairs. A considerable number of the genes in this region were regulatory genes. By considering only genomic regions that were significant for at least 50% of the 42 trait pairs, we observed more regions spread across the autosomes and the X-chromosome. All genomic regions identified were highly enriched for trait-specific QTL linked to sex-limited traits (percentage of normal sperm, metabolic weight, average daily gain, carcass weight, age at puberty, etc.). The gene list created from these identified regions was enriched for biological pathways that contribute to the observed differences between sexes. Our results demonstrate that genomic regions associated with male and female sex-limited reproductive traits are distributed across the genome. Yet, chromosome X seems to exert a relatively larger effect on the phenotypic variation observed between the sexes.

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.

Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations

Many of the world's agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

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.

Genome-Wide Association Study to Identify QTL for Carcass Traits in Korean Hanwoo Cattle

This study aimed to identify genetic associations with carcass traits in Hanwoo cattle using a genome-wide association study. A total of 9302 phenotypes were analyzed, and all animals were genotyped using the Illumina Bovine 50K v.3 SNP chip. Heritabilities for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) were estimated as 0.42, 0.36, 0.36, and 0.47, respectively, using the GBLUP model, and 0.47, 0.37, 0.36, and 0.42, respectively, using the Bayes B model. We identified 129 common SNPs using DGEBV and 118 common SNPs using GEBV on BTA6, BTA13, and BTA14, suggesting their potential association with the traits of interest. No common SNPs were found between the GBLUP and Bayes B methods when using residuals as a response variable in GWAS. The most promising candidate genes for CWT included SLIT2, PACRGL, KCNIP4, RP1, XKR4, LYN, RPS20, MOS, FAM110B, UBXN2B, CYP7A1, SDCBP, NSMAF, TOX, CA8, LAP3, FAM184B, and NCAPG. For EMA, the genes IBSP, LAP3, FAM184B, LCORL, NCAPG, SLC30A9, and BEND4 demonstrated significance. Similarly, CYP7B1, ARMC1, PDE7A, and CRH were associated with BF, while CTSZ, GNAS, VAPB, and RAB22A were associated with MS. This finding offers valuable insights into genomic regions and molecular mechanisms influencing Hanwoo carcass traits, aiding efficient breeding strategies.

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.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Assessing Population Structure and Signatures of Selection in Wanbei Pigs Using Whole Genome Resequencing Data

Wanbei pig (WBP) is one of the indigenous pig resources in China and has many germplasm characteristics. However, research on its genome is lacking. To assess the genomic variation, population structure, and selection signatures, we resequenced 18 WBP for the first time and performed a comprehensive analysis with resequenced data of 10 Asian wild boars. In total, 590.03 Gb of data and approximately 41 million variants were obtained. Polymorphism level (θπ) ratio and genetic differentiation (fixation index)-based cross approaches were applied, and 539 regions, which harbored 176 genes, were selected. Functional analysis of the selected genes revealed that they were associated with lipid metabolism (SCP2, APOA1, APOA4, APOC3, CD36, BCL6, ADCY8), backfat thickness (PLAG1, CACNA2D1), muscle (MYOG), and reproduction (CABS1). Overall, our results provide a valuable resource for characterizing the uniqueness of WBP and a basis for future breeding.

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.

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.

Proteomic Analysis of Hypothalamus and Pituitary Gland in Pre and Postpubertal Brahman Heifers

The hypothalamus and the pituitary gland are directly involved in the complex systemic changes that drive the onset of puberty in cattle. Here, we applied integrated bioinformatics to elucidate the critical proteins underlying puberty and uncover potential molecular mechanisms from the hypothalamus and pituitary gland of prepubertal (n = 6) and postpubertal (n = 6) cattle. Proteomic analysis in the hypothalamus and pituitary gland revealed 275 and 186 differentially abundant (DA) proteins, respectively (adjusted p-value < 0.01). The proteome profiles found herein were integrated with previously acquired transcriptome profiles. These transcriptomic studies used the same tissues harvested from the same heifers at pre- and post-puberty. This comparison detected a small number of matched transcripts and protein changes at puberty in each tissue, suggesting the need for multiple omics analyses for interpreting complex biological systems. In the hypothalamus, upregulated DA proteins at post-puberty were enriched in pathways related to puberty, including GnRH, calcium and oxytocin signalling pathways, whereas downregulated proteins were observed in the estrogen signalling pathway, axon guidance and GABAergic synapse. Additionally, this study revealed that ribosomal pathway proteins in the pituitary were involved in the pubertal development of mammals. The reported molecules and derived protein-protein networks are a starting point for future experimental approaches that might dissect with more detail the role of each molecule to provide new insights into the mechanisms of puberty onset in cattle.

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.

Sustainable Intensification of Beef Production in the Tropics: The Role of Genetically Improving Sexual Precocity of Heifers

Simple Summary

Tropical pasture-based beef production systems play a vital role in global food security. The importance of promoting sustainable intensification of such systems has been debated worldwide. Demand for beef is growing together with concerns over the impact of its production on the environment. Implementing sustainable livestock intensification programs relies on animal genetic improvement. In tropical areas, the lack of sexual precocity is a bottleneck for cattle efficiency, directly impacting the sustainability of production systems. In the present review we present and discuss the state of the art of genetic evaluation for sexual precocity in Bos indicus beef cattle, covering the definition of measurable traits, genetic parameter estimates, genomic analyses, and a case study of selection for sexual precocity in Nellore breeding programs.

Abstract

Increasing productivity through continued animal genetic improvement is a crucial part of implementing sustainable livestock intensification programs. In Zebu cattle, the lack of sexual precocity is one of the main obstacles to improving beef production efficiency. Puberty-related traits are complex, but large-scale data sets from different “omics” have provided information on specific genes and biological processes with major effects on the expression of such traits, which can greatly increase animal genetic evaluation. In addition, genetic parameter estimates and genomic predictions involving sexual precocity indicator traits and productive, reproductive, and feed-efficiency related traits highlighted the feasibility and importance of direct selection for anticipating heifer reproductive life. Indeed, the case study of selection for sexual precocity in Nellore breeding programs presented here show that, in 12 years of selection for female early precocity and improved management practices, the phenotypic means of age at first calving showed a strong decreasing trend, changing from nearly 34 to less than 28 months, with a genetic trend of almost −2 days/year. In this period, the percentage of early pregnancy in the herds changed from around 10% to more than 60%, showing that the genetic improvement of heifer’s sexual precocity allows optimizing the productive cycle by reducing the number of unproductive animals in the herd. It has a direct impact on sustainability by better use of resources. Genomic selection breeding programs accounting for genotype by environment interaction represent promising tools for accelerating genetic progress for sexual precocity in tropical beef cattle.

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.

Bos taurus haplotypes segregating in Nellore (Bos indicus) cattle

Brazil is the largest exporter of beef in the world, and most of that beef derives from Nellore cattle. Although considered a zebu breed (Bos indicus), the history of Nellore cattle in Brazil is marked by the importation of bulls from India, the use of a Creole taurine (Bos taurus) maternal lineage to quickly expand the herds and backcrossing to Nellore bulls to recover zebu ancestry. As a consequence, the current Brazilian Nellore population carries an average taurine ancestry of approximately 1%. Although that percentage seems small, some taurine variants deviate substantially from that average, with the better-known cases being the PLAG1-Q haplotype involved with body size variation and the Guarani (P_G ) polled variant producing hornless animals. Here, we report taurine haplotypes in 9074 Nellore animals genotyped for 539 657 imputed SNP markers. Apart from PLAG1-Q and P_G , our analysis further revealed common taurine haplotypes (>3%) spanning genes related to immunity, growth, reproduction and hair and skin phenotypes. Using data from 22 economically important traits, we showed that many of the major QTL previously reported in the breed are at least partially driven by taurine haplotypes. As B. taurus and B. indicus haplotypes are highly divergent, presenting widely different sets of functional variants, our results provide promising targets for future scrutiny in Nellore cattle.

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.

A high plane of nutrition during early life alters the hypothalamic transcriptome of heifer calves

The aim was to examine the effect of rapid body weight gain during early calfhood consistent with earlier sexual development on the transcriptional profile of the hypothalamus. Angus X Holstein-Friesian heifer calves (19 ± 5 days of age) were offered a high (HI, n = 14) or moderate (MOD, n = 15) plane of nutrition from 3 to 21 weeks of age to achieve a growth rate of 1.2 kg/d and 0.5 kg/d, respectively. Following euthanasia at 21 weeks, the arcuate nucleus (ARC) region was separated from the remainder of the hypothalamus and both were subjected to RNA-Seq. HI calves exhibited altered expression of 80 and 39 transcripts in the ARC and the remaining hypothalamus, respectively (P < 0.05) including downregulation of AGRP and NPY and upregulation of POMC, previously implicated in precocious sexual development. Stress-signaling pathways were amongst the most highly dysregulated. Organ morphology, reproductive system development and function, and developmental disorder were amongst the networks derived from differentially expressed genes (DEGs) in the ARC. Gene co-expression analysis revealed DEGs within the ARC (POMC, CBLN2, CHGA) and hypothalamus (PENK) as hub genes. In conclusion, enhanced nutrition during early calfhood alters the biochemical regulation of the hypothalamus consistent with advanced sexual development in the prepubertal heifer.

Genome-Wide Association Analysis of Growth Curve Parameters in Chinese Simmental Beef Cattle

The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight-age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R² = 0.954). The parameters' mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.

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.

Genome-wide association study identified genomic regions and putative candidate genes affecting meat color traits in Nellore cattle

Single and multiple-trait GWAS were conducted to detect genomic regions and candidate genes associated with meat color traits (L*, lightness; a*, redness; b*, yellowness) in Nellore cattle. Phenotypic records of 5000 animals, and 3794 genotypes for 614,274 SNPs were used. The BLUPF90 family programs were used through single step GWAS approach. The top 10 genomic regions from single-trait GWAS explained 13.64%, 15.12% and 13% of genetic variance of L*, a* and b*, which harbored 129, 70, and 84 candidate genes, respectively. Regarding multiple-trait GWAS, the top 10 SNP windows explained 17.46%, 18.98% and 13.74% of genetic variance of L*, a* and b*, and harbored 124, 86, and 82 candidate genes, respectively. Pleiotropic effects were evidenced by the overlapping regions detected on BTA 15 and 26 associated with L* and a* (genetic correlation of -0.53), and on BTA 18 associated with a* and b* (genetic correlation of 0.60). Similar genomic regions located on BTA 2, 5, 6, and 18 were detected through single and multi-trait GWAS. Overlapped regions harbored a total of 30 functional candidate genes involved in mitochondrial activity, structural integrity of muscles, lipid oxidation, anaerobic metabolism, and muscular pH.

Association analysis of loci implied in "buffering" epistasis

The existence of buffering mechanisms is an emerging property of biological networks, and this results in the buildup of robustness through evolution. So far, there are no explicit methods to find loci implied in buffering mechanisms. However, buffering can be seen as interaction with genetic background. Here we develop this idea into a tractable model for quantitative genetics, in which the buffering effect of one locus with many other loci is condensed into a single statistical effect, multiplicative on the total additive genetic effect. This allows easier interpretation of the results and simplifies the problem of detecting epistasis from quadratic to linear in the number of loci. Using this formulation, we construct a linear model for genome-wide association studies that estimates and declares the significance of multiplicative epistatic effects at single loci. The model has the form of a variance components, norm reaction model and likelihood ratio tests are used for significance. This model is a generalization and explanation of previous ones. We test our model using bovine data: Brahman and Tropical Composite animals, phenotyped for body weight at yearling and genotyped at high density. After association analysis, we find a number of loci with buffering action in one, the other, or both breeds; these loci do not have a significant statistical additive effect. Most of these loci have been reported in previous studies, either with an additive effect or as footprints of selection. We identify buffering epistatic SNPs present in or near genes reported in the context of signatures of selection in multi-breed cattle population studies. Prominent among these genes are those associated with fertility (INHBA, TSHR, ESRRG, PRLR, and PPARG), growth (MSTN, GHR), coat characteristics (KIT, MITF, PRLR), and heat resistance (HSPA6 and HSPA1A). In these populations, we found loci that have a nonsignificant statistical additive effect but a significant epistatic effect. We argue that the discovery and study of loci associated with buffering effects allow attacking the difficult problems, among others, of the release of maintenance variance in artificial and natural selection, of quick adaptation to the environment, and of opposite signs of marker effects in different backgrounds. We conclude that our method and our results generate promising new perspectives for research in evolutionary and quantitative genetics based on the study of loci that buffer effect of other loci.

Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color

Coat color is among the most distinctive phenotypes in cattle. Worldwide, several breeds share peculiar coat color features such as the presence of a fawn pigmentation of the calf at birth, turning over time to grey, and sexual dichromatism. The aim of this study was to search for polymorphisms under differential selection by contrasting grey cattle breeds displaying the above phenotype with non-grey cattle breeds, and to identify the underlying genes. Using medium-density SNP array genotype data, a multi-cohort F_ST-outlier approach was adopted for a total of 60 pair-wise comparisons of the 15 grey with 4 non-grey cattle breeds (Angus, Limousin, Charolais, and Holstein), with the latter selected as representative of solid and piebald phenotypes, respectively. Overall, more than 50 candidate genes were detected; almost all were either directly or indirectly involved in pigmentation, and some of them were already known for their role in phenotypes related with hair graying in mammals. Notably, 17 relevant genes, including SDR16C5, MOS, SDCBP, and NSMAF, were located in a signal on BTA14 convergently observed in all the four considered scenarios. Overall, the key stages of pigmentation (melanocyte development, melanogenesis, and pigment trafficking/transfer) were all represented among the pleiotropic functions of the candidate genes, suggesting the complex nature of the grey phenotype in cattle.

Estimation of Genomic Breed Composition for Purebred and Crossbred Animals Using Sparsely Regularized Admixture Models

A variety of statistical methods, such as admixture models, have been used to estimate genomic breed composition (GBC). These methods, however, tend to produce non-zero components to reference breeds that shared some genomic similarity with a test animal. These non-essential GBC components, in turn, offset the estimated GBC for the breed to which it belongs. As a result, not all purebred animals have 100% GBC of their respective breeds, which statistically indicates an elevated false-negative rate in the identification of purebred animals with 100% GBC as the cutoff. Otherwise, a lower cutoff of estimated GBC will have to be used, which is arbitrary, and the results are less interpretable. In the present study, three admixture models with regularization were proposed, which produced sparse solutions through suppressing the noise in the estimated GBC due to genomic similarities. The regularization or penalty forms included the L1 norm penalty, minimax concave penalty (MCP), and smooth clipped absolute deviation (SCAD). The performances of these regularized admixture models on the estimation of GBC were examined in purebred and composite animals, respectively, and compared to that of the non-regularized admixture model as the baseline model. The results showed that, given optimal values for λ, the three sparsely regularized admixture models had higher power and thus reduced the false-negative rate for the breed identification of purebred animals than the non-regularized admixture model. Of the three regularized admixture models, the two with a non-convex penalty outperformed the one with L1 norm penalty. In the Brangus, a composite cattle breed, estimated GBC were roughly comparable among the four admixture models, but all the four models underestimated the GBC for these composite animals when non-ancestral breeds were included as the reference. In conclusion, the admixture models with sparse regularization gave more parsimonious, consistent and interpretable results of estimated GBC for purebred animals than the non-regularized admixture model. Nevertheless, the utility of regularized admixture models for estimating GBC in crossbred or composite animals needs to be taken with caution.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

Use of whole-genome sequence data and novel genomic selection strategies to improve selection for age at puberty in tropically-adapted beef heifers

Background

In tropically-adapted beef heifers, application of genomic prediction for age at puberty has been limited due to low prediction accuracies. Our aim was to investigate novel methods of pre-selecting whole-genome sequence (WGS) variants and alternative analysis methodologies; including genomic best linear unbiased prediction (GBLUP) with multiple genomic relationship matrices (MGRM) and Bayesian (BayesR) analyses, to determine if prediction accuracy for age at puberty can be improved.

Methods

Genotypes and phenotypes were obtained from two research herds. In total, 868 Brahman and 960 Tropical Composite heifers were recorded in the first population and 3695 Brahman, Santa Gertrudis and Droughtmaster heifers were recorded in the second population. Genotypes were imputed to 23 million whole-genome sequence variants. Eight strategies were used to pre-select variants from genome-wide association study (GWAS) results using conditional or joint (COJO) analyses. Pre-selected variants were included in three models, GBLUP with a single genomic relationship matrix (SGRM), GBLUP MGRM and BayesR. Five-way cross-validation was used to test the effect of marker panel density (6 K, 50 K and 800 K), analysis model, and inclusion of pre-selected WGS variants on prediction accuracy.

Results

In all tested scenarios, prediction accuracies for age at puberty were highest in BayesR analyses. The addition of pre-selected WGS variants had little effect on the accuracy of prediction when BayesR was used. The inclusion of WGS variants that were pre-selected using a meta-analysis with COJO analyses by chromosome, fitted in a MGRM model, had the highest prediction accuracies in the GBLUP analyses, regardless of marker density. When the low-density (6 K) panel was used, the prediction accuracy of GBLUP was equal (0.42) to that with the high-density panel when only six additional sequence variants (identified using meta-analysis COJO by chromosome) were included.

Conclusions

While BayesR consistently outperforms other methods in terms of prediction accuracies, reasonable improvements in accuracy can be achieved when using GBLUP and low-density panels with the inclusion of a relatively small number of highly relevant WGS variants.

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.

Reaffirmation of known major genes and the identification of novel candidate genes associated with carcass-related metrics based on whole genome sequence within a large multi-breed cattle population

Background

The high narrow sense heritability of carcass traits suggests that the underlying additive genetic potential of an individual should be strongly correlated with both animal carcass quality and quantity, and therefore, by extension, carcass value. Therefore, the objective of the present study was to detect genomic regions associated with three carcass traits, namely carcass weight, conformation and fat cover, using imputed whole genome sequence in 28,470 dairy and beef sires from six breeds with a total of 2,199,926 phenotyped progeny.

Results

Major genes previously associated with carcass performance were identified, as well as several putative novel candidate genes that likely operate both within and across breeds. The role of MSTN in carcass performance was re-affirmed with the segregating Q204X mutation explaining 1.21, 1.11 and 5.95% of the genetic variance in carcass weight, fat and conformation, respectively in the Charolais population. In addition, a genomic region on BTA6 encompassing the NCAPG/LCORL locus, which is a known candidate locus associated with body size, was associated with carcass weight in Angus, Charolais and Limousin. Novel candidate genes identified included ZFAT in Angus, and SLC40A1 and the olfactory gene cluster on BTA15 in Charolais. Although the majority of associations were breed specific, associations that operated across breeds included SORCS1 on BTA26, MCTP2 on BTA21 and ARL15 on BTA20; these are of particular interest due to their potential informativeness in across-breed genomic evaluations. Genomic regions affecting all three carcass traits were identified in each of the breeds, although these were mainly concentrated on BTA2 and BTA6, surrounding MSTN and NCAPG/LCORL, respectively. This suggests that although major genes may be associated with all three carcass traits, the majority of genes containing significant variants (unadjusted p-value < 10^− 4) may be trait specific associations of small effect.

Conclusions

Although plausible novel candidate genes were identified, the proportion of variance explained by these candidates was minimal thus reaffirming that while carcass performance may be affected by major genes in the form of MSTN and NCAPG/LCORL, the majority of variance is attributed to the additive (and possibly multiplicative) effect of many polymorphisms of small effect.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6071-9) contains supplementary material, which is available to authorized users.

Genomic clues of the evolutionary history of Bos indicus cattle

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro-Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian-Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole-genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.

Across-breed validation study confirms and identifies new loci associated with sexual precocity in Brahman and Nellore cattle

The aim of this study was to identify candidate regions associated with sexual precocity in Bos indicus. Nellore and Brahman were set as validation and discovery populations, respectively. SNP selected in Brahman to validate in Nellore were from gene regions affecting reproductive traits (G1) and significant SNP (p ≤ 10^-3 ) from a meta-analysis (G2). In the validation population, early pregnancy (EP) and scrotal circumference (SC) were evaluated. To perform GWAS in validation population, we used regression and Bayes C. SNP with p ≤ 10^-3 in regression and Bayes factor ≥3 in Bayes C were deemed significant. Significant SNP (for EP or SC) or SNP in their ±250 Kb vicinity region, which were in at least one discovery set (G1 or G2), were considered validated. SNP identified in both G1 and G2 were considered candidate. For EP, 145 SNP were validated in G1 and 41 in G2, and for SC, these numbers were 14 and 2. For EP, 21 candidate SNP were detected (G1 and G2). For SC, no candidate SNP were identified. Validated SNP and their vicinity region were located close to quantitative trait loci or genes related to reproductive traits and were enriched in gene ontology terms related to reproductive success. These are therefore strong candidate regions for sexual precocity in Nellore and Brahman.

Genome-wide association study reveals candidate genes associated with body measurement traits in Chinese Wagyu beef cattle

We performed a genome-wide association study to identify candidate genes for body measurement traits in 463 Wagyu beef cattle typed with the Illumina Bovine HD 770K SNP array. At the genome-wide level, we detected 18, five and one SNPs associated with hip height, body height and body length respectively. In total, these SNPs are within or near 11 genes, six of which (PENK, XKR4, IMPAD1, PLAG1, CCND2 and SNTG1) have been reported previously and five of which (CSMD3, LAP3, SYN3, FAM19A5 and TIMP3) are novel candidate genes that we found to be associated with body measurement traits. Further exploration of these candidate genes will facilitate genetic improvement in Chinese Wagyu beef cattle.

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.

Towards multi-breed genomic evaluations for female fertility of tropical beef cattle

Developing accurate genomic evaluations of fertility for tropical beef cattle must deal with at least two major challenges (i) recording cow fertility traits in extensive production systems on large numbers of cows and (ii) the genomic evaluations should work across the breeds, crossbreds, and composites used in tropical beef production. Here, we assess accuracy of genomic evaluations for a trait which can be collected on a large scale in extensive conditions, corpus luteum score (CLscore), which is 1 if ovarian scanning indicates a heifer has cycled by 600 d and 0 if not, in a multi-breed population. A total of 3,696 heifers, including 979 Brahmans, 914 Droughtmasters, and 1,803 Santa Gertrudis in seven herds across 3-yr cohorts with CLscores, were genotyped for 24,211 SNPs. Genotypes were imputed to 728,785 SNPs. GBLUP and BayesR were used to predict GEBV. Accuracy of GEBV was evaluated with two validation strategies. In the first strategy, the last year cohort of heifers from each herd was used for validation, such that every herd had heifers in both reference and validation populations. In the second validation strategy, each herd in turn was removed in its entirety from the reference population, and was used for validation. For both validation strategies, accuracy of GEBV for single breed and multi-breed reference populations was assessed. For the first validation strategy, accuracy of GEBV ranged from 0.2 for Brahmans to 0.4 for Droughtmasters. Increasing marker density from 24K SNPs to 728K SNPs resulted in a small increase in accuracy, and including multiple-breeds in the reference did not help improve accuracy. These results suggest that provided a herd has animals in the reference population, the accuracy of the GEBV is largely determined by within herd (linkage) information. The situation was very different when entire herds were predicted in the second validation. In this case accuracy of GEBV using only 24K SNPs and only a within breed reference was close to zero for all breeds. Accuracy increased substantially when 728K SNPs, BayesR, and a multi-breed reference were used, from 0.15 for Brahmans to 0.35 for Santa Gertrudis. Given the second validation strategy is more likely to reflect the situation for many herds in tropical beef production (no animals in the reference), genomic evaluations for fertility in tropical beef cattle should be based on high-density markers (728K SNPs) and should be multi-breed.

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.

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.

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.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals

Stature is affected by many polymorphisms of small effect in humans ¹ . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes^2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10^-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

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.

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.

Genome-wide association study provides strong evidence of genes affecting the reproductive performance of Nellore beef cows

Reproductive traits are economically important for beef cattle production; however, these traits are still a bottleneck in indicine cattle since these animals typically reach puberty at older ages when compared to taurine breeds. In addition, reproductive traits are complex phenotypes, i.e., they are controlled by both the environment and many small-effect genes involved in different pathways. In this study, we conducted genome-wide association study (GWAS) and functional analyses to identify important genes and pathways associated with heifer rebreeding (HR) and with the number of calvings at 53 months of age (NC53) in Nellore cows. A total of 142,878 and 244,311 phenotypes for HR and NC53, respectively, and 2,925 animals genotyped with the Illumina Bovine HD panel (Illumina^®, San Diego, CA, USA) were used in GWAS applying the weighted single-step GBLUP (WssGBLUP) method. Several genes associated with reproductive events were detected in the 20 most important 1Mb windows for both traits. Significant pathways for HR and NC53 were associated with lipid metabolism and immune processes, respectively. MHC class II genes, detected on chromosome 23 (window 25-26Mb) for NC53, were significantly associated with pregnancy success of Nellore cows. These genes have been proved previously to be associated with reproductive traits such as mate choice in other breeds and species. Our results suggest that genes associated with the reproductive traits HR and NC53 may be involved in embryo development in mammalian species. Furthermore, some genes associated with mate choice may affect pregnancy success in Nellore cattle.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: , , , , and . These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy.

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.

Association between single nucleotide polymorphisms and sexual precocity in Nellore heifers

The aim of this study was to determine the extent (r²) of linkage disequilibrium (LD) in the genome of Nellore cattle, and to examine associations between single nucleotide polymorphisms (SNP) and age at first calving (AFC) and early pregnancy (EP) using a panel of high-density SNPs and data from 1182 Nellore females. A total of 13 contemporary groups (CG) were used consisting of farm, season, and year of birth. For genome-wide association analysis, SNPs with a minor allele frequency (MAF)<0.05 and animals with a call rate<0.90 were excluded, totaling 431,885 SNPs. For statistical analysis, a linear model was used for AFC and a threshold model for EP. To estimate the significance of the associations for the two traits, the model included the categorical fixed effects of CG, SNPs, and sire. In addition, the polygenic effect was included in the analysis. The additive effects and dominance deviations of Bonferroni-adjusted significant SNPs for AFC and EP were estimated using orthogonal contrasts. The average estimate of r² for all autosomes was 0.18 at a distance of 4.8kb and the mean MAF was 0.25±0.13. The LD decreased as the distance between markers increased: 0.35 (1kb) to 0.12 (100kb). Eleven significant associations were detected in seven different chromosomes. Seven SNPs were associated with AFC and four were associated with EP. Three SNPs were significant for both traits. The identification of SNPs associated with AFC and EP may contribute for selecting sexually precocious animals.