奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状,随着平衡育种理念的提出和发展,繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展,全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

A candidate gene association study for nine economically important traits in Italian Holstein cattle

We genotyped 58 single nucleotide polymorphisms (SNPs) in 25 candidate genes in about 800 Italian Holstein sires. Fifty‐six (minor allele frequency >0.02) were used to evaluate their association with single traits: milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP), milk protein percentage (PP), milk somatic cell count (MSCC); and complex indexes: longevity, fertility and productivity–functionality type (PFT), using deregressed proofs, after adjustment for familial relatedness. Thirty‐two SNPs were significantly associated (proportion of false positives <0.05) with different traits: 16 with MSCC, 15 with PY, 14 with MY, 12 with PFT, eight with longevity, eight with FY, eight with PP, five with FP and two with fertility. In particular, a SNP in the promoter region of the PRLR gene was associated with eight of nine traits. DGAT1 polymorphisms were highly associated with FP and FY. Casein gene markers were associated with several traits, confirming the role of the casein gene cluster in affecting milk yield, milk quality and health traits. Other SNPs in genes located on chromosome 6 were associated with PY, PP, PFT, MY (PPARGC1A) and MSCC (KIT). This latter association may suggest a biological link between the degree of piebaldism in Holstein and immunological functions affecting somatic cell count and mastitis resistance. Other significant SNPs were in the ACACA, CRH, CXCR1, FASN, GH1, LEP, LGB (also known as PAEP), MFGE8, SRC, TG, THRSP and TPH1 genes. These results provide information that can complement QTL mapping and genome‐wide association studies in Holstein.     

Genome-wide association study reveals five nucleotide sequence variants for carcass traits in beef cattle

Genetic associations of nucleotide sequence variants with carcass traits in beef cattle were investigated using a genome-wide single nucleotide polymorphism (SNP) assay. Three hundred and thirteen Korean cattle were genotyped with the Illumina BovineSNP50 BeadChip, and 39,129 SNPs from 311 animals were analysed for each carcass phenotype after filtering by quality assurance. Five sequence markers were associated with one of the meat quantity or quality traits; rs109593638 on chromosome 3 with marbling score, rs109821175 on chromosome 11 and rs110862496 on chromosome 13 with backfat thickness (BFT), and rs110228023 on chromosome 6 and rs110201414 on chromosome 16 with eye muscle area (EMA) (P < 1.27 × 10(-6) , Bonferonni P < 0.05). The ss96319521 SNP, located within a gene with functions of muscle development, dishevelled homolog 1 (DVL1), would be a desirable candidate marker. Individuals with genotype CC at this gene appeared to have increased both EMA and carcass weight. Fine-mapping would be required to refine each of the five association signals shown in the current study for future application in marker-assisted selection for genetic improvement of beef quality and quantity.     

Accuracies of genomic prediction for twenty economically important traits in Chinese Simmental beef cattle

Genomic prediction has been widely utilized to estimate genomic breeding values (GEBVs) in farm animals. In this study, we conducted genomic prediction for 20 economically important traits including growth, carcass and meat quality traits in Chinese Simmental beef cattle. Five approaches (GBLUP, BayesA, BayesB, BayesCπ and BayesR) were used to estimate the genomic breeding values. The predictive accuracies ranged from 0.159 (lean meat percentage estimated by BayesCπ) to 0.518 (striploin weight estimated by BayesR). Moreover, we found that the average predictive accuracies across 20 traits were 0.361, 0.361, 0.367, 0.367 and 0.378, and the averaged regression coefficients were 0.89, 0.86, 0.89, 0.94 and 0.95 for GBLUP, BayesA, BayesB, BayesCπ and BayesR respectively. The genomic prediction accuracies were mostly moderate and high for growth and carcass traits, whereas meat quality traits showed relatively low accuracies. We concluded that Bayesian regression approaches, especially for BayesR and BayesCπ, were slightly superior to GBLUP for most traits. Increasing with the sizes of reference population, these two approaches are feasible for future application of genomic selection in Chinese beef cattle.     

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 × 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.     

