Eigen decomposition expedites longitudinal genome-wide association studies for milk production traits in Chinese Holstein

Identification of regulatory loci and candidate genes related to body weight traits in broilers based on different models

BACKGROUND

Growth traits are crucial for the economic viability in broiler production, as they significantly contribute to the cost of rearing. Maximizing body weight (BW) while minimizing feed intake is key to enhancing the efficiency of broiler breeding. Identifying the genetic architecture associated with BW trait is therefore a critical step in enhancing breeding strategies.

RESULTS

We conducted a genome-wide association study (GWAS) using two statistical approaches: single-trait GWAS and longitudinal GWAS. The study was performed on the BW trait at five developmental stages (72, 81, 89, 113, and 120 days) and mid-test metabolic weight (MWT) across four growth cycles. Transcriptome sequencing analysis was also included to investigate the differential expression of candidate genes identified through the GWAS models, particularly linked to BW and MWT traits. Utilizing the chicken 55K single nucleotide polymorphism (SNP) array, we identified 52,060 SNPs in the genomic data of 4,493 Wenchang chickens. The single-trait GWAS model revealed 42 BW-associated SNPs, corresponding to 18 potential genes. For MWT, 47 SNPs were associated, mapping to 31 candidate genes. The longitudinal GWAS model identified 34 BW-linked SNPs, annotated with 22 candidate genes, and 21 MWT-linked SNPs, annotated with 10 candidate genes. Notably, 16 SNPs on chromosome 4 were associated with both BW and MWT, located within the 73.08Mb-76.82Mb region. Nine genes were annotated from this region, including STIM2, SEL1L3, SEPSECS, LGI2, SOD3, KCNIP4, NCAPG, FAM184B, LDB2. Notably, there are 32 overlapping SNPs identified in both the single-trait and longitudinal GWAS models, suggesting consistent associations for both BW and MWT. These overlapping SNPs represent robust loci that may influence both traits across different statistical approaches. Transcriptome sequencing indicated differential expression of LDB2 and SEL1L3 between high and low BW groups.

CONCLUSION

Our study has uncovered novel candidate genes that are potentially involved in growth traits, providing valuable insights for broiler breeding. The identified SNPs and genes could serve as genetic markers for selecting broilers with improved growth efficiency, which may lead to more cost-effective and productive broiler farming.

Genomic insights into growth traits in German Black Pied cattle: a dual-purpose breed at risk

The German Black Pied cattle (DSN) is an endangered dual-purpose breed valued for its genetic diversity and high milk fat and protein content. However, due to competition with higher-yielding dairy breeds, the DSN population has declined, leading to its designation as an endangered breed. While previous research has focused on the milk production traits of DSN, this study aims to address meat traits to further understand the genetic determination of the dual-purpose characteristics of the breed. We conducted genome-wide association studies on 669 DSN bulls to identify genetic loci associated with birth weight, BW, and BW gain at different growth stages. Using imputed whole-genome sequencing data, we identified 14 quantitative trait loci across ten chromosomes. Significant associations were found for birth weight on chromosomes 5 and 18, for body weight at 3 weeks (BW3w) on chromosomes 3 and 16, for body weight at 7 months (BW7m) on chromosomes 3 and 10, and for body weight gain from birth or 3 weeks to 18 months (BWG0d-18m, BWG3w-18m) on chromosomes 4 and 7. Key positional candidate genes influencing muscle and fat tissue development included RERGL and LMO3 (identified for birth weight), MET and CAPZA2 (identified for BWG0d-18m) which are essential for skeletal muscle development and actin filament regulation, respectively, TLN2 (identified for BW7m), MYO1F and ADAMTS10 (identified for BWG3w-18m) which are critical for actin filament assembly, cytoskeletal function, and skeletal development, respectively. Candidate genes such as CPT2 (identified or BW3w) and VPS13C (identified or BW7m) are involved in lipid metabolism and mitochondrial function. Additionally, candidate genes such as IGSF3 (identified for BW7m), KLRC1 and members of the C-type lectin family (identified for birth weight) are associated with immune regulation, and thus, suggest a potential interplay between metabolism, immune function, and growth efficiency. These findings highlight the distinct genetic mechanisms underlying growth at various developmental stages, underscoring the importance of breed-specific genetic evaluations. The identified loci also overlap with previously reported loci for meat and production traits in other cattle breeds, underscoring their relevance and potential utility in DSN breeding strategies. This study provides a foundation for conservation and genomic breeding strategies to maintain the dual-purpose characteristics of DSN through optimising both meat and milk production.

