Multi-Trait GWAS and New Candidate Genes Annotation for Growth Curve Parameters in Brahman 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.

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.

The growth curve determination and economic trait correlation for Qinchuan bull population

Carcass weight, as a measure of meat yield, and body measurements are directly correlated traits in livestock. However, longitudinally collected phenotype records of local breeds are not comprehensive. The research was performed on Qinchuan bull population to understand their growth and development, and data from Qinchuan bull that was weighed and measured at birth, 6, 12, 18, and 24 months of age was analyzed. Furthermore, Logistic, Brody, Gompertz, and Bertallanffy were used to fit the growth curves for weight and body size traits. The results showed that the four curve models have good fitting degrees for the weight and body size (R2 > 0.99), and the Bertallanffy model exhibited a good fit to the measured data of body weight, and the model estimated the inflection point of body weight as (5.43 months of age, 122.01 kg). Particularly, the limited mature body weight can reach 557.8 kg by the Brody model. Body weight was significantly positively correlated with body height, hip height, body length, chest circumference, abdominal girth, and calf girth (p < 0.0001), and the correlation between body weight and body length was the highest (r = 0.975). The regression equation predicting body weight was Y = -275.691 + 3.28 X3 + 1.311 X4 - 0.397 X5.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis

BACKGROUND

The genetic locus responsible for duck body size has been fully explained before, but the growth trait-related genetic basis is still waiting to be explored. For example, the genetic site related to growth rate, an important economic trait affecting marketing weight and feeding cost, is still unclear. Here, we performed genome wide association study (GWAS) to identify growth rate-associated genes and mutations.

RESULT

In the current study, the body weight data of 358 ducks were recorded every 10 days from hatching to 120 days of age. According to the growth curve, we evaluated the relative and absolute growth rates (RGR and AGR) of 5 stages during the early rapid growth period. GWAS results for RGRs identified 31 significant SNPs on autosomes, and these SNPs were annotated by 24 protein-coding genes. Fourteen autosomal SNPs were significantly associated with AGRs. In addition, 4 shared significant SNPs were identified as having an association with both AGR and RGR, which were Chr2: 11483045 C>T, Chr2: 13750217 G>A, Chr2: 42508231 G>A and Chr2: 43644612 C>T. Among them, Chr2: 11483045 C>T, Chr2: 42508231 G>A, and Chr2: 43644612 C>T were annotated by ASAP1, LYN and CABYR, respectively. ASAP1 and LYN have already been proven to play roles in the growth and development of other species. In addition, we genotyped every duck using the most significant SNP (Chr2: 42508231 G>A) and compared the growth rate difference among each genotype population. The results showed that the growth rates of individuals carrying the Chr2: 42508231 A allele were significantly lower than those without this allele. Moreover, the results of the Mendelian randomization (MR) analysis supported the idea that the growth rate and birth weight had a causal effect on the adult body weight, with the growth rate having a greater effect size.

CONCLUSION

In this study, 41 SNPs significantly related to growth rate were identified. In addition, we considered that the ASAP1 and LYN genes are essential candidate genes affecting the duck growth rate. The growth rate also showed the potential to be used as a reliable predictor of adult weight, providing a theoretical reference for preselection.

Early Growth and Development and Nonlinear Model Fitting Analysis of Ashidan Yak

Understanding animal growth plays an important role in improving animal genetics and breeding. In order to explore the early growth and development law of Ashidan yak, the body weight (BW), wither height (WH), body oblique length (BL) and chest girth (CG) of 260 female Ashidan yaks were measured. These individuals grew under grazing conditions, and growth traits were measured at 6, 12, 18 and 30 months of age. Then the absolute growth and relative growth of Ashidan yak were calculated, and five nonlinear models (Logistic model, Gompertz model, Brody model, von Bertalanffy model and Richards model) were used to fit the growth curve of Ashidan yak. The fitting effect of the model was evaluated according to MSE, AIC and BIC. The results showed that the growth rate of Ashidan yak was the fastest from 12 to 18 months old, and the growth was slow or even stagnant from 6 to 12 months old. The AIC and BIC values of the Richards model were the lowest among the five models, with an AIC value of 4543.98 and a BIC value of 4563.19. The Richards model estimated body weight at 155.642 kg. In summary, the growth rate of female Ashidan yak changes with the seasons, growing faster in warm seasons and slower in cold seasons. Richards model is the best model to describe the growth curve of female Ashidan yak in five nonlinear models.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

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

