A multi-trait, meta-analysis for detecting pleiotropic polymorphisms for stature, fatness and reproduction in beef cattle

Comparative analysis of inbreeding parameters and runs of homozygosity islands in 2 Italian autochthonous cattle breeds mainly raised in the Parmigiano-Reggiano cheese production region

Reggiana and Modenese are autochthonous cattle breeds, reared in the North of Italy, that can be mainly distinguished for their standard coat color (Reggiana is red, whereas Modenese is white with some pale gray shades). Almost all milk produced by these breeds is transformed into 2 mono-breed branded Parmigiano-Reggiano cheeses, from which farmers receive the economic incomes needed for the sustainable conservation of these animal genetic resources. After the setting up of their herd books in 1960s, these breeds experienced a strong reduction in the population size that was subsequently reverted starting in the 1990s (Reggiana) or more recently (Modenese) reaching at present a total of about 2,800 and 500 registered cows, respectively. Due to the small population size of these breeds, inbreeding is a very important cause of concern for their conservation programs. Inbreeding is traditionally estimated using pedigree data, which are summarized in an inbreeding coefficient calculated at the individual level (F_PED). However, incompleteness of pedigree information and registration errors can affect the effectiveness of conservation strategies. High-throughput SNP genotyping platforms allow investigation of inbreeding using genome information that can overcome the limits of pedigree data. Several approaches have been proposed to estimate genomic inbreeding, with the use of runs of homozygosity (ROH) considered to be the more appropriate. In this study, several pedigree and genomic inbreeding parameters, calculated using the whole herd book populations or considering genotyping information (GeneSeek GGP Bovine 150K) from 1,684 Reggiana cattle and 323 Modenese cattle, were compared. Average inbreeding values per year were used to calculate effective population size. Reggiana breed had generally lower genomic inbreeding values than Modenese breed. The low correlation between pedigree-based and genomic-based parameters (ranging from 0.187 to 0.195 and 0.319 to 0.323 in the Reggiana and Modenese breeds, respectively) reflected the common problems of local populations in which pedigree records are not complete. The high proportion of short ROH over the total number of ROH indicates no major recent inbreeding events in both breeds. ROH islands spread over the genome of the 2 breeds (15 in Reggiana and 14 in Modenese) identified several signatures of selection. Some of these included genes affecting milk production traits, stature, body conformation traits (with a main ROH island in both breeds on BTA6 containing the ABCG2, NCAPG, and LCORL genes) and coat color (on BTA13 in Modenese containing the ASIP gene). In conclusion, this work provides an extensive comparative analysis of pedigree and genomic inbreeding parameters and relevant genomic information that will be useful in the conservation strategies of these 2 iconic local cattle breeds.

Direct introgression of untapped diversity into elite wheat lines