Genome-Wide Association Study of Insertions and Deletions Identified Novel Loci Associated with Milk Production Traits in Dairy Cattle

Genome-wide association study (GWAS) have identified a large number of SNPs associated with milk production traits in dairy cattle. Behind SNPs, INDELs are the second most abundant genetic polymorphisms in the genome, which may exhibit an independent association with complex traits in humans and other species. However, there are no reports on GWASs of INDELs for milk production traits in dairy cattle. In this study, using imputed sequence data, we performed INDEL-based and SNP-based GWASs for milk production traits in a Holstein cattle population. We identified 58 unique significant INDELs for one or multiple traits. The majority of these INDELs are in considerable LD with nearby significant SNPs. However, through conditional association analysis, we identified nine INDELs which showed independent associations. Genomic annotations of these INDELs indicated some novel associated genes, i.e., TRNAG-CCC, EPPK1, PPM1K, PTDSS1, and mir-10163, which were not reported in previous SNP-based GWASs. Our findings suggest that INDEL-based GWASs could be valuable complement to SNP-based GWASs for milk production traits.

Genome-wide association study for milk yield, frame, and udder conformation traits of Gir dairy cattle

Genome-wide association studies (GWAS) are employed to identify genomic regions and candidate genes associated with several traits. The aim of this study was to perform a GWAS to identify causative variants and genes associated with milk yield, frame, and udder conformation traits in Gir dairy cattle. Body conformation traits were classified as "frame," and "udder" traits for this study. After genotyping imputation and quality control 42,105 polymorphisms were available for analyses and 24,489 cows with pedigree information had phenotypes. First, P-value was calculated based on the variance of the prediction error of the SNP-effects on the first iteration. After that, 2 more iterations were performed to carry out the weighted single-step genome-wide association methodology, performed using genomic moving windows defined based on linkage disequilibrium. The significant SNPs and top 10 windows explaining the highest percentage of additive genetic variance were selected and used for QTL and gene annotation. The variants identified in our work overlapped with QTLs from the animal QTL database on chromosomes 1 to 23, except for chromosome 4. The Gir breed is less studied than the Holstein breed and as such the animal QTL database is biased to Holstein results. Hence it is noteworthy that our GWAS had similarities with previously described QTLs. These previously known QTLs were related to milk yield, body height, rump angle, udder width, and udder depth. In total, 5 genes were annotated. Of these genes, FAM13A and CMSS1 had been previously related to bone and carcass weight in cattle.

Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens

The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Maize green leaf area index dynamics: genetic basis of a new secondary trait for grain yield in optimal and drought conditions

KEY MESSAGE

Green Leaf Area Index dynamics is a promising secondary trait for grain yield and drought tolerance. Multivariate GWAS is particularly well suited to identify the genetic determinants of the green leaf area index dynamics. Improvement of maize grain yield is impeded by important genotype-environment interactions, especially under drought conditions. The use of secondary traits, that are correlated with yield, more heritable and less prone to genotype-environment interactions, can increase breeding efficiency. Here, we studied the genetic basis of a new secondary trait: the green leaf area index (GLAI) dynamics over the maize life cycle. For this, we used an unmanned aerial vehicle to characterize the GLAI dynamics of a diverse panel in well-watered and water-deficient trials in two years. From the dynamics, we derived 24 traits (slopes, durations, areas under the curve), and showed that six of them were heritable traits representative of the panel diversity. To identify the genetic determinants of GLAI, we compared two genome-wide association approaches: a univariate (single-trait) method and a multivariate (multi-trait) method combining GLAI traits, grain yield, and precocity. The explicit modeling of correlation structure between secondary traits and grain yield in the multivariate mixed model led to 2.5 times more associations detected. A total of 475 quantitative trait loci (QTLs) were detected. The genetic architecture of GLAI traits appears less complex than that of yield with stronger-effect QTLs that are more stable between environments. We also showed that a subset of GLAI QTLs explains nearly one fifth of yield variability across a larger environmental network of 11 water-deficient trials. GLAI dynamics is a promising grain yield secondary trait in optimal and drought conditions, and the detected QTLs could help to increase breeding efficiency through a marker-assisted approach.