Background

Genotypic information produced from single nucleotide polymorphism (SNP) arrays has routinely been used to identify genomic regions associated with complex traits in beef and dairy cattle. Herein, we assembled a dataset consisting of 15,815 Red Angus beef cattle distributed across the continental U.S. and a union set of 836,118 imputed SNPs to conduct genome-wide association analyses (GWAA) for growth traits using univariate linear mixed models (LMM); including birth weight, weaning weight, and yearling weight. Genomic relationship matrix heritability estimates were produced for all growth traits, and genotype-by-environment (GxE) interactions were investigated.

Results

Moderate to high heritabilities with small standard errors were estimated for birth weight (0.51 ± 0.01), weaning weight (0.25 ± 0.01), and yearling weight (0.42 ± 0.01). GWAA revealed 12 pleiotropic QTL (BTA6, BTA14, BTA20) influencing Red Angus birth weight, weaning weight, and yearling weight which met a nominal significance threshold (P ≤ 1e-05) for polygenic traits using 836K imputed SNPs. Moreover, positional candidate genes associated with Red Angus growth traits in this study (i.e., LCORL, LOC782905, NCAPG, HERC6, FAM184B, SLIT2, MMRN1, KCNIP4, CCSER1, GRID2, ARRDC3, PLAG1, IMPAD1, NSMAF, PENK, LOC112449660, MOS, SH3PXD2B, STC2, CPEB4) were also previously associated with feed efficiency, growth, and carcass traits in beef cattle. Collectively, 14 significant GxE interactions were also detected, but were less consistent among the investigated traits at a nominal significance threshold (P ≤ 1e-05); with one pleiotropic GxE interaction detected on BTA28 (24 Mb) for Red Angus weaning weight and yearling weight.

Conclusions

Sixteen well-supported QTL regions detected from the GWAA and GxE GWAA for growth traits (birth weight, weaning weight, yearling weight) in U.S. Red Angus cattle were found to be pleiotropic. Twelve of these pleiotropic QTL were also identified in previous studies focusing on feed efficiency and growth traits in multiple beef breeds and/or their composites. In agreement with other beef cattle GxE studies our results implicate the role of vasodilation, metabolism, and the nervous system in the genetic sensitivity to environmental stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08667-6.

Genomic Prediction of Green Fraction Dynamics in Soybean Using Unmanned Aerial Vehicles Observations

With the widespread use of high-throughput phenotyping systems, growth process data are expected to become more easily available. By applying genomic prediction to growth data, it will be possible to predict the growth of untested genotypes. Predicting the growth process will be useful for crop breeding, as variability in the growth process has a significant impact on the management of plant cultivation. However, the integration of growth modeling and genomic prediction has yet to be studied in depth. In this study, we implemented new prediction models to propose a novel growth prediction scheme. Phenotype data of 198 soybean germplasm genotypes were acquired for 3 years in experimental fields in Tottori, Japan. The longitudinal changes in the green fractions were measured using UAV remote sensing. Then, a dynamic model was fitted to the green fraction to extract the dynamic characteristics of the green fraction as five parameters. Using the estimated growth parameters, we developed models for genomic prediction of the growth process and tested whether the inclusion of the dynamic model contributed to better prediction of growth. Our proposed models consist of two steps: first, predicting the parameters of the dynamics model with genomic prediction, and then substituting the predicted values for the parameters of the dynamics model. By evaluating the heritability of the growth parameters, the dynamic model was able to effectively extract genetic diversity in the growth characteristics of the green fraction. In addition, the proposed prediction model showed higher prediction accuracy than conventional genomic prediction models, especially when the future growth of the test population is a prediction target given the observed values in the first half of growth as training data. This indicates that our model was able to successfully combine information from the early growth period with phenotypic data from the training population for prediction. This prediction method could be applied to selection at an early growth stage in crop breeding, and could reduce the cost and time of field trials.

