School of Agriculture, Food and Ecosystem Sciences - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 102
-
ItemErratum to "Evaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle" (JDS Commun. 3:114-119).(American Dairy Science Association, 2022-09)[This corrects the article DOI: 10.3168/jdsc.2021-0150.].
-
ItemAdditive Genetic Variation in Schizophrenia Risk Is Shared by Populations of African and European Descent(CELL PRESS, 2013-09-05)To investigate the extent to which the proportion of schizophrenia's additive genetic variation tagged by SNPs is shared by populations of European and African descent, we analyzed the largest combined African descent (AD [n = 2,142]) and European descent (ED [n = 4,990]) schizophrenia case-control genome-wide association study (GWAS) data set available, the Molecular Genetics of Schizophrenia (MGS) data set. We show how a method that uses genomic similarities at measured SNPs to estimate the additive genetic correlation (SNP correlation [SNP-rg]) between traits can be extended to estimate SNP-rg for the same trait between ethnicities. We estimated SNP-rg for schizophrenia between the MGS ED and MGS AD samples to be 0.66 (SE = 0.23), which is significantly different from 0 (p(SNP-rg = 0) = 0.0003), but not 1 (p(SNP-rg = 1) = 0.26). We re-estimated SNP-rg between an independent ED data set (n = 6,665) and the MGS AD sample to be 0.61 (SE = 0.21, p(SNP-rg = 0) = 0.0003, p(SNP-rg = 1) = 0.16). These results suggest that many schizophrenia risk alleles are shared across ethnic groups and predate African-European divergence.
-
ItemGenetic variation in histone modifications and gene expression identifies regulatory variants in the mammary gland of cattle(BMC, 2022-12-08)BACKGROUND: Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. RESULTS: We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. CONCLUSIONS: We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels.
-
ItemSharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency(BMC, 2022-09-06)BACKGROUND: Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows. RESULTS: GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (rg) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase. CONCLUSIONS: The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.
-
ItemBayesR3 enables fast MCMC blocked processing for largescale multi-trait genomic prediction and QTN mapping analysis.(Springer Science and Business Media LLC, 2022-07-05)Bayesian methods, such as BayesR, for predicting the genetic value or risk of individuals from their genotypes, such as Single Nucleotide Polymorphisms (SNP), are often implemented using a Markov Chain Monte Carlo (MCMC) process. However, the generation of Markov chains is computationally slow. We introduce a form of blocked Gibbs sampling for estimating SNP effects from Markov chains that greatly reduces computational time by sampling each SNP effect iteratively n-times from conditional block posteriors. Subsequent iteration over all blocks m-times produces chains of length m × n. We use this strategy to solve large-scale genomic prediction and fine-mapping problems using the Bayesian MCMC mixed-effects genetic model, BayesR3. We validate the method using simulated data, followed by analysis of empirical dairy cattle data using high dimension milk mid infra-red spectra data as an example of "omics" data and show its use to increase the precision of mapping variants affecting milk, fat, and protein yields relative to a univariate analysis of milk, fat, and protein.
-
ItemMetaGS: an accurate method to impute and combine SNP effects across populations using summary statistics.(Springer Science and Business Media LLC, 2022-06-02)BACKGROUND: Meta-analysis describes a category of statistical methods that aim at combining the results of multiple studies to increase statistical power by exploiting summary statistics. Different industries that use genomic prediction do not share their raw data due to logistic or privacy restrictions, which can limit the size of their reference populations and creates a need for a practical meta-analysis method. RESULTS: We developed a meta-analysis, named MetaGS, that duplicates the results of multi-trait best linear unbiased prediction (mBLUP) analysis without accessing raw data. MetaGS exploits the correlations among different populations to produce more accurate population-specific single nucleotide polymorphism (SNP) effects. The method improves SNP effect estimations for a given population depending on its relations to other populations. MetaGS was tested on milk, fat and protein yield data of Australian Holstein and Jersey cattle and it generated very similar genomic estimated breeding values to those produced using the mBLUP method for all traits in both breeds. One of the major difficulties when combining SNP effects across populations is the use of different variants for the populations, which limits the applications of meta-analysis in practice. We solved this issue by developing a method to impute missing summary statistics without using raw data. Our results showed that imputing summary statistics can be done with high accuracy (r > 0.9) even when more than 70% of the SNPs were missing with a minimal effect on prediction accuracy. CONCLUSIONS: We demonstrated that MetaGS can replace the mBLUP model when raw data cannot be shared, which can lead to more flexible collaborations compared to the single-trait BLUP model.
-
ItemNo Preview AvailablePhantom epistasis between unlinked loci(NATURE PORTFOLIO, 2021-08-12)
-
ItemAuthor Correction: Assortative mating biases marker-based heritability estimators.(Springer Science and Business Media LLC, 2022-04-01)
-
ItemAssortative mating biases marker-based heritability estimators(NATURE PORTFOLIO, 2022-02-03)Many traits are subject to assortative mating, with recent molecular genetic findings confirming longstanding theoretical predictions that assortative mating induces long range dependence across causal variants. However, all marker-based heritability estimators implicitly assume mating is random. We provide mathematical and simulation-based evidence demonstrating that both method-of-moments and likelihood-based estimators are biased in the presence of assortative mating and derive corrected heritability estimators for traits subject to assortment. Finally, we demonstrate that the empirical patterns of estimates across methods and sample sizes for real traits subject to assortative mating are congruent with expected assortative mating-induced biases. For example, marker-based heritability estimates for height are 14% - 23% higher than corrected estimates using UK Biobank data.
-
ItemEvaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle.(American Dairy Science Association, 2022-03)Although selection for increased milk production traits has led to a genetic increase in body weight (BW), the genetic gain in milk production has exceeded the gain in BW, so gross feed efficiency has improved. Nonetheless, greater gains may be possible by directly selecting for a measure of feed efficiency. Australia first introduced Feed Saved (FS) estimated breeding value (EBV) in 2015. Feed Saved combines residual feed intake (RFI) genomic EBV and maintenance requirements calculated from mature BW EBV. The FS EBV was designed to enable the selection of cows for reduced energy requirements with similar milk production. In this study, we used a reference population of 3,711 animals in a multivariate analysis including Australian heifers (AUSh), Australian cows (AUSc), and overseas cows (OVEc) to update the Australian EBV for lifetime RFI (i.e., a breeding value that incorporated RFI in growing and lactating cows) and to recalculate the FS EBV in Australian Holstein bulls (AUSb). The estimates of genomic heritabilities using univariate (only AUSc or AUSh) to trivariate (including the OVEc) analyses were similar. Genomic heritabilities for RFI were estimated as 0.18 for AUSc, 0.27 for OVEc, and 0.36 for AUSh. The genomic correlation for RFI between AUSc and AUSh was 0.47 and that between AUSc and OVEc was 0.94, but these estimates were associated with large standard errors (range: 0.18-0.28). The reliability of lifetime RFI (a component of FS) in the trivariate analysis (i.e., including OVEc) increased from 11% to 20% compared with the 2015 model and was greater, by 12%, than in a bivariate analysis in which the reference population included only AUSc and AUSh. By applying the prediction equation of the 2020 model, the average reliability of the FS EBV in 20,816 AUSb that were born between 2010 and 2020 improved from 33% to 43%. Previous selection strategies-that is, using the predecessor of the Balanced Performance Index (Australian Profit Ranking index) that did not include FS-have resulted in an unfavorable genetic trend in FS. However, this unfavorable trend has stabilized since 2015, when FS was included in the Balanced Performance Index, and is expected to move in a favorable direction with selection on Balanced Performance Index or the Health Weighted Index. Doubling the reference population, particularly by incorporating international data for feed efficiency, has improved the reliability of the FS EBV. This could lead to increased genetic gain for feed efficiency in the Australian industry.