Longitudinal genome-wide association analysis using a single-step random regression model for height in Japanese Holstein cattle

Growth traits, such as body weight and height, are essential in the design of genetic improvement programs of dairy cattle due to their relationship with feeding efficiency, longevity, and health. We investigated genomic regions influencing height across growth stages in Japanese Holstein cattle using a single-step random regression model. We used 72,921 records from birth to 60 mo of age with 4,111 animals born between 2000 and 2016. The analysis included 1,410 genotyped animals with 35,319 single nucleotide polymorphisms, consisting of 883 females with records and 527 bulls, and 30,745 animals with pedigree information. A single genomic region at the 58.4 megabase pair on chromosome 18 was consistently identified across 6 age points of 10, 20, 30, 40, 50, and 60 mo after multiple testing corrections for the significance threshold. Twelve candidate genes, previously reported for longevity and gestation length, were found near the identified genomic region. Another location near the identified region was also previously associated with body conformation, fertility, and calving difficulty. Functional Gene Ontology enrichment analysis suggested that the candidate genes regulate dephosphorylation and phosphatase activity. Our findings show that further study of the identified candidate genes will contribute to a better understanding of the genetic basis of height in Japanese Holstein cattle.

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.

Longitudinal genome-wide association studies of milk production traits in Holstein cattle using whole-genome sequence data imputed from medium-density chip data

Longitudinal traits, such as milk production traits in dairy cattle, are featured by having phenotypic values at multiple time points, which change dynamically over time. In this study, we first imputed SNP chip (50-100K) data to whole-genome sequence (WGS) data in a Chinese Holstein population consisting of 6,470 cows. The imputation accuracies were 0.88 to 0.97 on average after quality control. We then performed longitudinal GWAS in this population based on a random regression test-day model using the imputed WGS data. The longitudinal GWAS revealed 16, 39, and 75 quantitative trait locus regions associated with milk yield, fat percentage, and protein percentage, respectively. We estimated the 95% confidence intervals (CI) for these quantitative trait locus regions using the logP drop method and identified 581 genes involved in these CI. Further, we focused on the CI that covered or overlapped with only 1 gene or the CI that contained an extremely significant top SNP. Twenty-eight candidate genes were identified in these CI. Most of them have been reported in the literature to be associated with milk production traits, such as DGAT1, HSF1, MGST1, GHR, ABCG2, ADCK5, and CSN1S1. Among the unreported novel genes, some also showed good potential as candidate genes, such as CCSER1, CUX2, SNTB1, RGS7, OSR2, and STK3, and are worth being further investigated. Our study provided not only new insights into the candidate genes for milk production traits, but also a general framework for longitudinal GWAS based on random regression test-day model using WGS data.

Exploring milk loss and variability during environmental perturbations across lactation stages as resilience indicators in Holstein cattle

Genetic selection for resilience is essential to improve the long-term sustainability of the dairy cattle industry, especially the ability of cows to maintain their level of production when exposed to environmental disturbances. Recording of daily milk yield provides an opportunity to develop resilience indicators based on milk losses and fluctuations in daily milk yield caused by environmental disturbances. In this context, our study aimed to explore milk loss traits and measures of variability in daily milk yield, including log-transformed standard deviation of milk deviations (Lnsd), lag-1 autocorrelation (Ra), and skewness of the deviations (Ske), as indicators of general resilience in dairy cows. The unperturbed dynamics of milk yield as well as milk loss were predicted using an iterative procedure of lactation curve modeling. Milk fluctuations were defined as a period of at least 10 successive days of negative deviations in which milk yield dropped at least once below 90% of the expected values. Genetic parameters of these indicators and their genetic correlation with economically important traits were estimated using single-trait and bivariate animal models and 8,935 lactations (after quality control) from 6,816 Chinese Holstein cows. In general, cows experienced an average of 3.73 environmental disturbances with a milk loss of 267 kg of milk per lactation. Each fluctuation lasted for 19.80 ± 11.46 days. Milk loss traits are heritable with heritability estimates ranging from 0.004 to 0.061. The heritabilities differed between Lnsd (0.135-0.250), Ra (0.008-0.058), and Ske (0.001-0.075), with the highest heritability estimate of 0.250 ± 0.020 for Lnsd when removing the first and last 10 days in milk in a lactation (Lnsd2). Based on moderate to high genetic correlations, lower Lnsd2 is associated with less milk losses, better reproductive performance, and lower disease incidence. These findings indicate that among the variables evaluated, Lnsd2 is the most promising indicator for breeding for improved resilience in Holstein cattle.