Genome-wide association study reveals candidate genes relevant to body weight in female turkeys (Meleagris gallopavo)

The underlying genetic mechanisms affecting turkey growth traits have not been widely investigated. Genome-wide association studies (GWAS) is a powerful approach to identify candidate regions associated with complex phenotypes and diseases in livestock. In the present study, we performed GWAS to identify regions associated with 18-week body weight in a turkey population. The data included body weight observations for 24,989 female turkeys genotyped based on a 65K SNP panel. The analysis was carried out using a univariate mixed linear model with hatch-week-year and the 2 top principal components fitted as fixed effects and the accumulated polygenic effect of all markers captured by the genomic relationship matrix as random. Thirty-three significant markers were observed on 1, 2, 3, 4, 7 and 12 chromosomes, while 26 showed strong linkage disequilibrium extending up to 410 kb. These significant markers were mapped to 37 genes, of which 13 were novel. Interestingly, many of the investigated genes are known to be involved in growth and body weight. For instance, genes AKR1D1, PARP12, BOC, NCOA1, ADCY3 and CHCHD7 regulate growth, body weight, metabolism, digestion, bile acid biosynthetic and development of muscle cells. In summary, the results of our study revealed novel candidate genomic regions and candidate genes that could be managed within a turkey breeding program and adapted in fine mapping of quantitative trait loci to enhance genetic improvement in this species.

Association Between Natriuretic Peptide Receptor 2 (NPR2) RS208158047 Polymorphism and Fattening Performance of Young Bulls

The objective of this study was to determine fattening performance data for Charolais, Limousin and Blonde d’Aquitaine beef cattle and associate these data with NPR2 gene 8:g.59961937 T>C (rs208158047) mutation. Experiments were conducted with 176 beef cattle (77 Charolais, 66 Limousin and 33 Blonde d’Aquitaine) at nine months of age. Experiments lasted for 9 months and animals were slaughtered at the age of 18 months. Cattle body weights were determined at four different periods: beginning of fattening (d0), 60th day of fattening (d60), 120th day of fattening (d120) and at the end of fattening (sw). In terms of rs208158047 mutation of Charolais, Limousin and Blonde d’Aquitaine breeds, TT and CT genotypes were identified, and CC genotype was not encountered. The association of average daily gain (ADG) in d0-d60, d0-d120 and d0-sw periods with the genotypes of rs208158047 mutation was found to be significant (P<0.05). Greater ADGs were observed in rs208158047-CT genotypes compared to rs208158047-TT genotypes. These results indicate that the selection of bovine NPR2 gene could be used to ensure the breeding direction for growth related traits of the beef cattle.

Development of DNA markers for improvement of meat quality in a Japanese Black cattle population in Hyogo Prefecture

The polymorphisms associated with economic traits in livestock animals provide useful information as genetic indicators for breeding improvement. Over the last two decades, several DNA markers have been developed in Japanese Black cattle; however, the effect of these markers differs across populations due to differences in their genetic structures and backgrounds. As such, there is a need to verify the effectiveness of these markers in each population. This review summarizes the effectiveness of previously reported markers on carcass traits and the development of novel DNA markers in a Japanese Black cattle population in Hyogo Prefecture. As result of genome wide association studies and resequencing analyses, two novel significant markers associated with meat quality-related traits (beef marbling and fatty acid composition) were developed. These findings will lead to the identification of responsible genes and polymorphisms and contribute to the development of novel DNA markers for numerous traits in various cattle populations.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Genome-Wide Association Studies for Growth Curves in Meat Rabbits Through the Single-Step Nonlinear Mixed Model