Heritability and genetic correlation of body weight and Kleiber ratio in Limousin and Charolais beef cattle breeds

Reducing the environmental impact of livestock production is now indispensable and genetic selection can be of great support for this purpose. Measures that can identify high body growth at low maintenance costs in production animals are particularly useful since resources have been increasingly limited. Therefore, the goal of this study was to estimate genetic parameters for BW and Kleiber ratio (KR) in 210-day-old and 365-day-old Charolais and Limousin breeds. A database comprising animals born from 1999 to 2018 was used in a multitrait model applying Bayesian inference. The heritability for BW is high in Charolais (0.39 and 0.42 for BW210 and BW365, respectively) and moderate in Limousin (0.22), indicating possible genetic gains for BW in both breeds. The genetic variability of KR should also allow satisfactory genetic gains. In addition, the genetic correlation between BW and KR ranged from low to moderate. Thus, selection over KR should have no effects on BW, showing that high body growth can be obtained without changes in efficiency.     

Genome-wide association study for growth traits in Nelore cattle

The objective of this study was to investigate the association of single nucleotide polymorphisms (SNPs) with birth weight, weight gain from birth to weaning and from weaning to yearling, yearling height and cow weight in Nelore cattle. Data from 5064 animals participating in the DeltaGen and PAINT breeding programs were used. The animals were genotyped with a panel of 777 962 SNPs (Illumina BovineHD BeadChip) and 412 993 SNPs remained after quality control analysis of the genomic data. A genome-wide association study was performed using a single-step methodology. The analyses were processed with the BLUPF90 family of programs. When applied to a genome-wide association studies, the single-step GBLUP methodology is an iterative process that estimates weights for the SNPs. The weights of SNPs were included in all analyses by iteratively applying the single-step GBLUP methodology and repeated twice so that the effect of the SNP and the effect of the animal were recalculated in order to increase the weight of SNPs with large effects and to reduce the weight of those with small effects. The genome-wide association results are reported based on the proportion of variance explained by windows of 50 adjacent SNPs. Considering the two iterations, only windows with an additive genetic variance >1.5% were presented in the results. Associations were observed with birth weight on BTA 14, with weight gain from birth to weaning on BTA 5 and 29, with weight gain from weaning to yearling on BTA 11, and with yearling height on BTA 8, showing the genes TMEM68 (transmembrane protein 8B) associated with birth weight and yearling height, XKR4 (XK, Kell blood group complex subunit-related family, member 4) associated with birth weight, NPR2 (natriuretic peptide receptor B) associated with yearling height, and REG3G (regenerating islet-derived 3-gamma) associated with weight gain from weaning to yearling. These genes play an important role in feed intake, weight gain and the regulation of skeletal growth.     

Genome-wide association mapping of three important traits using bread wheat elite breeding populations

The exponential development of molecular markers enables a more effective study of the genetic architecture of traits of economic importance, like test weight in wheat (Triticum aestivum L.), for which a high value is desired by most end-users. The association mapping (AM) method now allows more precise exploration of the entire genome. AM requires populations with substantial genetic variability of the traits of interest. The breeding lines at the end of a selection cycle, characterized for numerous traits, represent a potentially useful population for AM studies. Using three elite line populations, selected by several breeders and genotyped with about 2,500 Diversity Arrays Technology markers, several associations were identified between these markers and test weight, grain yield and heading date. To minimize spurious associations, we compared the general linear model and mixed linear model (MLM), which adjust for population structure and kinship differently. The MLM model with the kinship matrix was the most efficient. Finally, elite lines from several breeding programs had sufficient genetic variability to allow for the mapping of several chromosomal regions involved in the variation of three important traits.     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状, 随着平衡育种理念的提出和发展, 繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展, 全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

Observations on the nursing behaviour of beef cows suckling Charolais cross calves

Three separate 48 h observations were made of the suckling behaviour of

• 1.(1) Hereford × Friesian cows under grazing conditions,
• 2.(2) Blue Grey cows under grazing conditions, and
• 3.(3) Blue Grey cows housed.