Genome-wide association analysis of milk production, somatic cell score, and body conformation traits in Holstein cows

Milk production and body conformation traits are critical economic traits for dairy cows. To understand the basic genetic structure for those traits, a genome wide association study was performed on milk yield, milk fat yield, milk fat percentage, milk protein yield, milk protein percentage, somatic cell score, body form composite index, daily capacity composite index, feed, and leg conformation traits, based on the Illumina Bovine HD100k BeadChip. A total of 57, 12 and 26 SNPs were found to be related to the milk production, somatic cell score and body conformation traits in the Holstein cattle. Genes with pleiotropic effect were also found in this study. Seven significant SNPs were associated with multi-traits and were located on the PLEC, PLEKHA5, TONSL, PTGER4, and LCORL genes. In addition, some important candidate genes, like GPAT3, CEBPB, AGO2, SLC37A1, and FNDC3B, were found to participate in fat metabolism or mammary gland development. These results can be used as candidate genes for milk production, somatic cell score, and body conformation traits of Holstein cows, and are helpful for further gene function analysis to improve milk production and quality.

Genome-Wide Association Study for Lactation Performance in the Early and Peak Stages of Lactation in Holstein Dairy Cows

Simple Summary

Although genome-wide association studies (GWAS) have been carried out within a variety of cattle breeds, they are mainly based on the accumulated 305-day lactation yield traits estimated by summing the test-day recorded every day during the lactation period, or combining the weekly or monthly test-day records by linear interpolation. Since the additive genetic variance for milk yield and composition changes during lactation, the genetic effects of QTL related to these traits are not constant during the lactation period. Therefore, a better understanding of the genetic architecture of milk production traits in different lactation stages (e.g., beginning, peak, and end stages of lactation) is needed. The aim of this study was to detect genomic loci associated with lactation performance during 9 to 50 days in milk (DIM) in Holstein dairy cows. Candidate genes identified for milk production traits showed contrasting results between the EARLY and PEAK stages of lactation. Based on the results of this study, it can be concluded that in any genomic study it should be taken into account that the genetic effects of genes related to the lactation performance are not constant during the lactation period.

Abstract

This study aimed to detect genomic loci associated with the lactation performance during 9 to 50 days in milk (DIM) in Holstein dairy cows. Daily milk yield (MY), fat yield (FY), and protein yield (PY) during 9 to 50 DIM were recorded on 134 multiparous Holstein dairy cows distributed in four research herds. Fat- and protein-corrected milk (FPCM), fat-corrected milk (FCM), and energy-corrected milk (ECM) were predicted. The records collected during 9 to 25 DIM were put into the early stage of lactation (EARLY) and those collected during 26 to 50 DIM were put into the peak stage of lactation (PEAK). Then, the mean of traits in each cow included in each lactation stage (EARLY and PEAK) were estimated and used as phenotypic observations for the genome-wide association study. The included animals were genotyped with the Illumina BovineHD Genotyping BeadChip (Illumina Inc., San Diego, CA, USA) for a total of 777,962 single nucleotide polymorphisms (SNPs). After quality control, 585,109 variants were analyzed using GEMMA software in a mixed linear model. Although there was no SNP associated with traits included at the 5% genome-wide significance threshold, 18 SNPs were identified to be associated with milk yield and composition at the suggestive genome-wide significance threshold. Candidate genes identified for milk production traits showed contrasting results between the EARLY and PEAK stages of lactation. This suggests that differential sets of candidate genes underlie the phenotypic expression of the considered traits in the EARLY and PEAK stages of lactation. Although further functional studies are needed to validate our findings in independent populations, it can be concluded that in any genomic study it should be taken into account that the genetic effects of genes related to the lactation performance are not constant during the lactation period.