Growth is a complex trait with moderate to high heritability in livestock and must be described by the longitudinal data measured over multiple time points. Therefore, the used phenotype in genome-wide association studies (GWAS) of growth traits could be either the measures at the preselected time point or the fitted parameters of whole growth trajectory. A promising alternative approach was recently proposed that combined the fitting of growth curves and estimation of single-nucleotide polymorphism (SNP) effects into single-step nonlinear mixed model (NMM). In this study, we collected the body weights at 35, 42, 49, 56, 63, 70, and 84 days of age for 401 animals in a crossbred population of meat rabbits and compared five fitting models of growth curves (Logistic, Gompertz, Brody, Von Bertalanffy, and Richards). The logistic model was preferably selected and subjected to GWAS using the approach of single-step NMM, which was based on 87,704 genome-wide SNPs. A total of 45 significant SNPs distributed on five chromosomes were found to simultaneously affect the two growth parameters of mature weight (A) and maturity rate (K). However, no SNP was found to be independently associated with either A or K. Seven positional genes, including KCNIP4, GBA3, PPARGC1A, LDB2, SHISA3, GNA13, and FGF10, were suggested to be candidates affecting growth performances in meat rabbits. To the best of our knowledge, this is the first report of GWAS based on single-step NMM for longitudinal traits in rabbits, which also revealed the genetic architecture of growth traits that are helpful in implementing genome selection.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

Multi-trait GWAS using imputed high-density genotypes from whole-genome sequencing identifies genes associated with body traits in Nile tilapia

Background

Body traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits.

Results

A multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively.

Conclusions

The high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07341-z.

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.

Identification of new genes and quantitative trait locis associated with growth curve parameters in F2 chicken population using genome-wide association study

The markers which are correlated with the growth curve parameters help in understanding the characteristics of individual growth during the rearing of livestock. This study aimed to identify a set of biomarkers through a GWAS for growth curve parameters in crossbred chickens using the Illumnia 60K chicken SNP Beadchip. Growth data were collected from a total of 301 birds from cross of a broiler line and native chickens. Using the Gompertz-Laird model, two growth curve parameters, the instantaneous growth rate per day (L) and the coefficient of relative growth or maturing index (k), were estimated. The L and k were used to estimate five derived parameters, namely asymptotic (mature) body weight, body weight at inflection point, age at the inflection point, average growth rate and maximum growth rate. These parameters were considered as phenotypic values in the GWAS based on generalized linear models. The results of the GWAS indicated 21 significant markers, which were located near or within 46 genes. A number of these genes, such as GH, RET, GRB14, FTSJ3 and CCK, are important for growth and meat quality in chickens, and some of them are growth related in other species such as sheep and cattle (GPI, XIRP2, GALNTL6, BMS1, THSD4, TRHDE, SHISA9, ACSL6 and DYNC1LI2). The other genes are associated with developmental biological pathways. These genes are particuarly related to body weight, average daily gain and growth QTL. The results of this study can shed light on the genetic mechanism of biological functions of growth factors in broiler chickens, which is useful for developing management practices and accelerating genetic progress in breeding programs.

Connecting mathematical models to genomes: joint estimation of model parameters and genome-wide marker effects on these parameters

MOTIVATION

Parameters of mathematical models used in biology may be genotype-specific and regarded as new traits. Therefore, an accurate estimation of these parameters and the association mapping on the estimated parameters can lead to important findings regarding the genetic architecture of biological processes. In this study, a statistical framework for a joint analysis (JA) of model parameters and genome-wide marker effects on these parameters was proposed and evaluated.

RESULTS

In the simulation analyses based on different types of mathematical models, the JA inferred the model parameters and identified the responsible genomic regions more accurately than the independent analysis (IA). The JA of real plant data provided interesting insights into photosensitivity, which were uncovered by the IA.

AVAILABILITY AND IMPLEMENTATION