All the cows suckled Charolais cross calves. The mean suckling frequencies recorded were 5.5, 5.5 and 9.1 times per 24 h in the first, second and third study, respectively. A trend for older calves to be suckled less frequently was observed in all three studies, although this was not significant. Milk yield of the cow was considered a possible factor involved in determining suckling frequency. A relatively uniform distribution of suckling throughout the “day” was recorded in the second and third study, with three peaks of suckling activity apparent in the first study. The mean duration of suckling was 6.9, 10.7 and 10.6 min in the first, second and third study, respectively.     

Genome-wide association study for growth and fatness traits in Chinese Sujiang pigs

Growth and fatness traits are complex and economically important traits in the pig industry. The molecular basis underlying porcine growth and fatness traits remains largely unknown. To uncover genetic loci and candidate genes for these traits, we explored the GeneSeek GGP Porcine 80K SNP chip to perform a GWAS for seven growth and fatness traits in 365 individuals from the Sujiang pig, a recently developed breed in China. We identified two, 17, one and 11 SNPs surpassing the suggestively significant threshold (P < 1.86 × 10⁻⁵) for body weight, chest circumference, chest width and backfat thickness respectively. Of these SNPs, 20 represent novel genetic loci, and five and four SNPs were respectively associated with chest circumference and backfat thickness at a genome-wide significant threshold (P < 9.31 × 10⁻⁷). Eight SNPs had a pleiotropic effect on both chest circumference and backfat thickness. The most remarkable locus resided in a region between 72.95 and 76.27 Mb on pig chromosome 4, harboring a number of previously reported quantitative trait loci related to backfat deposition. In addition to two reported genes (PLAG1 and TAS2R38), we identified four genes including GABRB3, ZNF106, XKR4 and MGAM as novel candidates for body weight and backfat thickness at the mapped loci. Our findings provide insights into the genetic architecture of porcine growth and fatness traits and potential markers for selective breeding of Chinese Sujiang pigs.     

Genome-wide linkage and peak-wide association study of obesity-related quantitative traits in Caribbean Hispanics

Although obesity is more prevalent in Hispanics than non-Hispanic whites in the United States, little is known about the genetic etiology of the related traits in this population. To identify genetic loci influencing obesity in non-Mexican Hispanics, we performed a genome-wide linkage scan in 1,390 subjects from 100 Caribbean Hispanic families on six obesity-related quantitative traits: body mass index (BMI), body weight, waist circumference, waist-to-hip ratio, abdominal and average triceps skinfold thickness after adjusting for significant demographic and lifestyle factors. We then carried out an association analysis of the linkage peaks and the FTO gene in an independent community-based Hispanic subcohort (N = 652, 64% Caribbean Hispanics) from the Northern Manhattan Study. Evidence of linkage was strongest on 1q43 with multipoint LOD score of 2.45 (p = 0.0004) for body weight. Suggestive linkage evidence of LOD > 2.0 was also identified on 1q43 for BMI (LOD = 2.03), 14q32 for abdominal skinfold thickness (LOD = 2.17), 16p12 for BMI (LOD = 2.27) and weight (LOD = 2.26), and 16q23–24 for average triceps skinfold thickness (LOD = 2.32). In the association analysis of 6,440 single nucleotide polymorphisms (SNPs) under 1-LOD unit down regions of our linkage peaks on chromosome 1q43 and 16p12 as well as in the FTO gene, we found that two SNPs (rs6665519 and rs669231) on 1q43 and one FTO SNP (rs12447427) were significantly associated with BMI or body weight after adjustment for multiple testing. Our results suggest that in addition to FTO, multiple genetic loci, particularly those on 1q43 region, may contribute to the variations in obesity-related quantitative traits in Caribbean Hispanics.     

Genome-wide association study of reproductive traits in Nellore heifers using Bayesian inference

Background

