Data from: Multiple-trait genome-wide association study based on principal component analysis for residual covariance matrix

Given the drawbacks of implementing multivariate analysis for mapping multiple traits in genome-wide association study (GWAS), principal component analysis (PCA) has been widely used to generate independent ‘super traits’ from the original multivariate phenotypic traits for the univariate analysis. However, parameter estimates in this framework may not be the same as those from the joint analysis of all traits, leading to spurious linkage results. In this paper, we propose to perform the PCA for residual covariance matrix instead of the phenotypical covariance matrix, based on which multiple traits are transformed to a group of pseudo principal components. The PCA for residual covariance matrix allows analyzing each pseudo principal component separately. In addition, all parameter estimates are equivalent to those obtained from the joint multivariate analysis under a linear transformation. However, a fast least absolute shrinkage and selection operator (LASSO) for estimating the sparse oversaturated genetic model greatly reduces the computational costs of this procedure. Extensive simulations show statistical and computational efficiencies of the proposed method. We illustrate this method in a GWAS for 20 slaughtering traits and meat quality traits in beef cattle

Forward LASSO analysis for high-order interactions in genome-wide association study

Previous genome-wide association study (GWAS) focused on low-order interactions between pairwise single-nucleotide polymorphisms (SNPs) with significant main effects. Little is known how high-order interactions effect, especially one among the SNPs without main effects regulates quantitative traits. Within the frameworks of linear model and generalized linear model, the LASSO with coordinate descent step can be used to simultaneously analyze thousands and thousands of SNPs for normal and discrete traits. With consideration of high-order interactions among SNPs, a huge number of genetic effects make the LASSO failing to work under the presented condition of computation. Forward LASSO analysis is, therefore, proposed to shrink most of genetic effects to be zeros stage by stage. Simulation demonstrates that our proposed method could be used instead of the LASSO method for full model in mapping high-order interactions. Application of forward LASSO method is provided to GWAS for carcass traits and meat quality traits in beef cattle

Yang R: Forward LASSO analysis for high-order interactions in genome-wide association study. Brief Bioinform 2013. Jun 17. [Epub ahead of print

Abstract Previous genome-wide association study (GWAS) focused on low-order interactions between pairwise single-nucleotide polymorphisms (SNPs) with significant main effects. Little is known how high-order interactions effect, especially one among the SNPs without main effects regulates quantitative traits. Within the frameworks of linear model and generalized linear model, the LASSO with coordinate descent step can be used to simultaneously analyze thousands and thousands of SNPs for normal and discrete traits. With consideration of high-order interactions among SNPs, a huge number of genetic effects make the LASSO failing to work under the presented condition of computation. Forward LASSO analysis is, therefore, proposed to shrink most of genetic effects to be zeros stage by stage. Simulation demonstrates that our proposed method could be used instead of the LASSO method for full model in mapping high-order interactions. Application of forward LASSO method is provided to GWAS for carcass traits and meat quality traits in beef cattle

Generalized Linear Model for Mapping Discrete Trait Loci Implemented with LASSO Algorithm

Generalized estimating equation (GEE) algorithm under a heterogeneous residual variance model is an extension of the iteratively reweighted least squares (IRLS) method for continuous traits to discrete traits. In contrast to mixture model-based expectation-maximization (EM) algorithm, the GEE algorithm can well detect quantitative trait locus (QTL), especially large effect QTLs located in large marker intervals in the manner of high computing speed. Based on a single QTL model, however, the GEE algorithm has very limited statistical power to detect multiple QTLs because of ignoring other linked QTLs. In this study, the fast least absolute shrinkage and selection operator (LASSO) is derived for generalized linear model (GLM) with all possible link functions. Under a heterogeneous residual variance model, the LASSO for GLM is used to iteratively estimate the non-zero genetic effects of those loci over entire genome. The iteratively reweighted LASSO is therefore extended to mapping QTL for discrete traits, such as ordinal, binary, and Poisson traits. The simulated and real data analyses are conducted to demonstrate the efficiency of the proposed method to simultaneously identify multiple QTLs for binary and Poisson traits as examples

Genome-wide association studies using haplotypes and individual SNPs in Simmental cattle.

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components

Genome-wide detection of selective signatures in Simmental cattle

Artificial selection has greatly improved the beef production performance and changed its genetic basis. High-density SNP markers provide a way to track these changes and use selective signatures to search for the genes associated with artificial selection. In this study, we performed extended haplotype homozygosity (EHH) tests based on Illumina BovineSNP50 (54 K) Chip data from 942 Simmental cattle to identify significant core regions containing selective signatures, then verified the biological significance of these identified regions based on some commonly used bioinformatics analyses. A total of 224 regions over the whole genome in Simmental cattle showing the highest significance and containing some important functional genes, such as GHSR, TG and CANCNA2D1 were chosen. We also observed some significant terms in the enrichment analyses of second GO terms and KEGG pathways, indicating that these genes are associated with economically relevant cattle traits. This is the first detection of selection signature in Simmental cattle. Our findings significantly expand the selection signature map of the cattle genome, and identify functional candidate genes under positive selection for future genetic research

Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data

Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion–deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships

942

After quanlity control, this experimental population consists of 942 young Simmental bulls born between 2008 and 2011, which are originated from Ulgai, Xilingol league, Inner Mongolia of China. After weaning, the cattle were moved to Beijing Jinweifuren cattle farm and were fattened under the same feeding and management environment. Each individual was timely observed for growth and development traits until slaughtered from 16 to 18 months old. During the period of slaughter, carcass traits and meat quality traits were measured according to Institutional Meat Purchase Specifications for fresh beef guidelines

Genome-wide association study identifies loci and candidate genes for meat quality traits in Simmental beef cattle

Improving meat quality is the best way to enhance profitability and strengthen competitiveness in beef industry. Identification of genetic variants that control beef quality traits can help breeders design optimal breeding programs to achieve this goal. We carried out a genome-wide association study for meat quality traits in 1141 Simmental cattle using the Illumina Bovine HD 770K SNP array to identify the candidate genes and genomic regions associated with meat quality traits for beef cattle, including fat color, meat color, marbling score, longissimus muscle area, and shear force. In our study, we identified twenty significant single-nucleotide polymorphisms (SNPs) (p < 1.47 x 10(-6)) associated with these five meat quality traits. Notably, we observed several SNPs were in or near eleven genes which have been reported previously, including TMEM236, SORL1, TRDN, S100A10, AP2S1, KCTD16, LOC506594, DHX15, LAMA4, PREX1, and BRINP3. We identified a haplotype block on BTA13 containing five significant SNPs associated with fat color trait. We also found one of 19 SNPs was associated with multiple traits (shear force and longissimus muscle area) on BTA7. Our results offer valuable insights to further explore the potential mechanism of meat quality traits in Simmental beef cattle