The statistical framework is provided by the R package GenomeBasedModel available at https://github.com/Onogi/GenomeBasedModel. All R and C++ scripts used in this study are also available at the site.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic sequencing analysis reveals copy number variations and their associations with economically important traits in beef cattle

Copy number variation (CNV) represents a major source of genetic variation, which may have potentially large effects, including alternating gene regulation and dosage, as well as contributing to gene expression and risk for normal phenotypic variability. We carried out a comprehensive analysis of CNV based on whole genome sequencing in Chinese Simmental beef cattle. Totally, we found 9313 deletion and 234 duplication events, covering 147.5 Mb autosomal regions. Within them, 257 deletion events of high frequency overlapped with 193 known RefGenes. Among these genes, we observed several genes were related to economically important traits, like residual feed intake, immune responding, pregnancy rate and muscle differentiation. Using a locus-based analysis, we identified 11 deletions and 1 duplication, which were significantly associated with three traits including carcass weight, tenderloin and longissimus muscle area. Our sequencing-based study provided important insights into investigating the association of CNVs with important traits in beef cattle.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Genome-Wide Association Study Using Individual Single-Nucleotide Polymorphisms and Haplotypes for Erythrocyte Traits in Alpine Merino Sheep

Adaptation to high-altitude hypoxia is essential for domestic animals, such as yak, Tibetan chicken, and Tibetan sheep, living on high plateaus, as it ensures efficient oxygen absorption and utilization. Red blood cells are the primary medium for transporting oxygen in the blood. However, little is known about the genetic mechanism of erythrocyte traits. Genome-wide association studies (GWASs) based on single markers or haplotypes have identified potential mechanisms for genetic variation and quantitative traits. To identify loci associated with erythrocyte traits, we performed a GWAS based on the method of the single marker and haplotype in 498 Alpine Merino sheep for six erythrocyte traits: red blood cell count (RBC), hemoglobin (HGB), hematocrit (HCT), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), and RBC volume distribution width coefficient of variation (RWD_CV). Forty-two significant single-nucleotide polymorphisms (SNPs) associated with the six erythrocyte traits were detected by means of a single-marker GWAS, and 34 significant haplotypes associated with five erythrocyte traits were detected by means of haplotype analysis. We identified six genes (DHCR24, SPATA9, FLI1, PLCB1, EFNB2, and SH2B3) as potential genes of interest via gene function annotations, location, and expression variation. In particular, FLI1 and PLCB1 were associated with hematopoiesis and erythropoiesis, respectively. These results provide a theoretical basis for analyzing erythrocyte traits and high-altitude hypoxia adaptation in Alpine Merino sheep and will be a useful reference for future studies of plateau-dwelling livestock.

Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

Bovine mastitis is a common inflammatory disease caused by multiple factors in early lactation or dry period. Genome wide association studies (GWAS) can provide a convenient and effective strategy for understanding the biological basis of mastitis and better prevention. 2b-RADseq is a high-throughput sequencing technique that offers a powerful method for genome-wide genetic marker development and genotyping. In this study, single nucleotide polymorphisms (SNPs) of the immune-regulated gene correlative with mastitis were screened and identified by two stage association analysis via GWAS-2b-RADseq in Chinese Holstein cows. We have screened 10,058 high quality SNPs from 7,957,920 tags and calculated their allele frequencies. Twenty-seven significant SNPs were co-labeled in two GWAS analysis models [Bayesian (P < 0.001) and Logistic regression (P < 0.01)], and only three SNPs (rs75762330, C > T, PIC = 0.2999; rs88640083, A > G, PIC = 0.1676; rs20438858, G > A, PIC = 0.3366) were annotated to immune-regulated genes (PTK2B, SYK, and TNFRSF21). Identified three SNPs are located in non-coding regions with low or moderate genetic polymorphisms. However, independent sample population validation (Case-control study) data showed that three important SNPs (rs75762330, P < 0.025, OR > 1; rs88640083, P < 0.005, OR > 1; rs20438858, P < 0.001, OR < 1) were significantly associated with clinical mastitis trait. Importantly, PTK2B and SYK expression was down-regulated in both peripheral blood leukocytes (PBLs) of clinical mastitis cows and in vitro LPS (E. coli)-stimulated bovine mammary epithelial cells, while TNFRSF21 was up-regulated. Under the same conditions, expression of Toll-like receptor 4 (TLR4), AKT1, and pro-inflammatory factors (IL-1β and IL-8) were also up-regulated. Interestingly, network analysis indicated that PTK2B and SYK are co-expressed in innate immune signaling pathway of Chinese Holstein. Taken together, these results provided strong evidence for the study of SNPs in bovine mastitis, and revealed the role of SYK, PTK2B, and TNFRSF21 in bovine mastitis susceptibility/tolerance.

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.