Bayesian analysis of longitudinal traits in the Korea Association Resource (KARE) cohort

Various methodologies for the genetic analysis of longitudinal data have been proposed and applied to data from large-scale genome-wide association studies (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with traits of interest and to detect SNP-time interactions. We recently proposed a grid-based Bayesian mixed model for longitudinal genetic data and showed that our Bayesian method increased the statistical power compared to the corresponding univariate method and well detected SNP-time interactions. In this paper, we further analyze longitudinal obesity-related traits such as body mass index, hip circumference, waist circumference, and waist-hip ratio from Korea Association Resource data to evaluate the proposed Bayesian method. We first conducted GWAS analyses of cross-sectional traits and combined the results of GWAS analyses through a meta-analysis based on a trajectory model and a random-effects model. We then applied our Bayesian method to a subset of SNPs selected by meta-analysis to further discover SNPs associated with traits of interest and SNP-time interactions. The proposed Bayesian method identified several novel SNPs associated with longitudinal obesity-related traits, and almost 25% of the identified SNPs had significant p-values for SNP-time interactions.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

The Candidate Chromosomal Regions Responsible for Milk Yield of Cow: A GWAS Meta-Analysis

Milk yield (MY) is highly heritable and an economically important trait in dairy livestock species. To increase power to detect candidate genomic regions for this trait, we carried out a meta-analysis of genome-wide association studies (GWAS). In the present study, we identified 19 studies in PubMed for the meta-analysis. After review of the studies, 16 studies passed the filters for meta-analysis, and the number of chromosomes, detected markers and their positions, number of animals, and p-values were extracted from these studies and recorded. The final data set based on 16 GWAS studies had 353,698 cows and 3950 markers and was analyzed using METAL software. Our findings revealed 1712 significant (p-value < 2.5 × 10−6) genomic loci related to MY, with markers associated with MY found on all autosomes and sex chromosomes and the majority of them found on chromosome 14. Furthermore, gene ontology (GO) annotation was used to explore biological functions of the genes associated with MY; therefore, different regions of this chromosome may be suitable as genomic regions for further research into gene expression.

Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle

Simple Summary

Genome-wide association study (GWAS) has become the main approach for detecting functional genes that affects complex traits. For growth traits, the conventional GWAS method can only deal with the single-record traits observed at specific time points, rather than the longitudinal traits measured at multiple time points. Previous studies have reported the random regression model (RRM) for longitudinal data could overcome the limitation of the traditional GWAS model. Here, we present an association analysis based on RRM (GWAS-RRM) for 808 Chinese Simmental beef cattle at four stages of age. Ultimately, 37 significant single-nucleotide polymorphisms (SNPs) and several important candidate genes were screened to be associated with the body weight. Enrichment analysis showed these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. This study not only offers a further understanding of the genetic basis for growth traits in beef cattle, but also provides a robust analytics tool for longitudinal traits in various species.

Abstract

Body weight (BW) is an important longitudinal trait that directly described the growth gain of bovine in production. However, previous genome-wide association study (GWAS) mainly focused on the single-record traits, with less attention paid to longitudinal traits. Compared with traditional GWAS models, the association studies based on the random regression model (GWAS-RRM) have better performance in the control of the false positive rate through considering time-stage effects. In this study, the BW trait data were collected from 808 Chinese Simmental beef cattle aged 0, 6, 12, and 18 months, then we performed a GWAS-RRM to fit the time-varied SNP effect. The results showed a total of 37 significant SNPs were associated with BW. Gene functional annotation and enrichment analysis indicated FGF4, ANGPT4, PLA2G4A, and ITGA5 were promising candidate genes for BW. Moreover, these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. These findings will provide prior molecular information for bovine gene-based selection, as well as facilitate the extensive application of GWAS-RRM in domestic animals.

Hierarchical mixed-model expedites genome-wide longitudinal association analysis

A hierarchical random regression model (Hi-RRM) was extended into a genome-wide association analysis for longitudinal data, which significantly reduced the dimensionality of repeated measurements. The Hi-RRM first modeled the phenotypic trajectory of each individual using a RRM and then associated phenotypic regressions with genetic markers using a multivariate mixed model (mvLMM). By spectral decomposition of genomic relationship and regression covariance matrices, the mvLMM was transformed into a multiple linear regression, which improved computing efficiency while implementing mvLMM associations in efficient mixed-model association expedited (EMMAX). Compared with the existing RRM-based association analyses, the statistical utility of Hi-RRM was demonstrated by simulation experiments. The method proposed here was also applied to find the quantitative trait nucleotides controlling the growth pattern of egg weights in poultry data.

Integrative analysis of genome-wide DNA methylation and gene expression profiles reveals important epigenetic genes related to milk production traits in dairy cattle

Epigenetic modification plays a critical role in establishing and maintaining cell differentiation, embryo development, tumorigenesis and many complex diseases. However, little is known about the epigenetic regulatory mechanisms for milk production in dairy cattle. Here, we conducted an epigenome-wide study, together with gene expression profiles to identify important epigenetic candidate genes related to the milk production traits in dairy cattle. Whole-genome bisulphite sequencing and RNA sequencing were employed to detect differentially methylated genes (DMG) and differentially expressed genes (DEG) in blood samples in dry period and lactation period between two groups of cows with extremely high and low milk production performance. A total of 10,877 and 6,617 differentially methylated regions were identified between the two groups in the two periods, which corresponded to 3,601 and 2,802 DMGs, respectively. Furthermore, 156 DEGs overlap with DMGs in comparison of the two groups, and 131 DEGs overlap with DMGs in comparison of the two periods. By integrating methylome, transcriptome and GWAS data, some potential candidate genes for milk production traits in dairy cattle were suggested, such as DOCK1, PTK2 and PIK3R1. Our studies may contribute to a better understanding of epigenetic modification on milk production traits of dairy cattle.

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.

Bayesian Analysis of the Association between Casein Complex Haplotype Variants and Milk Yield, Composition, and Curve Shape Parameters in Murciano-Granadina Goats

Considering casein haplotype variants rather than SNPs may maximize the understanding of heritable mechanisms and their implication on the expression of functional traits related to milk production. Effects of casein complex haplotypes on milk yield, milk composition, and curve shape parameters were used using a Bayesian inference for ANOVA. We identified 48 single nucleotide polymorphisms (SNPs) present in the casein complex of 159 unrelated individuals of diverse ancestry, which were organized into 86 haplotypes. The Ali and Schaeffer model was chosen as the best fitting model for milk yield (Kg), protein, fat, dry matter, and lactose (%), while parabolic yield-density was chosen as the best fitting model for somatic cells count (SCC × 103 sc/mL). Peak and persistence for all traits were computed respectively. Statistically significant differences (p < 0.05) were found for milk yield and components. However, no significant difference was found for any curve shape parameter except for protein percentage peak. Those haplotypes for which higher milk yields were reported were the ones that had higher percentages for protein, fat, dry matter, and lactose, while the opposite trend was described by somatic cells counts. Conclusively, casein complex haplotypes can be considered in selection strategies for economically important traits in dairy goats.

Integrating High-Throughput Phenotyping and Statistical Genomic Methods to Genetically Improve Longitudinal Traits in Crops

The rapid development of remote sensing in agronomic research allows the dynamic nature of longitudinal traits to be adequately described, which may enhance the genetic improvement of crop efficiency. For traits such as light interception, biomass accumulation, and responses to stressors, the data generated by the various high-throughput phenotyping (HTP) methods requires adequate statistical techniques to evaluate phenotypic records throughout time. As a consequence, information about plant functioning and activation of genes, as well as the interaction of gene networks at different stages of plant development and in response to environmental stimulus can be exploited. In this review, we outline the current analytical approaches in quantitative genetics that are applied to longitudinal traits in crops throughout development, describe the advantages and pitfalls of each approach, and indicate future research directions and opportunities.

Genome-wide association study for genotype by lactation stage interaction of milk production traits in dairy cattle

Substantial evidence demonstrates that the genetic background of milk production traits changes during lactation. However, most GWAS for milk production traits assume that genetic effects are constant during lactation and therefore might miss those quantitative trait loci (QTL) whose effects change during lactation. The GWAS for genotype by lactation stage interaction are aimed at explicitly detecting the QTL whose effects change during lactation. The purpose of this study was to perform GWAS for genotype by lactation stage interaction for milk yield, lactose yield, lactose content, fat yield, fat content, protein yield, and somatic cell score to detect QTL with changing effects during lactation. For this study, 19,286 test-day records of 1,800 first-parity Dutch Holstein cows were available and cows were genotyped using a 50K SNP panel. A total of 7 genomic regions with effects that change during lactation were detected in the GWAS for genotype by lactation stage interaction. Two regions on Bos taurus autosome (BTA)14 and BTA19 were also significant based on a GWAS that assumed constant genetic effects during lactation. Five regions on BTA4, BTA10, BTA11, BTA16, and BTA23 were only significant in the GWAS for genotype by lactation stage interaction. The biological mechanisms that cause these changes in genetic effects are still unknown, but negative energy balance and effects of pregnancy may play a role. These findings increase our understanding of the genetic background of lactation and may contribute to the development of better management indicators based on milk composition.

Genome-wide association for milk production and lactation curve parameters in Holstein dairy cows

The aim of this study was to identify genomic regions associated with 305‐day milk yield and lactation curve parameters on primiparous (n = 9,910) and multiparous (n = 11,158) Holstein cows. The SNP solutions were estimated using a weighted single‐step genomic BLUP approach and imputed high‐density panel (777k) genotypes. The proportion of genetic variance explained by windows of 50 consecutive SNP (with an average of 165 Kb) was calculated, and regions that accounted for more than 0.50% of the variance were used to search for candidate genes. Estimated heritabilities were 0.37, 0.34, 0.17, 0.12, 0.30 and 0.19, respectively, for 305‐day milk yield, peak yield, peak time, ramp, scale and decay for primiparous cows. Genetic correlations of 305‐day milk yield with peak yield, peak time, ramp, scale and decay in primiparous cows were 0.99, 0.63, 0.20, 0.97 and −0.52, respectively. The results identified three windows on BTA14 associated with 305‐day milk yield and the parameters of lactation curve in primi‐ and multiparous cows. Previously proposed candidate genes for milk yield supported by this work include GRINA, CYHR1, FOXH1, TONSL, PPP1R16A, ARHGAP39, MAF1, OPLAH and MROH1, whereas newly identified candidate genes are MIR2308, ZNF7, ZNF34, SLURP1, MAFA and KIFC2 (BTA14). The protein lipidation biological process term, which plays a key role in controlling protein localization and function, was identified as the most important term enriched by the identified genes.

Invited review: Advances and applications of random regression models: From quantitative genetics to genomics

An important goal in animal breeding is to improve longitudinal traits; that is, traits recorded multiple times during an individual's lifetime or physiological cycle. Longitudinal traits were first genetically evaluated based on accumulated phenotypic expression, phenotypic expression at specific time points, or repeatability models. Until now, the genetic evaluation of longitudinal traits has mainly focused on using random regression models (RRM). Random regression models enable fitting random genetic and environmental effects over time, which results in higher accuracy of estimated breeding values compared with other statistical approaches. In addition, RRM provide insights about temporal variation of biological processes and the physiological implications underlying the studied traits. Despite the fact that genomic information has substantially contributed to increase the rates of genetic progress for a variety of economically important traits in several livestock species, less attention has been given to longitudinal traits in recent years. However, including genomic information to evaluate longitudinal traits using RRM is a feasible alternative to yield more accurate selection and culling decisions, because selection of young animals may be based on the complete pattern of the production curve with higher accuracy compared with the use of traditional parent average (i.e., without genomic information). Moreover, RRM can be used to estimate SNP effects over time in genome-wide association studies. Thus, by analyzing marker associations over time, regions with higher effects at specific points in time are more likely to be identified. Despite the advances in applications of RRM in genetic evaluations, more research is needed to successfully combine RRM and genomic information. Future research should provide a better understanding of the temporal variation of biological processes and their physiological implications underlying the longitudinal traits.

Genome-wide association studies for genetic effects that change during lactation in dairy cattle

Genetic effects on milk production traits in dairy cattle might change during lactation. However, most genome-wide association studies (GWAS) for milk production traits assume that genetic effects are constant during lactation. This assumption might lead to missing these quantitative trait loci (QTL) whose effects change during lactation. This study aimed to screen the whole genome specifically for QTL whose effects change during lactation. For this purpose, 4 different GWAS approaches were performed using test-day milk protein content records: (1) separate GWAS for specific lactation stages, (2) GWAS for estimated Wilmink lactation curve parameters, (3) a GWAS using a repeatability model where SNP effects are assumed constant during lactation, and (4) a GWAS for genotype by lactation stage interaction using a repeatability model and accounting for changing genetic effects during lactation. Separate GWAS for specific lactation stages suggested that the detection power greatly differs between lactation stages and that genetic effects of some QTL change during lactation. The GWAS for estimated Wilmink lactation curve parameters detected many chromosomal regions for Wilmink parameter a (protein content level), whereas 2 regions for Wilmink parameter b (decrease in protein content toward nadir) and no regions for Wilmink parameter c (increase in protein content after nadir) were detected. Twenty chromosomal regions were detected with effects on milk protein content; however, there was no evidence that their effects changed during lactation. For 5 chromosomal regions located on chromosomes 3, 9, 10, 14, and 27, significant evidence was observed for a genotype by lactation stage interaction and thus their effects on milk protein content changed during lactation. Three of these 5 regions were only identified using a GWAS for genotype by lactation stage interaction. Our study demonstrated that GWAS for genotype by lactation stage interaction offers new possibilities to identify QTL involved in milk protein content. The performed approaches can be applied to other milk production traits. Identification of QTL whose genetic effects change during lactation will help elucidate the genetic and biological background of milk production.

Efficient multivariate analysis algorithms for longitudinal genome-wide association studies

Motivation

Current dynamic phenotyping system introduces time as an extra dimension to genome-wide association studies (GWAS), which helps to explore the mechanism of dynamical genetic control for complex longitudinal traits. However, existing methods for longitudinal GWAS either ignore the covariance among observations of different time points or encounter computational efficiency issues.

Results

We herein developed efficient genome-wide multivariate association algorithms for longitudinal data. In contrast to existing univariate linear mixed model analyses, the proposed method has improved statistic power for association detection and computational speed. In addition, the new method can analyze unbalanced longitudinal data with thousands of individuals and more than ten thousand records within a few hours. The corresponding time for balanced longitudinal data is just a few minutes.

Availability and implementation

A software package to implement the efficient algorithm named GMA (https://github.com/chaoning/GMA) is available freely for interested users in relevant fields.

Supplementary information

Supplementary data are available at Bioinformatics online.

Short communication: Replication of genome-wide association studies for milk production traits in Chinese Holstein by an efficient rotated linear mixed model

Milk is regarded as an important nutrient for humans, and Chinese Holstein cows provide high-quality milk for billions of Chinese people. Therefore, detecting quantitative trait nucleotides (QTN) or candidate genes for milk production traits in Chinese Holstein is important. In this study, we performed genome-wide association studies (GWAS) in a Chinese Holstein population of 6,675 cows and 71,633 SNP using deregressed proofs (DRP) as phenotypes to replicate our previous study in a population of 1,815 cows and 39,163 SNP using estimated breeding values (EBV) as phenotypes. The associations between 3 milk production traits-milk yield (MY), fat percentage (FP), and protein percentage (PP)-and the SNP were determined by using an efficient rotated linear mixed model, which benefits from linear transformations of genomic estimated values and Eigen decomposition of the genomic relationship matrix algorithm. In total, we detected 94 SNP that were significantly associated with one or more milk production traits, including 7 SNP for MY, 76 for FP, and 36 for PP; 87% of these SNP were distributed across Bos taurus autosomes 14 and 20. In total, 83 SNP were found to be located within the reported quantitative trait loci (QTL) regions, and one novel segment (between 1.41 and 1.49 Mb) on chromosome 14 was significantly associated with FP, which could be an important candidate QTL region. In addition, the detected intervals were narrowed down from the reported regions harboring causal variants. The top significant SNP for the 3 traits was ARS-BFGL-NGS-4939, which is located within the DGAT1 gene. Five detected genes (CYHR1, FOXH1, OPLAH, PLEC, VPS28) have effects on all 3 traits. Our study provides a suite of QTN, candidate genes, and a novel QTL associated with milk production traits, and thus forms a solid basis for genomic selection and molecular breeding for milk production traits in Chinese Holstein.