An important goal of Zebu breeding programs is to improve reproductive performance. A major problem faced with the genetic improvement of reproductive traits is that recording the time for an animal to reach sexual maturity is costly. Another issue is that accurate estimates of breeding values are obtained only a long time after the young bulls have gone through selection. An alternative to overcome these problems is to use traits that are indicators of the reproductive efficiency of the herd and are easier to measure, such as age at first calving. Another problem is that heifers that have conceived once may fail to conceive in the next breeding season, which increases production costs. Thus, increasing heifer’s rebreeding rates should improve the economic efficiency of the herd. Response to selection for these traits tends to be slow, since they have a low heritability and phenotypic information is provided only later in the life of the animal. Genome-wide association studies (GWAS) are useful to investigate the genetic mechanisms that underlie these traits by identifying the genes and metabolic pathways involved.

Results

Data from 1853 females belonging to the Agricultural Jacarezinho LTDA were used. Genotyping was performed using the BovineHD BeadChip (777 962 single nucleotide polymorphisms (SNPs)) according to the protocol of Illumina - Infinium Assay II ® Multi-Sample HiScan with the unit SQ ™ System. After quality control, 305 348 SNPs were used for GWAS. Forty-two and 19 SNPs had a Bayes factor greater than 150 for heifer rebreeding and age at first calving, respectively. All significant SNPs for age at first calving were significant for heifer rebreeding. These 42 SNPs were next or within 35 genes that were distributed over 18 chromosomes and comprised 27 protein-encoding genes, six pseudogenes and two miscellaneous noncoding RNAs.

Conclusions

The use of Bayes factor to determine the significance of SNPs allowed us to identify two sets of 42 and 19 significant SNPs for heifer rebreeding and age at first calving, respectively, which explain 11.35 % and 6.42 % of their phenotypic variance, respectively. These SNPs provide relevant information to help elucidate which genes affect these traits.     

Genome-wide association study for calving performance using high-density genotypes in dairy and beef cattle

Background

Calving difficulty and perinatal mortality are prevalent in modern-day cattle production systems. It is well-established that there is a genetic component to both traits, yet little is known about their underlying genomic architecture, particularly in beef breeds. Therefore, we performed a genome-wide association study using high-density genotypes to elucidate the genomic architecture of these traits and to identify regions of the bovine genome associated with them.

Results

Genomic regions associated with calving difficulty (direct and maternal) and perinatal mortality were detected using two statistical approaches: (1) single-SNP (single nucleotide polymorphism) regression and (2) a Bayesian approach. Data included high-density genotypes on 770 Holstein-Friesian, 927 Charolais and 963 Limousin bulls. Several novel or previously identified genomic regions were detected but associations differed by breed. For example, two genomic associations, one each on chromosomes 18 and 2 explained 2.49 % and 3.13 % of the genetic variance in direct calving difficulty in the Holstein-Friesian and Charolais populations, respectively. Imputed Holstein-Friesian sequence data was used to refine the genomic regions responsible for significant associations. Several candidate genes on chromosome 18 were identified and four highly significant missense variants were detected within three of these genes (SIGLEC12, CTU1, and ZNF615). Nevertheless, only CTU1 contained a missense variant with a putative impact on direct calving difficulty based on SIFT (0.06) and Polyphen (0.95) scores. Using imputed sequence data, we refined a genomic region on chromosome 4 associated with maternal calving difficulty in the Holstein-Friesian population and found the strongest association with an intronic variant in the PCLO gene. A meta-analysis was performed across the three breeds for each calving performance trait to identify common variants associated with these traits in the three breeds. Our results suggest that a portion of the genetic variation in calving performance is common to all three breeds.

Conclusion

The genomic architecture of calving performance is complex and mainly influenced by many polymorphisms of small effect. We identified several associations of moderate effect size but the majority were breed-specific, indicating that breed-specific alleles exist for calving performance or that the linkage phase between genotyped allele and causal mutation varies between breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0126-4) contains supplementary material, which is available to authorized users.     

Genome-wide association study identifies candidate genes for male fertility traits in humans