Development of a structural growth curve model that considers the causal effect of initial phenotypes

BACKGROUND

Growth curves have been widely used in genetic analyses to gain insights into the growth characteristics of both animals and plants. However, several questions remain unanswered, including how the initial phenotypes affect growth and what is the duration of any such impact. For beef cattle production in Japan, calves are procured from farms that specialize in reproduction and then moved to other farms where they are fattened to achieve their market/purchase value. However, the causal effect of growth, while calves are on the reproductive farms, on their growth during fattening remains unclear. To investigate this, we developed a model that combines a structural equation with a growth curve model. The causal effect was modeled with B-splines, which allows inference of the effect as a curve. We fitted the proposed structural growth curve model to repeated measures of body weight from a Japanese beef cattle population (n = 3831) to estimate the curve of the causal effect of the calves' initial weight on their trajectory of growth when they are on fattening farms.

RESULTS

Maternal and reproduction farm effects explained 26% of the phenotypic variance of initial weight at fattening farms. The structural growth curve model was fitted to remove the effects of these factors in growth curve analysis at fattening farms. The estimated curve of causal effects remained at approximately 0.8 for 200 d after the calves entered the fattening farms, which means that 64% of the phenotypic variance was explained by the initial weight. Then, the effect decreased linearly and disappeared approximately 620 d after entering the fattening farms, which corresponded to an average age of 871.5 d.

CONCLUSIONS

The proposed model is expected to provide more accurate estimates of genetic values for growth patterns because the confounding causal factors such as maternal and reproduction farm effects are removed. Moreover, examination of the inferred curve of the causal effect enabled us to estimate the effect of a calf's initial weight at arbitrary times during growth, which could provide suitable information for decision-making when shifting the time of slaughter, building models for genetic evaluation, and selecting calves for market.

Use of nonlinear mixed models for describing testicular volume growth curve in Nellore bulls

This study aimed to describe longitudinal testicular volume (TV) data of Nellore bulls by using nonlinear mixed models. Dataset consisted of 2,294 TV measurements from 505 bulls with ages ranging from 563 to 4,307 days. Nine nonlinear models were evaluated: Brody, Gompertz, Hill, Logistic I and II, Meloun, Michaelis-Menten, Mitscherlich and von Bertalanffy. Goodness of fit was evaluated by Akaike's information criterion (AIC), Bayesian information criterion (BIC), adjusted R², percentage of convergence, error sum of squares (ESS), mean absolute deviation (MAD), and average prediction error (APE). These criteria were used to select the best model, then the absolute growth rate (AGR) for TV was estimated by the first derivate of the adjusted model related to time (∂Y/∂t). The values of adjusted R², ESS, MAD and APE were similar among models. Percentage of convergence was higher for the Logistic I (76.8%), Logistic II (75.5%) and Mitscherlich (78.6%) models, but Logistic I and II showed the lowest values of AIC and BIC, indicating a better fit, so the Logistic I model was chosen for subsequent analyses. The TV growth occurred at a high rate until the inflection point, which was estimated at approximately 22 months of age; it stabilized and reached a plateau at approximately 2,500 days of age. This may suggest that TV is more related to sexual maturity than precocity. Additionally, the Logistic I model was used to estimate the growth curve and the AGR of the testicular length and width. As a result, testicular length increased at a higher rate than testicular width until approximately 1,600 days of age, indicating that the testes become longer with increasing age.