The effective utilization of natural variation has become essential in addressing the challenges that climate change and population growth pose to global food security. Currently adopted protracted approaches to introgress exotic alleles into elite cultivars need substantial transformation. Here, through a strategic three-way crossing scheme among diverse exotics and the best historical elites (exotic/elite1//elite2), 2,867 pre-breeding lines were developed, genotyped and screened for multiple agronomic traits in four mega-environments. A meta-genome-wide association study, selective sweeps and haplotype-block-based analyses unveiled selection footprints in the genomes of pre-breeding lines as well as exotic-specific associations with agronomic traits. A simulation with a neutrality assumption demonstrated that many pre-breeding lines had significant exotic contributions despite substantial selection bias towards elite genomes. National breeding programmes worldwide have adopted 95 lines for germplasm enhancement, and 7 additional lines are being advanced in varietal release trials. This study presents a great leap forwards in the mobilization of GenBank variation to the breeding pipelines.

New loci and neuronal pathways for resilience to heat stress in cattle

While understanding the genetic basis of heat tolerance is crucial in the context of global warming's effect on humans, livestock, and wildlife, the specific genetic variants and biological features that confer thermotolerance in animals are still not well characterized. We used dairy cows as a model to study heat tolerance because they are lactating, and therefore often prone to thermal stress. The data comprised almost 0.5 million milk records (milk, fat, and proteins) of 29,107 Australian Holsteins, each having around 15 million imputed sequence variants. Dairy animals often reduce their milk production when temperature and humidity rise; thus, the phenotypes used to measure an individual's heat tolerance were defined as the rate of milk production decline (slope traits) with a rising temperature-humidity index. With these slope traits, we performed a genome-wide association study (GWAS) using different approaches, including conditional analyses, to correct for the relationship between heat tolerance and level of milk production. The results revealed multiple novel loci for heat tolerance, including 61 potential functional variants at sites highly conserved across 100 vertebrate species. Moreover, it was interesting that specific candidate variants and genes are related to the neuronal system (ITPR1, ITPR2, and GRIA4) and neuroactive ligand-receptor interaction functions for heat tolerance (NPFFR2, CALCR, and GHR), providing a novel insight that can help to develop genetic and management approaches to combat heat stress.

Genome-wide association studies for production, respiratory disease, and immune-related traits in Landrace pigs

Identification of a quantitative trait locus (QTL) related to a chronic respiratory disease such as Mycoplasmal pneumonia of swine (MPS) and immune-related traits is important for the genetic improvement of disease resistance in pigs. The objective of this study was to detect a novel QTL for a total of 22 production, respiratory disease, and immune-related traits in Landrace pigs. A total of 874 Landrace purebred pigs, which were selected based on MPS resistance, were genotyped using the Illumina PorcineSNP60 BeadChip. We performed single nucleotide polymorphism (SNP)-based and haplotype-based genome-wide association studies (GWAS) to detect a novel QTL and to evaluate the possibility of a pleiotropic QTL for these traits. SNP-based GWAS detected a total of six significant regions in backfat thickness, ratio of granular leucocytes to lymphatic cells, plasma concentration of cortisol at different ages, and complement alternative pathway activity in serum. The significant region detected by haplotype-based GWAS was overlapped across the region detected by SNP-based GWAS. Most of these detected QTL regions were novel regions with some candidate genes located in them. With regard to a pleiotropic QTL among traits, only three of these detected QTL regions overlapped among traits, and many detected regions independently affected the traits.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Meta-GWAS for quantitative trait loci identification in soybean

We report a meta-Genome Wide Association Study involving 73 published studies in soybean [Glycine max L. (Merr.)] covering 17,556 unique accessions, with improved statistical power for robust detection of loci associated with a broad range of traits. De novo GWAS and meta-analysis were conducted for composition traits including fatty acid and amino acid composition traits, disease resistance traits, and agronomic traits including seed yield, plant height, stem lodging, seed weight, seed mottling, seed quality, flowering timing, and pod shattering. To examine differences in detectability and test statistical power between single- and multi-environment GWAS, comparison of meta-GWAS results to those from the constituent experiments were performed. Using meta-GWAS analysis and the analysis of individual studies, we report 483 peaks at 393 unique loci. Using stringent criteria to detect significant marker-trait associations, 59 candidate genes were identified, including 17 agronomic traits loci, 19 for seed-related traits, and 33 for disease reaction traits. This study identified potentially valuable candidate genes that affect multiple traits. The success in narrowing down the genomic region for some loci through overlapping mapping results of multiple studies is a promising avenue for community-based studies and plant breeding applications.

A conditional multi-trait sequence GWAS discovers pleiotropic candidate genes and variants for sheep wool, skin wrinkle and breech cover traits

Background

Imputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to “fly strike” (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle.

Results

A stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed.

Conclusions

The wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00651-0.

Meta-analysis of genome-wide association studies reveal common loci controlling agronomic and quality traits in a wide range of normal and heat stressed environments

Several stable QTL were detected using metaGWAS analysis for different agronomic and quality traits under 26 normal and heat stressed environments. Heat stress, exacerbated by global warming, has a negative influence on wheat production worldwide and climate resilient cultivars can help mitigate these impacts. Selection decisions should therefore depend on multi-environment experiments representing a range of temperatures at critical stages of development. Here, we applied a meta-genome wide association analysis (metaGWAS) approach to detect stable QTL with significant effects across multiple environments. The metaGWAS was applied to 11 traits scored in 26 trials that were sown at optimal or late times of sowing (TOS1 and TOS2, respectively) at five locations. A total of 2571 unique wheat genotypes (13,959 genotypes across all environments) were included and the analysis conducted on TOS1, TOS2 and both times of sowing combined (TOS1&2). The germplasm was genotyped using a 90 k Infinium chip and imputed to exome sequence level, resulting in 341,195 single nucleotide polymorphisms (SNPs). The average accuracy across all imputed SNPs was high (92.4%). The three metaGWAS analyses revealed 107 QTL for the 11 traits, of which 16 were detected in all three analyses and 23 were detected in TOS1&2 only. The remaining QTL were detected in either TOS1 or TOS2 with or without TOS1&2, reflecting the complex interactions between the environments and the detected QTL. Eight QTL were associated with grain yield and seven with multiple traits. The identified QTL provide an important resource for gene enrichment and fine mapping to further understand the mechanisms of gene × environment interaction under both heat stressed and unstressed conditions.

Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono‐breed Parmigiano–Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome‐wise level the within‐breed detected variability of three pigmentation‐related traits (intensity of red coat colour, based on three classes – light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour – pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two‐thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single‐marker and haplotype‐based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG–LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.

Detection of Genomic Regions with Pleiotropic Effects for Growth and Carcass Quality Traits in the Rubia Gallega Cattle Breed

Simple Summary

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. We have developed a ssGWAS by backsolving the SNP effects after implementing a ssGBLUP. The results showed an apparent heterogeneity of the additive genetic variance across the genome. Some of the genomic regions explaining the most of this additive variance were shared across traits, indicating the presence of pleiotropic effects, which were reflected in their genetic correlations.

Abstract

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. In recent years, genomic evaluation has been implemented by using a ssGBLUP (single-step Genomic Best Linear Unbiased Prediction). This procedure can reparameterized to perform ssGWAS (single-step Genome Wide Association Studies) by backsolving the SNP (single nucleotide polymorphisms) effects. Therefore, the objective of this study was to identify genomic regions associated with the genetic variability in growth and carcass quality traits. We implemented a ssGBLUP by using a database that included records for Birth Weight (BW-327,350 records-), Weaning Weight (WW-83,818-), Cold Carcass Weight (CCW-91,621-), Fatness (FAT-91,475-) and Conformation (CON-91,609-). The pedigree included 464,373 individuals, 2449 of which were genotyped. After a process of filtering, we ended up using 43,211 SNP markers. We used the GBLUP and SNPBLUP model equivalences to obtain the effects of the SNPs and then calculated the percentage of variance explained by the regions of the genome between 1 Mb. We identified 7 regions of the genome for CCW; 8 regions for BW, WW, FAT and 9 regions for CON, which explained the percentage of variance above 0.5%. Furthermore, a number of the genome regions had pleiotropic effects, located at: BTA1 (131–132 Mb), BTA2 (1–11 Mb), BTA3 (32–33 Mb), BTA6 (36–38 Mb), BTA16 (24–26 Mb), and BTA 21 (56–57 Mb). These regions contain, amongst others, the following candidate genes: NCK1, MSTN, KCNA3, LCORL, NCAPG, and RIN3.

Genomic selection and genetic gain for nut yield in an Australian macadamia breeding population

Background

Improving yield prediction and selection efficiency is critical for tree breeding. This is vital for macadamia trees with the time from crossing to production of new cultivars being almost a quarter of a century. Genomic selection (GS) is a useful tool in plant breeding, particularly with perennial trees, contributing to an increased rate of genetic gain and reducing the length of the breeding cycle. We investigated the potential of using GS methods to increase genetic gain and accelerate selection efficiency in the Australian macadamia breeding program with comparison to traditional breeding methods. This study evaluated the prediction accuracy of GS in a macadamia breeding population of 295 full-sib progeny from 32 families (29 parents, reciprocals combined), along with a subset of parents. Historical yield data for tree ages 5 to 8years were used in the study, along with a set of 4113 SNP markers. The traits of focus were average nut yield from tree ages 5 to 8years and yield stability, measured as the standard deviation of yield over these 4 years. GBLUP GS models were used to obtain genomic estimated breeding values for each genotype, with a five-fold cross-validation method and two techniques: prediction across related populations and prediction across unrelated populations.

Results

Narrow-sense heritability of yield and yield stability was low (h²=0.30 and 0.04, respectively). Prediction accuracy for yield was 0.57 for predictions across related populations and 0.14 when predicted across unrelated populations. Accuracy of prediction of yield stability was high (r=0.79) for predictions across related populations. Predicted genetic gain of yield using GS in related populations was 474g/year, more than double that of traditional breeding methods (226g/year), due to the halving of generation length from 8 to 4years.

Conclusions

The results of this study indicate that the incorporation of GS for yield into the Australian macadamia breeding program may accelerate genetic gain due to reduction in generation length, though the cost of genotyping appears to be a constraint at present.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07694-z.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Multivariate Analysis Identifies Eight Novel Loci Associated with Meat Productivity Traits in Sheep

Despite their economic value, sheep remain relatively poorly studied animals in terms of the number of known loci and genes associated with commercially important traits. This gap in our knowledge can be filled in by performing new genome-wide association studies (GWAS) or by re-analyzing previously documented data using novel powerful statistical methods. This study is focused on the search for new loci associated with meat productivity and carcass traits in sheep. With a multivariate approach applied to publicly available GWAS results, we identified eight novel loci associated with the meat productivity and carcass traits in sheep. Using an in silico follow-up approach, we prioritized 13 genes in these loci. One of eight novel loci near the FAM3C and WNT16 genes has been replicated in an independent sample of Russian sheep populations (N = 108). The novel loci were added to our regularly updated database increasing the number of known loci to more than 140.

A gene co-association network regulating gut microbial communities in a Duroc pig population

BACKGROUND

Analyses of gut microbiome composition in livestock species have shown its potential to contribute to the regulation of complex phenotypes. However, little is known about the host genetic control over the gut microbial communities. In pigs, previous studies are based on classical "single-gene-single-trait" approaches and have evaluated the role of host genome controlling gut prokaryote and eukaryote communities separately.

RESULTS

In order to determine the ability of the host genome to control the diversity and composition of microbial communities in healthy pigs, we undertook genome-wide association studies (GWAS) for 39 microbial phenotypes that included 2 diversity indexes, and the relative abundance of 31 bacterial and six commensal protist genera in 390 pigs genotyped for 70 K SNPs. The GWAS results were processed through a 3-step analytical pipeline comprised of (1) association weight matrix; (2) regulatory impact factor; and (3) partial correlation and information theory. The inferred gene regulatory network comprised 3561 genes (within a 5 kb distance from a relevant SNP-P < 0.05) and 738,913 connections (SNP-to-SNP co-associations). Our findings highlight the complexity and polygenic nature of the pig gut microbial ecosystem. Prominent within the network were 5 regulators, PRDM15, STAT1, ssc-mir-371, SOX9 and RUNX2 which gathered 942, 607, 588, 284 and 273 connections, respectively. PRDM15 modulates the transcription of upstream regulators of WNT and MAPK-ERK signaling to safeguard naive pluripotency and regulates the production of Th1- and Th2-type immune response. The signal transducer STAT1 has long been associated with immune processes and was recently identified as a potential regulator of vaccine response to porcine reproductive and respiratory syndrome. The list of regulators was enriched for immune-related pathways, and the list of predicted targets includes candidate genes previously reported as associated with microbiota profile in pigs, mice and human, such as SLIT3, SLC39A8, NOS1, IL1R2, DAB1, TOX3, SPP1, THSD7B, ELF2, PIANP, A2ML1, and IFNAR1. Moreover, we show the existence of host-genetic variants jointly associated with the relative abundance of butyrate producer bacteria and host performance.

CONCLUSIONS

Taken together, our results identified regulators, candidate genes, and mechanisms linked with microbiome modulation by the host. They further highlight the value of the proposed analytical pipeline to exploit pleiotropy and the crosstalk between bacteria and protists as significant contributors to host-microbiome interactions and identify genetic markers and candidate genes that can be incorporated in breeding program to improve host-performance and microbial traits. Video Abstract.

Investigating the Effect of Imputed Structural Variants from Whole-Genome Sequence on Genome-Wide Association and Genomic Prediction in Dairy Cattle

Simple Summary

Structural variants are large changes to the DNA sequences that differ from individual to individual. We discovered and quality-controlled a set of 24,908 structural variants and used a technique called imputation to infer them into 35,588 Holstein and Jersey cattle. We then investigated whether the structural variants affected key dairy cattle traits such as milk production, fertility and overall conformation. Structural variants explained generally less than 10 percent of the phenotypic variation in these traits. Four of the structural variants were significantly associated with dairy cattle production traits. However, the inclusion of the structural variants in the genomic prediction model did not increase genomic prediction accuracy.

Abstract

Structural variations (SVs) are large DNA segments of deletions, duplications, copy number variations, inversions and translocations in a re-sequenced genome compared to a reference genome. They have been found to be associated with several complex traits in dairy cattle and could potentially help to improve genomic prediction accuracy of dairy traits. Imputation of SVs was performed in individuals genotyped with single-nucleotide polymorphism (SNP) panels without the expense of sequencing them. In this study, we generated 24,908 high-quality SVs in a total of 478 whole-genome sequenced Holstein and Jersey cattle. We imputed 4489 SVs with R2 > 0.5 into 35,568 Holstein and Jersey dairy cattle with 578,999 SNPs with two pipelines, FImpute and Eagle2.3-Minimac3. Genome-wide association studies for production, fertility and overall type with these 4489 SVs revealed four significant SVs, of which two were highly linked to significant SNP. We also estimated the variance components for SNP and SV models for these traits using genomic best linear unbiased prediction (GBLUP). Furthermore, we assessed the effect on genomic prediction accuracy of adding SVs to GBLUP models. The estimated percentage of genetic variance captured by SVs for production traits was up to 4.57% for milk yield in bulls and 3.53% for protein yield in cows. Finally, no consistent increase in genomic prediction accuracy was observed when including SVs in GBLUP.

Revisiting the genome-wide significance threshold for common variant GWAS

Over the last decade, GWAS meta-analyses have used a strict P-value threshold of 5 × 10-8 to classify associations as significant. Here, we use our current understanding of frequently studied traits including lipid levels, height, and BMI to revisit this genome-wide significance threshold. We compare the performance of studies using the P = 5 × 10-8 threshold in terms of true and false positive rate to other multiple testing strategies: (1) less stringent P-value thresholds, (2) controlling the FDR with the Benjamini-Hochberg and Benjamini-Yekutieli procedure, and (3) controlling the Bayesian FDR with posterior probabilities. We applied these procedures to re-analyze results from the Global Lipids and GIANT GWAS meta-analysis consortia and supported them with extensive simulation that mimics the empirical data. We observe in simulated studies with sample sizes ∼20,000 and >120,000 that relaxing the P-value threshold to 5 × 10-7 increased discovery at the cost of 18% and 8% of additional loci being false positive results, respectively. FDR and Bayesian FDR are well controlled for both sample sizes with a few exceptions that disappear under a less stringent definition of true positives and the two approaches yield similar results. Our work quantifies the value of using a relaxed P-value threshold in large studies to increase their true positive discovery but also show the excess false positive rates due to such actions in modest-sized studies. These results may guide investigators considering different thresholds in replication studies and downstream work such as gene-set enrichment or pathway analysis. Finally, we demonstrate the viability of FDR-controlling procedures in GWAS.

Genome-wide fine-mapping identifies pleiotropic and functional variants that predict many traits across global cattle populations

The difficulty in finding causative mutations has hampered their use in genomic prediction. Here, we present a methodology to fine-map potentially causal variants genome-wide by integrating the functional, evolutionary and pleiotropic information of variants using GWAS, variant clustering and Bayesian mixture models. Our analysis of 17 million sequence variants in 44,000+ Australian dairy cattle for 34 traits suggests, on average, one pleiotropic QTL existing in each 50 kb chromosome-segment. We selected a set of 80k variants representing potentially causal variants within each chromosome segment to develop a bovine XT-50K genotyping array. The custom array contains many pleiotropic variants with biological functions, including splicing QTLs and variants at conserved sites across 100 vertebrate species. This biology-informed custom array outperformed the standard array in predicting genetic value of multiple traits across populations in independent datasets of 90,000+ dairy cattle from the USA, Australia and New Zealand.

Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits

Background

Variants that regulate transcription, such as expression quantitative trait loci (eQTL), have shown enrichment in genome-wide association studies (GWAS) for mammalian complex traits. However, no study has reported eQTL in sheep, although it is an important agricultural species for which many GWAS of complex meat traits have been conducted. Using RNA sequence data produced from liver and muscle from 149 sheep and imputed whole-genome single nucleotide polymorphisms (SNPs), our aim was to dissect the genetic architecture of the transcriptome by associating sheep genotypes with three major molecular phenotypes including gene expression (geQTL), exon expression (eeQTL) and RNA splicing (sQTL). We also examined these three types of eQTL for their enrichment in GWAS of multi-meat traits and fatty acid profiles.

Results

Whereas a relatively small number of molecular phenotypes were significantly heritable (h² > 0, P < 0.05), their mean heritability ranged from 0.67 to 0.73 for liver and from 0.71 to 0.77 for muscle. Association analysis between molecular phenotypes and SNPs within ± 1 Mb identified many significant cis-eQTL (false discovery rate, FDR < 0.01). The median distance between the eQTL and transcription start sites (TSS) ranged from 68 to 153 kb across the three eQTL types. The number of common variants between geQTL, eeQTL and sQTL within each tissue, and the number of common variants between liver and muscle within each eQTL type were all significantly (P < 0.05) larger than expected by chance. The identified eQTL were significantly (P < 0.05) enriched in GWAS hits associated with 56 carcass traits and fatty acid profiles. For example, several geQTL in muscle mapped to the FAM184B gene, hundreds of sQTL in liver and muscle mapped to the CAST gene, and hundreds of sQTL in liver mapped to the C6 gene. These three genes are associated with body composition or fatty acid profiles.

Conclusions

We detected a large number of significant eQTL and found that the overlap of variants between eQTL types and tissues was prevalent. Many eQTL were also QTL for meat traits. Our study fills a gap in the knowledge on the regulatory variants and their role in complex traits for the sheep model.

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.

GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci

BACKGROUND

The development of high-throughput sequencing and genotyping methodologies has enabled the identification of thousands of genomic regions associated with several complex traits. The integration of multiple sources of biological information is a crucial step required to better understand patterns regulating the development of these traits.

FINDINGS

Genomic Annotation in Livestock for positional candidate LOci (GALLO) is an R package developed for the accurate annotation of genes and quantitative trait loci (QTLs) located in regions identified in common genomic analyses performed in livestock, such as genome-wide association studies and transcriptomics using RNA sequencing. Moreover, GALLO allows the graphical visualization of gene and QTL annotation results, data comparison among different grouping factors (e.g., methods, breeds, tissues, statistical models, studies), and QTL enrichment in different livestock species such as cattle, pigs, sheep, and chickens.

CONCLUSIONS

Consequently, GALLO is a useful package for annotation, identification of hidden patterns across datasets, and data mining previously reported associations, as well as the efficient examination of the genetic architecture of complex traits in livestock.

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.

A multi-breed GWAS for morphometric traits in four Beninese indigenous cattle breeds reveals loci associated with conformation, carcass and adaptive traits

BACKGROUND

Specific adaptive features including disease resistance and growth abilities in harsh environments are attributed to indigenous cattle breeds of Benin, but these breeds are endangered due to crossbreeding. So far, there is a lack of systematic trait recording, being the basis for breed characterizations, and for structured breeding program designs aiming on conservation. Bridging this gap, own phenotyping for morphological traits considered measurements for height at withers (HAW), sacrum height (SH), heart girth (HG), hip width (HW), body length (BL) and ear length (EL), including 449 cattle from the four indigenous Benin breeds Lagune, Somba, Borgou and Pabli. In order to utilize recent genomic tools for breed characterizations and genetic evaluations, phenotypes for novel traits were merged with high-density SNP marker data. Multi-breed genetic parameter estimations and genome-wide association studies (GWAS) for the six morphometric traits were carried out. Continuatively, we aimed on inferring genomic regions and functional loci potentially associated with conformation, carcass and adaptive traits.

RESULTS

SNP-based heritability estimates for the morphometric traits ranged between 0.46 ± 0.14 (HG) and 0.74 ± 0.13 (HW). Phenotypic and genetic correlations ranged from 0.25 ± 0.05 (HW-BL) to 0.89 ± 0.01 (HAW-SH), and from 0.14 ± 0.10 (HW-BL) to 0.85 ± 0.02 (HAW-SH), respectively. Three genome-wide and 25 chromosome-wide significant SNP positioned on different chromosomes were detected, located in very close chromosomal distance (±25 kb) to 15 genes (or located within the genes). The genes PIK3R6 and PIK3R1 showed direct functional associations with height and body size. We inferred the potential candidate genes VEPH1, CNTNAP5, GYPC for conformation, growth and carcass traits including body weight and body fat deposition. According to their functional annotations, detected potential candidate genes were associated with stress or immune response (genes PTAFR, PBRM1, ADAMTS12) and with feed efficiency (genes MEGF11 SLC16A4, CCDC117).

CONCLUSIONS

Accurate measurements contributed to large SNP heritabilities for some morphological traits, even for a small mixed-breed sample size. Multi-breed GWAS detected different loci associated with conformation or carcass traits. The identified potential candidate genes for immune response or feed efficiency indicators reflect the evolutionary development and adaptability features of the breeds.

Exploring natural diversity of Medicago truncatula reveals physiotypes and loci associated with the response of seedling performance to nitrate supply

Seedling pre-emergence is a critical phase of development for successful crop establishment because of its susceptibility to environmental conditions. In a context of reduced use of inorganic fertilizers, the genetic bases of the response of seedlings to nitrate supply received little attention. This issue is important even in legumes where nitrate absorption starts early after germination, before nodule development. Natural variation of traits characterizing seedling growth in the absence or presence of nitrate was investigated in a core collection of 192 accessions of Medicago truncatula. Plasticity indexes to the absence of nitrate were calculated. The genetic determinism of the traits was dissected by genome-wide association study (GWAS). The absence of nitrate affected seed biomass mobilization and root/shoot length ratio. However, the large range of genetic variability revealed different seedling performances within natural diversity. A principal component analysis (PCA) carried out with plasticity indexes highlighted four physiotypes of accessions differing in relationships between seedling elongation and seed biomass partitioning traits in response to the absence of nitrate. Finally, GWAS revealed 45 associations with single or combined traits corresponding to coordinates of accessions on PCA, as well as two clusters of genes encoding sugar transporters and glutathione transferases surrounding loci associated with seedling elongation traits.

Quantitative and Qualitative Role of Antagonistic Heterogeneity in Genetics of Blood Lipids

Prevailing strategies in genome-wide association studies (GWAS) mostly rely on principles of medical genetics emphasizing one gene, one function, one phenotype concept. Here, we performed GWAS of blood lipids leveraging a new systemic concept emphasizing complexity of genetic predisposition to such phenotypes. We focused on total cholesterol, low- and high-density lipoprotein cholesterols, and triglycerides available for 29,902 individuals of European ancestry from seven independent studies, men and women combined. To implement the new concept, we leveraged the inherent heterogeneity in genetic predisposition to such complex phenotypes and emphasized a new counter intuitive phenomenon of antagonistic genetic heterogeneity, which is characterized by misalignment of the directions of genetic effects and the phenotype correlation. This analysis identified 37 loci associated with blood lipids but only one locus, FBXO33, was not reported in previous top GWAS. We, however, found strong effect of antagonistic heterogeneity that leaded to profound (quantitative and qualitative) changes in the associations with blood lipids in most, 25 of 37 or 68%, loci. These changes suggested new roles for some genes, which functions were considered as well established such as GCKR, SIK3 (APOA1 locus), LIPC, LIPG, among the others. The antagonistic heterogeneity highlighted a new class of genetic associations emphasizing beneficial and adverse trade-offs in predisposition to lipids. Our results argue that rigorous analyses dissecting heterogeneity in genetic predisposition to complex traits such as lipids beyond those implemented in current GWAS are required to facilitate translation of genetic discoveries into health care.

Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

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.

Gene Networks Driving Genetic Variation in Milk and Cheese-Making Traits of Spanish Assaf Sheep

Most of the milk produced by sheep is used for the production of high-quality cheese. Consequently, traits related to milk coagulation properties and cheese yield are economically important to the Spanish dairy industry. The present study aims to identify candidate genes and their regulators related to 14 milk and cheese-making traits and to develop a low-density panel of markers that could be used to predict an individual's genetic potential for cheese-making efficiency. In this study, we performed a combination of the classical genome-wide association study (GWAS) with a stepwise regression method and a pleiotropy analysis to determine the best combination of the variants located within the confidence intervals of the potential candidate genes that may explain the greatest genetic variance for milk and cheese-making traits. Two gene networks related to milk and cheese-making traits were created using the genomic relationship matrices built through a stepwise multiple regression approach. Several co-associated genes in these networks are involved in biological processes previously found to be associated with milk synthesis and cheese-making efficiency. The methodology applied in this study enabled the selection of a co-association network comprised of 374 variants located in the surrounding of genes showing a potential influence on milk synthesis and cheese-making efficiency.

Integrating Genome-Wide CNVs Into QTLs and High Confidence GWAScore Regions Identified Positional Candidates for Sheep Economic Traits

Copy number variations (CNVs) are important source of genetic variation, which can affect diverse economic traits through a variety of mechanisms. In addition, genome scan can identify many quantitative trait loci (QTLs) for the economic traits, while genome-wide association studies (GWAS) can localize genetic variants associated with the phenotypic variations. Here, we developed a method called GWAScore which collected GWAS summary data to identify potential candidates, and integrated CNVs into QTLs and high confidence GWAScore regions to detect crucial CNV markers for sheep growth traits. We got 197 candidate genes which were overlapping with the candidate CNVs. Some crucial genes (MYLK3, TTC29, HERC6, ABCG2, RUNX1, etc.) showed significantly elevated GWAScore peaks than other candidate genes. In this study, we developed the GWAScore method to excavate the potential value of candidate genes as markers for the sheep molecular breeding.

Distinguishing pleiotropy from linked QTL between milk production traits and mastitis resistance in Nordic Holstein cattle

BACKGROUND

Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that affect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Effect-bearing variants often affect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic effects of variants in candidate genes. However, large sample sizes are required to achieve sufficient power. Multi-trait meta-analysis is one approach to deal with this problem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance.

RESULTS

For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis compared with the subjective assessment of overlapping of single-trait QTL confidence intervals. Pleiotropic effects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confirmed by bivariate association analysis. The previously reported pleiotropic effects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic effects of variants in GHR on milk yield and fat yield were confirmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic effects on milk production traits. For milk yield and mastitis resistance, we identified possible pleiotropic effects of variants in two genes, GC and DGAT1.

CONCLUSIONS

Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk production traits and identifies variants with pleiotropic effects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correlation between milk yield and mastitis.

Genome-wide association studies for yield component traits in a macadamia breeding population

Background

Breeding for new macadamia cultivars with high nut yield is expensive in terms of time, labour and cost. Most trees set nuts after four to five years, and candidate varieties for breeding are evaluated for at least eight years for various traits. Genome-wide association studies (GWAS) are promising methods to reduce evaluation and selection cycles by identifying genetic markers linked with key traits, potentially enabling early selection through marker-assisted selection. This study used 295 progeny from 32 full-sib families and 29 parents (18 phenotyped) which were planted across four sites, with each tree genotyped for 4113 SNPs. ASReml-R was used to perform association analyses with linear mixed models including a genomic relationship matrix to account for population structure. Traits investigated were: nut weight (NW), kernel weight (KW), kernel recovery (KR), percentage of whole kernels (WK), tree trunk circumference (TC), percentage of racemes that survived from flowering through to nut set, and number of nuts per raceme.

Results

Seven SNPs were significantly associated with NW (at a genome-wide false discovery rate of < 0.05), and four with WK. Multiple regression, as well as mapping of markers to genome assembly scaffolds suggested that some SNPs were detecting the same QTL. There were 44 significant SNPs identified for TC although multiple regression suggested detection of 16 separate QTLs.

Conclusions

These findings have important implications for macadamia breeding, and highlight the difficulties of heterozygous populations with rapid LD decay. By coupling validated marker-trait associations detected through GWAS with MAS, genetic gain could be increased by reducing the selection time for economically important nut characteristics. Genomic selection may be a more appropriate method to predict complex traits like tree size and yield.

Effect direction meta-analysis of GWAS identifies extreme, prevalent and shared pleiotropy in a large mammal

In genome-wide association studies (GWAS), variants showing consistent effect directions across populations are considered as true discoveries. We model this information in an Effect Direction MEta-analysis (EDME) to quantify pleiotropy using GWAS of 34 Cholesky-decorrelated traits in 44,000+ cattle with sequence variants. The effect-direction agreement between independent bull and cow datasets was used to quantify the false discovery rate by effect direction (FDRed) and the number of affected traits for prioritised variants. Variants with multi-trait p < 1e–6 affected 1∼22 traits with an average of 10 traits. EDME assigns pleiotropic variants to each trait which informs the biology behind complex traits. New pleiotropic loci are identified, including signals from the cattle FTO locus mirroring its bystander effects on human obesity. When validated in the 1000-Bull Genome database, the prioritized pleiotropic variants consistently predicted expected phenotypic differences between dairy and beef cattle. EDME provides robust approaches to control GWAS FDR and quantify pleiotropy.

Xiang et al. developed an Effect Direction Meta-analysis (EDME) approach to identify true pleiotropy. They used Cholesky-transformation to decorrelate the traits and identified many pleiotropic variants that consistently predicted phenotypic differences in cattle.

Genomic Regions Associated With Skeletal Type Traits in Beef and Dairy Cattle Are Common to Regions Associated With Carcass Traits, Feed Intake and Calving Difficulty

Linear type traits describing the skeletal characteristics of an animal are moderately to strongly genetically correlated with a range of other performance traits in cattle including feed intake, reproduction traits and carcass merit; thus, type traits could also provide useful insights into the morphological differences among animals underpinning phenotypic differences in these complex traits. The objective of the present study was to identify genomic regions associated with five subjectively scored skeletal linear traits, to determine if these associated regions are common in multiple beef and dairy breeds, and also to determine if these regions overlap with those proposed elsewhere to be associated with correlated performance traits. Analyses were carried out using linear mixed models on imputed whole genome sequence data separately in 1,444 Angus, 1,129 Hereford, 6,433 Charolais, 8,745 Limousin, 1,698 Simmental, and 4,494 Holstein-Friesian cattle, all scored for the linear type traits. There was, on average, 18 months difference in age at assessment of the beef versus the dairy animals. While the majority of the identified quantitative trait loci (QTL), and thus genes, were both trait-specific and breed-specific, a large-effect pleiotropic QTL on BTA6 containing the NCAPG and LCORL genes was associated with all skeletal traits in the Limousin population and with wither height in the Angus. Other than that, little overlap existed in detected QTLs for the skeletal type traits in the other breeds. Only two QTLs overlapped the beef and dairy breeds; both QTLs were located on BTA5 and were associated with height in both the Angus and the Holstein-Friesian, despite the difference in age at assessment. Several detected QTLs in the present study overlapped with QTLs documented elsewhere that are associated with carcass traits, feed intake, and calving difficulty. While most breeding programs select for the macro-traits like carcass weight, carcass conformation, and feed intake, the higher degree of granularity with selection on the individual linear type traits in a multi-trait index underpinning the macro-level goal traits, presents an opportunity to help resolve genetic antagonisms among morphological traits in the pursuit of the animal with optimum performance metrics.

Genetic correlations of polygenic disease traits: from theory to practice

The genetic correlation describes the genetic relationship between two traits and can contribute to a better understanding of the shared biological pathways and/or the causality relationships between them. The rarity of large family cohorts with recorded instances of two traits, particularly disease traits, has made it difficult to estimate genetic correlations using traditional epidemiological approaches. However, advances in genomic methodologies, such as genome-wide association studies, and widespread sharing of data now allow genetic correlations to be estimated for virtually any trait pair. Here, we review the definition, estimation, interpretation and uses of genetic correlations, with a focus on applications to human disease.

Genome-wide association and genotype by environment interactions for growth traits in U.S. Gelbvieh cattle

BACKGROUND

Single nucleotide polymorphism (SNP) arrays have facilitated discovery of genetic markers associated with complex traits in domestic cattle; thereby enabling modern breeding and selection programs. Genome-wide association analyses (GWAA) for growth traits were conducted on 10,837 geographically diverse U.S. Gelbvieh cattle using a union set of 856,527 imputed SNPs. Birth weight (BW), weaning weight (WW), and yearling weight (YW) were analyzed using GEMMA and EMMAX (via imputed genotypes). Genotype-by-environment (GxE) interactions were also investigated.

RESULTS

GEMMA and EMMAX produced moderate marker-based heritability estimates that were similar for BW (0.36-0.37, SE = 0.02-0.06), WW (0.27-0.29, SE = 0.01), and YW (0.39-0.41, SE = 0.01-0.02). GWAA using 856K imputed SNPs (GEMMA; EMMAX) revealed common positional candidate genes underlying pleiotropic QTL for Gelbvieh growth traits on BTA6, BTA7, BTA14, and BTA20. The estimated proportion of phenotypic variance explained (PVE) by the lead SNP defining these QTL (EMMAX) was larger and most similar for BW and YW, and smaller for WW. Collectively, GWAAs (GEMMA; EMMAX) produced a highly concordant set of BW, WW, and YW QTL that met a nominal significance level (P ≤ 1e-05), with prioritization of common positional candidate genes; including genes previously associated with stature, feed efficiency, and growth traits (i.e., PLAG1, NCAPG, LCORL, ARRDC3, STC2). Genotype-by-environment QTL were not consistent among traits at the nominal significance threshold (P ≤ 1e-05); although some shared QTL were apparent at less stringent significance thresholds (i.e., P ≤ 2e-05).

CONCLUSIONS

Pleiotropic QTL for growth traits were detected on BTA6, BTA7, BTA14, and BTA20 for U.S. Gelbvieh beef cattle. Seven QTL detected for Gelbvieh growth traits were also recently detected for feed efficiency and growth traits in U.S. Angus, SimAngus, and Hereford cattle. Marker-based heritability estimates and the detection of pleiotropic QTL segregating in multiple breeds support the implementation of multiple-breed genomic selection.

A genome-wide association study of reproduction traits in four pig populations with different genetic backgrounds

Objective

Genome-wide association study and two meta-analysis based on GWAS performed to explore the genetic mechanism underlying variation in pig number born alive (NBA) and total number born (TNB).

Methods

Single trait GWAS and two meta-analysis (single-trait meta analysis and multi-trait meta analysis) were used in our study for NBA and TNB on 3,121 Yorkshires from 4 populations, including three different American Yorkshire populations (n = 2,247) and one British Yorkshire populations (n = 874).

Results

The result of single trait GWAS showed that no significant associated single nucleotide polymorphisms (SNPs) were identified. Using single-trait meta analysis and multi-trait meta analysis within populations, 11 significant loci were identified associated with target traits. Spindlin 1, vascular endothelial growth factor A, forkhead box Q1, msh homeobox 1, and LHFPL tetraspan submily member 3 are five functionally plausible candidate genes for NBA and TNB. Compared to the single population GWAS, single-trait Meta analysis can improve the detection power to identify SNPs by integrating information of multiple populations. The multiple-trait analysis reduced the power to detect trait-specific loci but enhanced the power to identify the common loci across traits.

Conclusion

In total, our findings identified novel genes to be validated as candidates for NBA and TNB in pigs. Also, it enabled us to enlarge population size by including multiple populations with different genetic backgrounds and increase the power of GWAS by using meta analysis.

Identification of muscle-specific candidate genes in Simmental beef cattle using imputed next generation sequencing

Genome-wide association studies (GWAS) have commonly been used to identify candidate genes that control economically important traits in livestock. Our objective was to detect potential candidate genes associated mainly with muscle development traits related to dimension of hindquarter in cattle. A next generation sequencing (NGS) dataset to imputed to 12 million single nucleotide polymorphisms (SNPs) (from 1252 Simmental beef cattle) were used to search for genes affecting hindquarter traits using a linear, mixed model approach. We also used haplotype and linkage disequilibrium blocks to further support our identifications. We identified 202 significant SNPs in the bovine BTA4 chromosome region associated with width of hind leg, based on a stringent statistical threshold (p = 0.05/ effective number of SNPs identified). After exploring the region around these SNPs, we found candidate genes that were potentially related to the associated markers. More importantly, we identified a region of approximately 280 Kb on the BTA4 chromosome that harbored several muscle specific candidate genes, genes to be in a potential region for muscle development. However, we also found candidate gene SLC13A1 on BTA4, which seems to be associated with bone disorders (such as chondrodysplasia) in Simmental beef cattle.

Quantifying the contribution of sequence variants with regulatory and evolutionary significance to 34 bovine complex traits

Many genome variants shaping mammalian phenotype are hypothesized to regulate gene transcription and/or to be under selection. However, most of the evidence to support this hypothesis comes from human studies. Systematic evidence for regulatory and evolutionary signals contributing to complex traits in a different mammalian model is needed. Sequence variants associated with gene expression (expression quantitative trait loci [eQTLs]) and concentration of metabolites (metabolic quantitative trait loci [mQTLs]) and under histone-modification marks in several tissues were discovered from multiomics data of over 400 cattle. Variants under selection and evolutionary constraint were identified using genome databases of multiple species. These analyses defined 30 sets of variants, and for each set, we estimated the genetic variance the set explained across 34 complex traits in 11,923 bulls and 32,347 cows with 17,669,372 imputed variants. The per-variant trait heritability of these sets across traits was highly consistent (r > 0.94) between bulls and cows. Based on the per-variant heritability, conserved sites across 100 vertebrate species and mQTLs ranked the highest, followed by eQTLs, young variants, those under histone-modification marks, and selection signatures. From these results, we defined a Functional-And-Evolutionary Trait Heritability (FAETH) score indicating the functionality and predicted heritability of each variant. In additional 7,551 cattle, the high FAETH-ranking variants had significantly increased genetic variances and genomic prediction accuracies in 3 production traits compared to the low FAETH-ranking variants. The FAETH framework combines the information of gene regulation, evolution, and trait heritability to rank variants, and the publicly available FAETH data provide a set of biological priors for cattle genomic selection worldwide.

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.

RAPID COMMUNICATION: Multi-breed validation study unraveled genomic regions associated with puberty traits segregating across tropically adapted breeds1

An efficient strategy to improve QTL detection power is performing across-breed validation studies. Variants segregating across breeds are expected to be in high linkage disequilibrium (LD) with causal mutations affecting economically important traits. The aim of this study was to validate, in a Tropical Composite cattle (TC) population, QTL associations identified for sexual precocity traits in a Nellore and Brahman meta-analysis genome-wide association study. In total, 2,816 TC, 8,001 Nellore, and 2,210 Brahman animals were available for the analysis. For that, genomic regions significantly associated with puberty traits in the meta-analysis study were validated for the following sexual precocity traits in TC: age at first corpus luteum (AGECL), first postpartum anestrus interval (PPAI), and scrotal circumference at 18 months of age (SC). We considered validated QTL those underpinned by significant markers from the Nellore and Brahman meta-analysis (P ≤ 10-4) that were also significant for a TC trait, i.e., presenting a P-value of ≤10-3 for AGECL, PPAI, or SC. We also considered as validated QTL those regions where significant markers in the reference population were at ±250 kb from significant markers in the validation population. Using this criteria, 49 SNP were validated for AGECL, 4 for PPAI, and 14 for SC, from which 5 were in common with AGECL, totaling 62 validated SNP for these traits and 30 candidate genes surrounding them. Considering just candidate genes closest to the top SNP of each chromosome, for AGECL 8 candidate genes were identified: COL8A1, PENK, ENSBTAG00000047425, BPNT1, ADAMTS17, CCHCR1, SUFU, and ENSBTAG00000046374. For PPAI, 3 genes emerged as candidates (PCBP3, KCNK10, and MRPS5), and for SC 8 candidate genes were identified (SNORA70, TRAC, ASS1, BPNT1, LRRK1, PKHD1, PTPRM, and ENSBTAG00000045690). Several candidate regions presented here were previously associated with puberty traits in cattle. The majority of emerging candidate genes are related to biological processes involved in reproductive events, such as maintenance of gestation, and some are known to be expressed in reproductive tissues. Our results suggested that some QTL controlling early puberty seem to be segregating across cattle breeds adapted to tropical conditions.

Multi-trait meta-analyses reveal 25 quantitative trait loci for economically important traits in Brown Swiss cattle

Background

Little is known about the genetic architecture of economically important traits in Brown Swiss cattle because only few genome-wide association studies (GWAS) have been carried out in this breed. Moreover, most GWAS have been performed for single traits, thus not providing detailed insights into potentially existing pleiotropic effects of trait-associated loci.

Results

To compile a comprehensive catalogue of large-effect quantitative trait loci (QTL) segregating in Brown Swiss cattle, we carried out association tests between partially imputed genotypes at 598,016 SNPs and daughter-derived phenotypes for more than 50 economically important traits, including milk production, growth and carcass quality, body conformation, reproduction and calving traits in 4578 artificial insemination bulls from two cohorts of Brown Swiss cattle (Austrian-German and Swiss populations). Across-cohort multi-trait meta-analyses of the results from the single-trait GWAS revealed 25 quantitative trait loci (QTL; P < 8.36 × 10^− 8) for economically relevant traits on 17 Bos taurus autosomes (BTA). Evidence of pleiotropy was detected at five QTL located on BTA5, 6, 17, 21 and 25. Of these, two QTL at BTA6:90,486,780 and BTA25:1,455,150 affect a diverse range of economically important traits, including traits related to body conformation, calving, longevity and milking speed. Furthermore, the QTL at BTA6:90,486,780 seems to be a target of ongoing selection as evidenced by an integrated haplotype score of 2.49 and significant changes in allele frequency over the past 25 years, whereas either no or only weak evidence of selection was detected at all other QTL.

Conclusions

Our findings provide a comprehensive overview of QTL segregating in Brown Swiss cattle. Detected QTL explain between 2 and 10% of the variation in the estimated breeding values and thus may be considered as the most important QTL segregating in the Brown Swiss cattle breed. Multi-trait association testing boosts the power to detect pleiotropic QTL and assesses the full spectrum of phenotypes that are affected by trait-associated variants.

Electronic supplementary material

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

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 and evolutionary analyses reveal the formation of swine facial wrinkles in Chinese Erhualian pigs

Wrinkles are uneven concave-convex folds, ridges or creases in skin. Facial wrinkles appear in head, typically increasing along with aging. However in several Chinese indigenous pigs, such as Erhualian pigs, rich facial wrinkles have been generated during the growth stages as one of their breed characteristics. To investigate the genetic basis underlying the development of swine facial wrinkles, we estimated the folding extent of facial wrinkles in a herd of Erhualian pigs (n=332), and then conducted genome-wide association studies and multi-trait meta-analysis for facial wrinkles using 60K porcine chips. We found that facial wrinkles had high heritability estimates of ~0.7 in Erhualian pigs. Notably, only one genome-wide significant QTL was detected at 34.8 Mb on porcine chromosome 7. The most significant SNP rs80983858 located at the 3255-bp downstream of candidate gene GRM4, and the G allele was of benefit to increase facial wrinkles. Evolutionary and selection analyses suggested that the haplotypes containing G allele were under artificial selection, which was consistent with local animal sacrificial custom praying for longevity. Our findings made important clues for further deciphering the molecular mechanism of swine facial wrinkles formation, and shed light on the research of skin wrinkle development in human or other mammals.

Using imputed whole-genome sequence data to improve the accuracy of genomic prediction for parasite resistance in Australian sheep

Background

This study aimed at (1) comparing the accuracies of genomic prediction for parasite resistance in sheep based on whole-genome sequence (WGS) data to those based on 50k and high-density (HD) single nucleotide polymorphism (SNP) panels; (2) investigating whether the use of variants within quantitative trait loci (QTL) regions that were selected from regional heritability mapping (RHM) in an independent dataset improved the accuracy more than variants selected from genome-wide association studies (GWAS); and (3) comparing the prediction accuracies between variants selected from WGS data to variants selected from the HD SNP panel.

Results

The accuracy of genomic prediction improved marginally from 0.16 ± 0.02 and 0.18 ± 0.01 when using all the variants from 50k and HD genotypes, respectively, to 0.19 ± 0.01 when using all the variants from WGS data. Fitting a GRM from the selected variants alongside a GRM from the 50k SNP genotypes improved the prediction accuracy substantially compared to fitting the 50k SNP genotypes alone. The gain in prediction accuracy was slightly more pronounced when variants were selected from WGS data compared to when variants were selected from the HD panel. When sequence variants that passed the GWAS \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- log_{10} (p\,value)$$\end{document}-log10(pvalue) threshold of 3 across the entire genome were selected, the prediction accuracy improved by 5% (up to 0.21 ± 0.01), whereas when selection was limited to sequence variants that passed the same GWAS \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- log_{10} (p\,value)$$\end{document}-log10(pvalue) threshold of 3 in regions identified by RHM, the accuracy improved by 9% (up to 0.25 ± 0.01).

Conclusions

Our results show that through careful selection of sequence variants from the QTL regions, the accuracy of genomic prediction for parasite resistance in sheep can be improved. These findings have important implications for genomic prediction in sheep.

Functional annotation and Bayesian fine-mapping reveals candidate genes for important agronomic traits in Holstein bulls

A hundred years of data collection in dairy cattle can facilitate powerful studies of complex traits. Cattle GWAS have identified many associated genomic regions. With increasing numbers of cattle sequenced, fine-mapping of causal variants is becoming possible. Here we imputed selected sequence variants to 27,214 Holstein bulls that have highly reliable phenotypes for 35 production, reproduction, and body conformation traits. We performed single-marker scans for the 35 traits and multi-trait tests of the three trait groups, revealing 282 candidate QTL for fine-mapping. We developed a Bayesian Fine-MAPping approach (BFMAP) to integrate fine-mapping with functional enrichment analysis. Our fine-mapping identified 69 promising candidate genes, including ABCC9, VPS13B, MGST1, SCD, MKL1, CSN1S1 for production, CHEK2, GC, KALRN for reproduction, and TMTC2, ARRDC3, ZNF613, CCND2, FGF6 for conformation traits. Collectively, these results demonstrated the utility of BFMAP, identified candidate genes, and enhanced our understanding of the genetic basis of cattle complex traits.

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.