Despite the fact that hundreds of genes are known to affect fertility in animal models, relatively little is known about genes that influence natural fertility in humans. To broadly survey genes contributing to variation in male fertility, we conducted a genome-wide association study (GWAS) of two fertility traits (family size and birth rate) in 269 married men who are members of a founder population of European descent that proscribes contraception and has large family sizes. Associations between ～250,000 autosomal SNPs and the fertility traits were examined. A total of 41 SNPs with p ≤ 1 × 10(-4) for either trait were taken forward to a validation study of 123 ethnically diverse men from Chicago who had previously undergone semen analyses. Nine (22%) of the SNPs associated with reduced fertility in the GWAS were also associated with one or more of the ten measures of reduced sperm quantity and/or function, yielding 27 associations with p values < 0.05 and seven with p values < 0.01 in the validation study. On the basis of 5,000 permutations of our data, the probabilities of observing this many or more small p values were 0.0014 and 5.6 × 10(-4), respectively. Among the nine associated loci, outstanding candidates for male fertility genes include USP8, an essential deubiquitinating enzyme that has a role in acrosome assembly; UBD and EPSTI1, which have potential roles in innate immunity; and LRRC32, which encodes a latent transforming growth factor β (TGF-β) receptor on regulatory T cells. We suggest that mutations in these genes that are more severe may account for some of the unexplained infertility (or subfertility) in the general population.     

Salt effects on functional traits in model and in economically important Lotus species

A common stress on plants is NaCl‐derived soil salinity. Genus Lotus comprises model and economically important species, which have been studied regarding physiological responses to salinity. Leaf area ratio (LAR), root length ratio (RLR) and their components, specific leaf area (SLA) and leaf mass fraction (LMF) and specific root length (SRL) and root mass fraction (RMF) might be affected by high soil salinity. We characterised L. tenuis, L. corniculatus, L. filicaulis, L. creticus, L. burtii and L. japonicus grown under different salt concentrations (0, 50, 100 and 150 mm NaCl) on the basis of SLA, LMF, SRL and RMF using PCA. We also assessed effects of different salt concentrations on LAR and RLR in each species, and explored whether changes in these traits provide fitness benefit. Salinity (150 mm NaCl) increased LAR in L. burtii and L. corniculatus, but not in the remaining species. The highest salt concentration caused a decrease of RLR in L. japonicus Gifu, but not in the remaining species. Changes in LAR and RLR would not be adaptive, according to adaptiveness analysis, with the exception of SLA changes in L. corniculatus. PCA revealed that under favourable conditions plants optimise surfaces for light and nutrient acquisition (SLA and SRL), whereas at higher salt concentrations they favour carbon allocation to leaves and roots (LMF and RMF) in detriment to their surfaces. PCA also showed that L. creticus subjected to saline treatment was distinguished from the remaining Lotus species. We suggest that augmented carbon partitioning to leaves and roots could constitute a salt‐alleviating mechanism through toxic ion dilution.     

Genome-wide association study reveals the genetic architecture of 27 agronomic traits in tomato

Tomato (Solanum lycopersicum) is a highly valuable fruit crop, and yield is one of the most important agronomic traits. However, the genetic architecture underlying tomato yield-related traits has not been fully addressed. Based on ∼4.4 million single nucleotide polymorphisms obtained from 605 diverse accessions, we performed a comprehensive genome-wide association study for 27 agronomic traits in tomato. A total of 239 significant associations corresponding to 129 loci, harboring many previously reported and additional genes related to vegetative and reproductive development, were identified, and these loci explained an average of ∼8.8% of the phenotypic variance. A total of 51 loci associated with 25 traits have been under selection during tomato domestication and improvement. Furthermore, a candidate gene, Sl-ACTIVATED MALATE TRANSPORTER15, that encodes an aluminum-activated malate transporter was functionally characterized and shown to act as a pivotal regulator of leaf stomata formation, thereby affecting photosynthesis and drought resistance. This study provides valuable information for tomato genetic research and breeding.

A large-scale genome-wide association study sheds light on genetic and genomic bases underlying important yield-related traits in tomato.