Evaluation of feed efficiency traits for genetic improvement in Japanese Black cattle

We evaluated the genetic relationships (1) among feed efficiency traits with different fattening periods, (2) between feed efficiency traits and growth traits, and (3) between feed efficiency traits and carcass traits, to determine the influence of genetic factors on feed efficiency traits. In total, 4,578 Japanese Black cattle from a progeny testing program were used. Residual feed intake (RFI), residual BW gain (RG), and residual intake and BW gain (RIG) were defined as feed efficiency traits, and were measured for the first half (approximately 9 to 15 months of age), latter half (approximately 15 to 21 months of age), and total period of fattening (approximately 9 to 21 months of age). A single-trait animal model for estimating heritability and a two-trait animal model for estimating genetic and phenotypic correlations were used. The heritability estimates for RFI, RG, and RIG were different in each fattening period, ranging from 0.36 to 0.46, 0.19 to 0.28, and 0.28 to 0.34, respectively, and the heritability estimates for the total fattening period were greater than those for the first and latter halves separately. RIG showed the greatest preferred genetic correlation, with a greater feed conversion ratio than the other feed efficiency traits (ranging from -0.84 to -0.96). RG in the first and latter halves of the fattening period had different genetic correlations with the growth starting point (0.82 and -0.06, respectively) and maturity rate (0.49 and -0.51, respectively) of the Gompertz growth curve parameters, and is strongly dependent on the different fattening periods. Feed efficiency traits in different fattening periods had low genetic correlations with the carcass traits (from -0.05 to 0.19 for RFI; from 0.02 to 0.31 for RG; and from -0.11 to 0.20 for RIG). This study indicated the possibility for genetic improvement through the selection of high-RIG animals to decrease feed intake and increase BW gain without any unfavorable correlated responses affecting mature (asymptotic) weight and carcass grade.

Multiple-trait genomewide mapping and gene network analysis for scrotal circumference growth curves in Brahman cattle

Fertility traits are economically important in cattle breeding programs. Scrotal circumference (SC) measures are repeatable, easily obtained, highly heritable, and positively correlated with female fertility traits and sperm quality traits in males. A useful approach to summarize SC measures over time is using nonlinear models, which summarize specific measures of SC in a few parameters with biological interpretation. This approach facilitates the selection of bulls with larger SC and maturity index (K), that is, early maturing animals. Because SC is a sex-limited trait, identifying the underlying genomics of growth curve parameters will allow selection across both males and females. We reported the first multitrait genomewide association study (GWAS) of estimated growth curve parameters for SC data in Brahman cattle. Five widely used nonlinear models were tested to fit a total of 3,612 SC records, measured at 6, 12, 18, and 24 mo of age. The von Bertalanffy model, individually fitted for each animal, best fit this SC data. Parameter estimates SC at maturity (A) and K as well as SC at all ages were jointly analyzed in a GWAS to identify 1-Mb regions most strongly associated with each trait. Heritabilities were 0.25 for K and 0.32 for A and ranged from 0.51 to 0.72 for SC at 6 (SC6), 12 (SC12), 18 (SC18), and 24 mo of age (SC24). An overlapping window on chromosome 14 explaining around 0.8% of genetic variance for K, SC12, SC18, and SC24 was observed. The major positional candidate genes within 1 Mb upstream and downstream of this overlapping window were , , , and . Windows of 1 Mb explaining more than 0.4% of each trait on chromosomes 1, 3, 6, 7, 14, 17, 18, 24, 25, and 26 were identified. Pathways and net-work analyses were indicated through transcription factors playing a role on fertility traits: , , , , , , and . Further validation studies on larger populations or other breeds are required to validate these findings and to improve our understanding of the biology and complex genetic architecture of traits associated with scrotal growth and male fertility in cattle.

Genome-wide association study for birth, weaning and yearling weight in Colombian Brahman cattle

Genotypic and phenotypic data of 1,562 animals were analyzed to find genomic regions that potentially influence the birth weight (BW), weaning weight at seven months of age (WW) and yearling weight (YW) of Colombian Brahman cattle, with genotyping conducted using Illumina Bead chip array with 74,669 SNPs. A Single Step Genomic BLUP (ssGBLP), approach was used to estimate the proportion of variance explained by each marker. Multiple regions scattered across the genome were found to influence weights at different ages, also dependent on the trait component (direct or maternal). The most interesting regions were connected to previously identified QTLs and genes, such as ADAMTSL3, CAPN2, CAPN2, FABP6, ZEB2 influencing growth and weight traits. The identified regions will contribute to the development and refinement of genomic selection programs for Zebu Brahman cattle in Colombia.

Genome-Wide Association Study Reveals the PLAG1 Gene for Knuckle, Biceps and Shank Weight in Simmental Beef Cattle

Carcass traits of beef cattle have been genetically improved to increase yield of high quality meat. Genome-wide association study (GWAS) is a powerful method to identify genetic variants associated with carcass traits. For the 770K genotyped SNPs from 1141 Chinese Simmental cattle, we used the compressed mixed linear model (CMLM) to perform a genome-wide association study for knuckle, biceps and shank of beef carcass traits. Seventeen significantly associated SNPs were found, which are located on BTA6, BTA14 and BTA15. Interestingly, one pleiotropic quantitative trait nucleotide (QTN), named BovineHD1400007259 (p < 10⁻⁸) within the well-known gene region PLAG1-CHCHD7 on BTA14, was found to govern variation of the knuckle, biceps and shank traits. The QTN accounted for 8.6% of phenotypic variance for biceps. In addition, 16 more SNPs distributed on BTA14 were detected as being associated with the carcass traits.

Rab28 is a TBC1D1/TBC1D4 substrate involved in GLUT4 trafficking

The Rab-GTPase-activating proteins (GAPs) TBC1D1 and TBC1D4 play important roles in the insulin-stimulated translocation of the glucose transporter GLUT4 from intracellular vesicles to the plasma membrane in muscle cells and adipocytes. We identified Rab28 as a substrate for the GAP domains of both TBC1D1 and TBC1D4 in vitro. Rab28 is expressed in adipose cells and skeletal muscle, and its GTP-binding state is acutely regulated by insulin. We found that in intact isolated mouse skeletal muscle, siRNA-mediated knockdown of Rab28 decreases basal glucose uptake. Conversely, in primary rat adipose cells, overexpression of Rab28-Q72L, a constitutively active mutant, increases basal cell surface levels of an epitope-tagged HA-GLUT4. Our results indicate that Rab28 is a novel GTPase involved in the intracellular retention of GLUT4 in insulin target cells.

Comparison on accuracy of different nonlinear models in predicting growth of Podolica bulls

Animal growth does not follow a linear pattern, being explained mathematically by functions that have parameters with biological meaning. These parameters are used to estimate the expected weight of animals at specific ages. Several nonlinear models have been used to describe growth. This study was carried out to estimate the parameters of logistic, Gompertz, Richards and von Bertalanffy growth curve models in a sample of Podolica young bulls to determine the goodness of fit. Animals were weighed every 3 months from birth to 810 days of age. The results indicate that all the growth models used were easily fitted to the observed data with Gompertz and logistic functions presenting less computational difficulty in terms of number of iterations to achieve convergence. Moreover, logistic and Richards equations provided the best overall fit being useful to describe the growth of Podolica bulls. Considering that the literature lacks information on growth curves in Podolica breed, the study of a mathematical model for growth describing the developmental pattern of a specific population within a peculiar environment is a useful tool to improve Podolica breed production.

Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH ) with pedigree-based inbreeding coefficients (F_PED ). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.