A genome-wide association study of direct gestation length in US Holstein and Italian Brown populations

Direct gestation length influences economically important traits in dairy cattle that are related to birth and peri-natal survival of the calf. The objective of this study was to identify single nucleotide polymorphisms (SNPs) that are significantly associated with direct gestation length through a genome-wide association study. Data used in the analysis included 7,308,194 cow gestation lengths from daughters of 4743 United States Holstein sires in the Cooperative Dairy DNA Repository population and 580,157 gestation lengths from 749 sires in the Italian Brown population. Association analysis included 36,768 and 35,082 SNPs spanning all autosomes for Holstein and Brown Swiss, respectively. Multiple shrinkage Bayesian was employed. Estimates of heritability for both populations were moderate, with values of 0.32 (±0.03) and 0.29 (±0.02) for Holstein and Brown Swiss, respectively. A panel of SNPs was identified, which included SNPs that have significant effects on direct gestation length, of which the strongest candidate region is located on chromosome 18. Two regions not previously linked to direct calving ease and calf survival were identified on chromosome 7 and 28, corresponding to regions that contain genes related to embryonic development and foetal development. SNPs were also identified in regions that have been previously mapped for calving difficulty and longevity. This study identifies target regions for the investigation of direct foetal effects, which are a significant factor in determining the ease of calving.     

Body traits,carcass characteristics and price of cull cows as affected by farm type,breed, age and calving to culling interval

Beef production from cull cows is an additional source of income for dairy farms and greatly contributes to red meat production, but the sources of variation of live animal characteristics and the carcass traits of cull cows have rarely been examined. This study investigated the effects of the farm type, breed, age at slaughter (AGE) and calving to culling interval (Calv_Cull) on the body traits and carcass characteristics of dairy and dual-purpose cull cows. Data from 555 cull cows from 182 herds belonging to five farm types, characterised by a combination of housing and feeding systems, were recorded and analysed. Dairy breeds, such as Holstein Friesian and Brown Swiss, and dual-purpose breeds (Simmental, Rendena) were included in the trait assessments. The day before slaughter, the cows were weighed and scored for body condition (BCS) and fleshiness, and then, their heart girth and wither height were measured. At the slaughterhouse, the carcass weight (CW), dressing percentage (DP), carcass conformation and fatness scores, carcass price per kg and carcass total value were obtained. On average, the cows were slaughtered at nearly 71±27 months of age, 285±187 days after the last calving; 615±95 kg BW; and provided a 257±51 kg CW. Nearly 50% of the cows fell within the BCS range of 2.75 to 3.50, and the carcasses were mostly graded in the lowest class of conformation and fatness scores. Cull cows from free-stall farms had a higher DP, carcass conformation score and price than those from traditional tie-stall farms. The breed influenced the AGE, live animal characteristics and carcass traits. Cows from dairy breeds were younger at slaughter, had a lower BCS and fleshiness, and greater body measurements, but a lower DP and carcass price than those from dual-purpose breeds, although differences between the breeds were found within both groups. The age of the cows at slaughter influenced the Calv_Cull and increased the BW, body measurements and CW, but not the fleshiness and fatness appreciation (both in vivo and postmortem) or carcass price. The increasing Calv_Cull improved the BW, BCS, fleshiness, CW and carcass conformation and fatness. In conclusion, the decision to cull dairy cows should also take into account the factors that affect their carcass value in regards to improving the carcass price of cows.     

Unifying Genetic Canalization,Genetic Constraint,and Genotype-by-Environment Interaction: QTL by Genomic Background by Environment Interaction of Flowering Time in Boechera stricta

Natural populations exhibit substantial variation in quantitative traits. A quantitative trait is typically defined by its mean and variance, and to date most genetic mapping studies focus on loci altering trait means but not (co)variances. For single traits, the control of trait variance across genetic backgrounds is referred to as genetic canalization. With multiple traits, the genetic covariance among different traits in the same environment indicates the magnitude of potential genetic constraint, while genotype-by-environment interaction (GxE) concerns the same trait across different environments. While some have suggested that these three attributes of quantitative traits are different views of similar concepts, it is not yet clear, however, whether they have the same underlying genetic mechanism. Here, we detect quantitative trait loci (QTL) influencing the (co)variance of phenological traits in six distinct environments in Boechera stricta, a close relative of Arabidopsis. We identified nFT as the QTL altering the magnitude of phenological trait canalization, genetic constraint, and GxE. Both the magnitude and direction of nFT''s canalization effects depend on the environment, and to our knowledge, this reversibility of canalization across environments has not been reported previously. nFT''s effects on trait covariance structure (genetic constraint and GxE) likely result from the variable and reversible canalization effects across different traits and environments, which can be explained by the interaction among nFT, genomic backgrounds, and environmental stimuli. This view is supported by experiments demonstrating significant nFT by genomic background epistatic interactions affecting phenological traits and expression of the candidate gene, FT. In contrast to the well-known canalization gene Hsp90, the case of nFT may exemplify an alternative mechanism: Our results suggest that (at least in traits with major signal integrators such as flowering time) genetic canalization, genetic constraint, and GxE may have related genetic mechanisms resulting from interactions among major QTL, genomic backgrounds, and environments.     

Genetic parameters for feed efficiency in Romane rams and responses to single-generation selection

Feeding costs represent one of the highest expenditures in animal production systems. Breeding efficient animals that express their growth potential while eating less is therefore a major objective for breeders. We estimated the genetic parameters for feed intake, feed efficiency traits (residual feed intake (RFI) and feed conversion ratio (FCR)), growth and body composition traits in the Romane meat sheep breed. In these traits, selection responses to single-generation divergent selection on RFI were evaluated. From 2009 to 2016, a total of 951 male lambs were tested for 8 weeks starting from 3 months of age. They were weighed at the beginning and at the end of the testing period. Backfat thickness and muscle depth were recorded at the end of the testing period through ultrasound measurements. Feed intake was continuously recorded over the testing period using the automatic concentrate feeders. The heritability of RFI was estimated at 0.45 ± 0.08, which was higher than the heritability of FCR (0.30 ± 0.08). No significant genetic correlations were observed between RFI and growth traits. A favourable low negative genetic correlation was estimated between RFI and muscle depth (−0.30 ± 0.15), though additional data are needed to confirm these results. The selection of low RFI sires based on their breeding values led to the production of lambs eating significantly less concentrate (3% decrease in the average daily feed intake), but with the same growth as lambs from sires selected based on high RFI breeding values. We concluded that in meat sheep, RFI is a heritable trait that is genetically independent of post-weaning growth and body composition traits. A one-generation divergent selection based on RFI breeding values highlighted that substantial gains in feeding costs can be expected in selection schemes for meat sheep breeds.     

Genome-Wide Contribution of Genotype by Environment Interaction to Variation of Diabetes-Related Traits

While genome-wide association studies (GWAS) and candidate gene approaches have identified many genetic variants that contribute to disease risk as main effects, the impact of genotype by environment (GxE) interactions remains rather under-surveyed. To explore the importance of GxE interactions for diabetes-related traits, a tool for Genome-wide Complex Trait Analysis (GCTA) was used to examine GxE variance contribution of 15 macronutrients and lifestyle to the total phenotypic variance of diabetes-related traits at the genome-wide level in a European American population. GCTA identified two key environmental factors making significant contributions to the GxE variance for diabetes-related traits: carbohydrate for fasting insulin (25.1% of total variance, P-nominal = 0.032) and homeostasis model assessment of insulin resistance (HOMA-IR) (24.2% of total variance, P-nominal = 0.035), n-6 polyunsaturated fatty acid (PUFA) for HOMA-β-cell-function (39.0% of total variance, P-nominal = 0.005). To demonstrate and support the results from GCTA, a GxE GWAS was conducted with each of the significant dietary factors and a control E factor (dietary protein), which contributed a non-significant GxE variance. We observed that GxE GWAS for the environmental factor contributing a significant GxE variance yielded more significant SNPs than the control factor. For each trait, we selected all significant SNPs produced from GxE GWAS, and conducted anew the GCTA to estimate the variance they contributed. We noted the variance contributed by these SNPs is higher than that of the control. In conclusion, we utilized a novel method that demonstrates the importance of genome-wide GxE interactions in explaining the variance of diabetes-related traits.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

Effectiveness of genomic prediction on milk flow traits in dairy cattle

Background

Milkability, primarily evaluated by measurements of milking speed and time, has an economic impact in milk production of dairy cattle. Recently the Italian Brown Swiss Breeders Association has included milking speed in genetic evaluations. The main objective of this study was to investigate the possibility of implementing genomic selection for milk flow traits in the Italian Brown Swiss population and thereby evaluate the potential of genomic selection for novel traits in medium-sized populations. Predicted breeding values and reliabilities based on genomic information were compared with those obtained from traditional breeding values.

Methods

Milk flow measures for total milking time, ascending time, time of plateau, descending time, average milk flow and maximum milk flow were collected on 37 213 Italian Brown Swiss cows. Breeding values for genotyped sires (n = 1351) were obtained from standard BLUP and genome-enhanced breeding value techniques utilizing two-stage and single-step methods. Reliabilities from a validation dataset were estimated and used to compare accuracies obtained from parental averages with genome-enhanced predictions.

Results

Genome-enhanced breeding values evaluated using two-stage methods had similar reliabilities with values ranging from 0.34 to 0.49 for the different traits. Across two-stage methods, the average increase in reliability from parental average was approximately 0.17 for all traits, with the exception of descending time, for which reliability increased to 0.11. Combining genomic and pedigree information in a single-step produced the largest increases in reliability over parent averages: 0.20, 0.24, 0.21, 0.14, 0.20 and 0.21 for total milking time, ascending time, time of plateau, descending time, average milk flow and maximum milk flow, respectively.

Conclusions

Using genomic models increased the accuracy of prediction compared to traditional BLUP methods. Our results show that, among the methods used to predict genome-enhanced breeding values, the single-step method was the most successful at increasing the reliability for most traits. The single-step method takes advantage of all the data available, including phenotypes from non-genotyped animals, and can easily be incorporated into current breeding evaluations.     

Portion control for the treatment of obesity in the primary care setting

Background

Marbling defined by the amount and distribution of intramuscular fat, so-called Shimofuri, is an economically important trait of beef cattle in Japan. Our previous study detected 3 single nucleotide polymorphisms (SNPs), g.231054C > T, g.3109537C > T and c.*188G > A, respectively, in the 5' flanking region of the titin (TTN), the 5' flanking region of the ribosomal protein L27a (RPL27A) and the 3' untranslated region of the akirin 2 genes (AKIRIN2), which have been considered as positional functional candidates for the genes responsible for marbling, and showed association of these SNPs with marbling in Japanese Black beef cattle. In the present study, we investigated the allele frequency distribution of the 3 SNPs among the 5 cattle breeds, Japanese Black, Japanese Brown, Japanese Shorthorn, Holstein and Brown Swiss breeds.

Findings

We genotyped the TTN g.231054C > T, RPL27A g.3109537C > T and AKIRIN2 c.*188G > A SNPs by polymerase chain reaction-restriction fragment length polymorphism method, using 101 sires and 1,705 paternal half sib progeny steers from 8 sires for Japanese Black, 86 sires and 27 paternal half sib progeny steers from 3 sires for Japanese Brown, 79 sires and 264 paternal half sib progeny steers from 14 sires for Japanese Shorthorn, 119 unrelated cows for Holstein, and 118 unrelated cows for Brown Swiss breeds. As compared to the frequencies of the g.231054C > T T, g.3109537C > T T and c.*188G > A A alleles, associated with high marbling, in Japanese Black breed that has been subjected to a strong selection for high marbling, those in the breeds, Japanese Shorthorn, Holstein and Brown Swiss breeds, that have not been selected for high marbling were null or lower. The Japanese Brown breed selected slightly for high marbling showed lower frequency than Japanese Black breed in the g.3109537C > T T allele, whereas no differences were detected between the 2 breeds in the frequencies of the g.231054C > T T and c.*188G > A A alleles.

Conclusions

Based on this finding, we hypothesized that the pressure of the strong selection for high marbling in Japanese Black breed has increased the frequencies of the T, T and A alleles at the TTN g.231054C > T, RPL27A g.3109537C > T and AKIRIN2 c.*188G > A SNPs, respectively. This study, together with the previous association studies, suggested that the 3 SNPs may be useful for effective marker-assisted selection to increase the levels of marbling.     

Genome-wide interaction of genotype by erythrocyte n-3 fatty acids contributes to phenotypic variance of diabetes-related traits

Background

Little is known about the interplay between n-3 fatty acids and genetic variants for diabetes-related traits at the genome-wide level. The present study aimed to examine variance contributions of genotype by environment (GxE) interactions for different erythrocyte n-3 fatty acids and genetic variants for diabetes-related traits at the genome-wide level in a non-Hispanic white population living in the U.S.A. (n = 820). A tool for Genome-wide Complex Trait Analysis (GCTA) was used to estimate the genome-wide GxE variance contribution of four diabetes-related traits: HOMA-Insulin Resistance (HOMA-IR), fasting plasma insulin, glucose and adiponectin. A GxE genome-wide association study (GWAS) was conducted to further elucidate the GCTA results. Replication was conducted in the participants of the Boston Puerto Rican Health Study (BPRHS) without diabetes (n = 716).

Results

In GOLDN, docosapentaenoic acid (DPA) contributed the most significant GxE variance to the total phenotypic variance of both HOMA-IR (26.5%, P-nominal = 0.034) and fasting insulin (24.3%, P-nominal = 0.042). The ratio of arachidonic acid to eicosapentaenoic acid + docosahexaenoic acid contributed the most significant GxE variance to the total variance of fasting glucose (27.0%, P-nominal = 0.023). GxE variance of the arachidonic acid/eicosapentaenoic acid ratio showed a marginally significant contribution to the adiponectin variance (16.0%, P-nominal = 0.058). None of the GCTA results were significant after Bonferroni correction (P < 0.001). For each trait, the GxE GWAS identified a far larger number of significant single-nucleotide polymorphisms (P-interaction ≤ 10E-5) for the significant E factor (significant GxE variance contributor) than a control E factor (non-significant GxE variance contributor). In the BPRHS, DPA contributed a marginally significant GxE variance to the phenotypic variance of HOMA-IR (12.9%, P-nominal = 0.068) and fasting insulin (18.0%, P-nominal = 0.033).

Conclusion

Erythrocyte n-3 fatty acids contributed a significant GxE variance to diabetes-related traits at the genome-wide level.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-781) contains supplementary material, which is available to authorized users.     

Indirect genetic effects on the relationships between production and feeding behaviour traits in growing Duroc pigs

Performance and feeding behaviour traits in growing pigs could be affected by social interaction effects when animals are raised in group. So, properly knowing the genetic correlations between direct and social interaction effects among performance and feeding behaviour traits could improve the accuracy of the genetic evaluations. Our aim was to explore the role of feeding behaviour traits (FBT) and indirect genetic effects (IGEs) in the genetic evaluations of growing pigs. Thus, genetic parameters were estimated for production traits (PT): average daily gain, average daily feed consumption, feed conversion ratio and backfat thickness; as well as for FBT: average daily feeding rate, average daily feeding frequency, average daily occupation time and average daily time between consecutive visits. Traits were recorded in 1144 Duroc pigs during the fattening period. Two bivariate models were fitted: classic animal model and an animal model fitting IGE. Estimations were done following Bayesian procedures. Heritability estimates obtained with classic animal model for all studied traits were medium-high. The additional heritable variation captured by IGE supposed that the ratios of total genetic variance to phenotypic variance (T²) were higher than the heritability estimates obtained with the classic model, except for occupation time trait, when a lower value (0.20 ± 0.19) was estimated. This is due to a high and negative correlation between IGE and direct genetic effects (DGEs) of this particular trait (−0.78 ± 0.27). Results from classic animal model do not evidence a clear role of FBT to improve the accuracy of breeding value predictions for PT; only average daily feeding rate seems to show a positive correlation (around 0.50 to 0.60) with average daily gain, average daily feed consumption and backfat thickness. However, when IGE model was fitted, the number of estimates of genetic correlations between FBT and PT showing a relevant magnitude increased, generally for the correlations between IGE of FBT and DGE of PT; or particularly for the correlations between IGE of average daily feeding frequency, and the IGE of all the PT, except average daily gain. Thus, in evaluations using the animal model with IGE fitted, the inclusion of FBT could aid the improvement of the accuracy of breeding value predictions for PT. This is a consequence of the improved genetic relationships between traits that can be fitted when considering such models.     

Heritability estimates and adjustment factors for the effects of bull age and age of dam on yearling testicular size in breeds of bulls

Scrotal circumference, testicular length and body weight were measured in 3,090 yearling bulls of 12 breed groups finishing growth performance tests during a 5 yr period. Breeds were Limousin, Hereford, Charolais, Angus, Red Poll, Simmental, Pinzgauer, Brown Swiss, Gelbvieh and three crossbred breed groups. All bulls were born during a 60-d calving season starting in late March of each year and were subjected to similar management and environmental influences during the study. The bulls were the progeny of 307 sires averaging approximately 26 sires per breed group and 10 sons per sire. Breed group, sire within breed group, year and age-of-dam effects were important (P<0.01) for all testicular traits at both a constant age (354 d) and constant body weight (418 kg). Paternal half-sib estimates of heritability were 0.41 +/- 0.06 for ageconstant scrotal circumference, 0.34 +/- 0.06 for testicular length and 0.37 +/- 0.06 for calculated paired testicular volume, respectively. Age-constant genetic correlations between yearling body weight and testicular traits were small, indicating that testicular growth and body growth rates are largely independent, regardless of breed. Testicular size of bulls from 2-yr-old dams was smaller than that of bulls from older dams. Most of these age-of-dam effects on testicular size were removed when testicular size was adjusted for the effects of body weight, suggesting that age-of-dam effects on testicular size are primarily the result of age-of-dam effects on body weight. Age adjustment factors for yearling scrotal circumference did not differ (P>0.20) among breed groups and averaged 0.032 cm per day of age. Adjustment factors for age of dam were +1.3, +0.8, +0.4, and +0.0 cm for sons of 2-, 3-, 4- and 5-yr-old dams, respectively.     

Genetic parameters for production,health, fertility and longevity traits in dairy cows

Milk production, fertility, longevity and health records, were extracted from databases of two milk recording organisations in the United Kingdom for the first three lactations of the Holstein–Friesian breed. These included data related to health events (mastitis and lameness), voluntarily recorded on a proportion of farms. The data were analysed to calculate disease incidence levels and to estimate genetic parameters for health traits and their relationships with production and other functional traits. The resulting dataset consisted of 124 793 lactations from 75 137 animals of 1586 sires, recorded in 2434 herds. Incidence of health events increased with parity. The overall incidence of mastitis (MAS) and lameness (LAM), defined as binary traits, were 17% and 16%, respectively. Heritability estimates for MAS and LAM were 0.04 and 0.02, respectively, obtained from repeatability linear sire models. Heritability estimates of mastitis and lameness as count traits were slightly higher, 0.05 and 0.03, respectively. Genetic correlations were obtained by bivariate analyses of all pair-wise combinations between milk 305-day yield (MY), protein 305-day yield (PY), fat 305-day yield (FY), lactation average log_e transformed lactation average somatic cell count (SCS), calving interval (CI), days to first service (DFS), non-return at 56 days (NR56), number of inseminations (NINS), mastitis (MAS), number of mastitis episodes (NMAS), lameness (LAM), number of lameness episodes (NLAM) and lifespan score (LS). As expected, MAS was correlated most strongly with SCS (0.69), which supports the use of SCS as an indicator trait for mastitis. Genetic correlations between MAS and yield and fertility traits were of similar magnitude ranging from 0.27 to 0.33. Genetic correlations between MAS with LAM and LS were 0.38 and −0.59, respectively. Not all genetic correlations between LAM and other traits were significant because of fewer numbers of lameness records. LAM had significant genetic correlations with MY (0.38), PY (0.28), CI (0.35), NINS (0.38) and LS (−0.53). The heritability estimates of mastitis and lameness were low; therefore, genetic gain through direct selection alone would be slow, yet still positive and cumulative. Direct selection against mastitis and lameness as additional traits should reduce incidence of both diseases, and simultaneously improve fertility and longevity. However, both health traits had antagonistic relationships with production traits, thus genetic gain in production would be slower.     

Genome-wide association study of economically important traits in Charolais and Limousin beef cows

Genomic selection has proven effective for advancing genetic gain for key profit traits in dairy cattle production systems. However, its impact to-date on genetic improvement programs for beef cattle has been less effective. Despite this, the technology is thought to be particularly useful for low heritability traits such as those associated with reproductive efficiency. The objective of this study was to identify genetic variants associated with key determinants of reproductive and overall productive efficiency in beef cows. The analysis employed a large dataset derived from the national genetic evaluation program in Ireland for two of the most predominant beef breeds, viz. Charolais (n = 5 244 cows) and Limousin (n = 7 304 cows). Single nucleotide polymorphisms (SNPs) were identified as being statistically significantly associated (adj. P < 0.05) with both reproductive and productive traits for both breed types. However, there was little across breed commonality, with only two SNPs (rs110240246 and rs110344317; adj. P < 0.05) located within the genomic regions of the LCORL and MSTN genes respectively, identified in both Charolais and Limousin populations, associated with traits including carcass weight, cull-cow weight and live-weight. Significant SNPs within the MSTN gene were also associated with both reproduction and production related traits within each breed. Finally, traits including calving difficulty, calf mortality and calving interval were associated with SNPs within genomic regions comprising genes involved in cellular growth and lipid metabolism. Genetic variants identified as associated with both important reproductive efficiency and production related traits from this study warrant further analyses for their potential incorporation into breeding programmes to support the sustainability of beef cattle production.     

Selection for litter size and litter birthweight in Large White pigs: Maximum,mean and variability of reproduction traits

Gradually increasing trend of litter size poses a challenge to pig farmers in terms of managing larger litters. Therefore, it seems that a balanced approach that optimises litter size, litter birthweight, and uniformity of those traits is needed in order to address animal welfare and farm management concerns. This study aimed to investigate this issue by defining several traits for total number born (TNB), number born alive (NBA) and litter birthweight (LW). First, the highest value from at least five records per sow was selected as maximum (max) value for each reproduction trait. Second, a mean (mean) for each reproduction trait was calculated per sow. Last, the variability of reproduction traits between parities of the sow was calculated as log-transformed variance of residuals of all observations per sow for each reproduction trait (LnVar). In total, 23 193 Large White sows from Topigs Norsvin with 152 282 litter records were used for analysis in ASReml 4.1. Also, a simulation of breeding schemes was performed with the use of SelAction 2.1 and estimates from genetic analysis. Maximum value of reproductive traits had a much higher heritability than repeated observations or mean of reproduction traits, e.g., 0.31 for maxTNB vs. 0.12 for TNB and 0.07 for meanTNB, which allows for a faster response under selection. The maximum value traits, however, were found to carry more risks, i.e. higher ratio of stillborn (not for maxNBA) and increased variability of traits. Thus, using them in breeding programme should be carefully considered. The genetic coefficient of variation on SD level estimated to indicate the genetic magnitude for variability phenotypes indicated a maximum change of 6–9% in genetic SD of TNB, NBA and LW. The genetic correlations between mean and corresponding variability traits varied from 0.66 to 0.74, whereas the correlation between other mean and variability traits ranged from 0.33 to 0.99. The simulation indicated that even with selection targeted against the variability of reproduction traits, a very limited change should be expected due to a complex genetic and phenotypic relationship between the traits. In the scenarios with selection against LnVarTNB and LnVarLW, this was a decrease of 0.1–0.6% per year, whereas in scenario with selection against LnVarNBA, the range was 0.6–1.1% per year. It is still possible to increase litter size and birthweight further, however, a balance between mean and variability of reproduction traits is required, which can be obtained only by a very well designed breeding programme.     

Quantitative trait loci mapping for conjugated linoleic acid,vaccenic acid and ∆9‐desaturase in Italian Brown Swiss dairy cattle using selective DNA pooling

A selective DNA pooling approach was applied to identify QTL for conjugated linoleic acid (CLA), vaccenic acid (VA) and Δ⁹‐desaturase (D9D) milk content in Italian Brown Swiss dairy cattle. Milk samples from 60 animals with higher values (after correction for environmental factors) and 60 animals with lower values for each of these traits from each of five half‐sib families were pooled separately. The pools were genotyped using the Illumina BovineSNP50 BeadChip. Sire allele frequencies were compared between high and low tails at the sire and marker level for SNPs for which the sires were heterozygous. An r procedure was implemented to perform data analysis in a selective DNA pooling design. A correction for multiple tests was applied using the proportion of false positives among all test results. BTA 19 showed the largest number of markers in association with CLA. Associations between SNPs and the VA and Δ⁹‐desaturase traits were found on several chromosomes. A bioinformatics survey identified genes with an important role in pathways for milk fat and fatty acids metabolism within 1 Mb of SNP markers associated with fatty acids contents.     

Associations between the breeding values of Holstein-Friesian bulls and longevity and culling reasons of their daughters

Taking into account functional traits in the breeding practice should lead to a longer productive life of cows. However, despite the increased contribution of these traits in bull selection indices, their daughters are frequently culled as early as the 2nd or 3rd lactation. The problem is whether and to what extent the genetic potential of animals is realized in the production practice. Therefore, the purpose of this study was to determine the associations between the breeding value (BV) of bulls and their daughters for cow longevity and culling reasons in the Holstein-Friesian cattle population in Poland. Data for 532 062 cows culled in 2012, 2015, and 2018 were analyzed. A majority of 5 045 cow sires originated from Poland, Germany, France, the Netherlands, and the United States. The highest variation in the contribution of culling reasons was for the cows culled at the age of 2–4 years. The contribution of the culling reasons, analyzed in relation to the cow culling age, remained similar and the only exception was culling because of old age, for which a significant increase was observed only for the culling age of at least 9 years (13.8%), which was reached by only 7.3% of the cows. The sires were characterized by generally high BV for conformation and reproductive traits. However, they had, at most, the average genetic potential for functional longevity. There were a number of beneficial associations found between the BV of bulls and the distribution of culling reasons in their daughters. For example, it concerns relations between the somatic cell score in milk and culling due to udder diseases and low milk yield, between the interval from calving to first insemination and low milk yield, between the protein yield and old age, or between the BV for certain conformation traits (size, udder) and cow culling due to age. In these cases, as the BV increased for a given trait, the contribution of the corresponding cow culling reason tended to decrease. Our study showed that it seems reasonable to consider Holstein-Friesian cows aged at least 9 years at culling to be long-living animals. This is primarily evidenced by the rapid increase in the culling due to old age in relation to younger cows. Nowadays the above age limit can be suggested as a criterion of longevity for Holstein-Friesian cows but the criterion should be updated to the relation genotype-environment-economy that tends to change over time.     

Genetic architecture controlling variation in grain carotenoid composition and concentrations in two maize populations

Key message

Genetic control of maize grain carotenoid profiles is coordinated through several loci distributed throughout three secondary metabolic pathways, most of which exhibit additive, and more importantly, pleiotropic effects.

Abstract

The genetic basis for the variation in maize grain carotenoid concentrations was investigated in two F_2:3 populations, DEexp × CI7 and A619 × SC55, derived from high total carotenoid and high β-carotene inbred lines. A comparison of grain carotenoid concentrations from population DEexp × CI7 grown in different environments revealed significantly higher concentrations and greater trait variation in samples harvested from a subtropical environment relative to those from a temperate environment. Genotype by environment interactions was significant for most carotenoid traits. Using phenotypic data in additive, environment-specific genetic models, quantitative trait loci (QTL) were identified for absolute and derived carotenoid traits in each population, including those specific to the isomerization of β-carotene. A multivariate approach for these correlated traits was taken, using carotenoid trait principal components (PCs) that jointly accounted for 97 % or more of trait variation. Component loadings for carotenoid PCs were interpreted in the context of known substrate-product relationships within the carotenoid pathway. Importantly, QTL for univariate and multivariate traits were found to cluster in close proximity to map locations of loci involved in methyl-erythritol, isoprenoid and carotenoid metabolism. Several of these genes, including lycopene epsilon cyclase, carotenoid cleavage dioxygenase1 and beta-carotene hydroxylase, were mapped in the segregating populations. These loci exhibited pleiotropic effects on α-branch carotenoids, total carotenoid profile and β-branch carotenoids, respectively. Our results confirm that several QTL are involved in the modification of carotenoid profiles, and suggest genetic targets that could be used for the improvement of total carotenoid and β-carotene in future breeding populations.     

Unraveling the actual background of second litter syndrome in pigs: based on Large White data

Second litter syndrome (SLS) in sows is when fertility performance is lower in the second parity than in the first parity. The causes of SLS have been associated with lactation weight loss, premature first insemination, short lactation length, short weaning to insemination interval, season, and farm of farrowing. There is little known about the genetic background of SLS or if it is a real biological problem or just a statistical issue. Thus, we aimed to evaluate risk factors, investigate genetic background of SLS, and estimate the probability of SLS existing due to the statistical properties of the trait. The records of 246 799 litters (total number born, TNB) from 46 218 Large White sows were used. A total of 15 398 sows had SLS. Two traits were defined: first a binominal trait if a sow had SLS or not (biSLS) and second a continuous trait (Range) created by subtracting the total number of piglets born in the first parity (TNB1) from the piglets born in the second parity (TNB2). Lactation length, farm, and season of the farrowing had significant effects on SLS traits when tested as fixed effects in the genetic model. These effects are farm management-related factors. The age at first insemination and weaning to insemination interval were significant only for other reproduction traits (e.g., TNB1, TNB2, litter weight in parity 1 and 2). The heritability of biSLS was 0.05 (on observed scale), whereas heritability of Range was 0.03. To verify the existence of SLS data with records of 50 000 sows and 9 parities was simulated. The simulations showed that the average expected frequency of SLS across all the parities was 0.49 (±0.05) while the observed frequency in the actual data was 0.46 (±0.04). We compared this to SLS frequencies in 67 farms and only 2 farms had more piglets born in the first parity compared to the second. Therefore, on the individual sow level SLS is likely due to statistical properties of the trait, whereas on the farm level SLS is likely due to farm management. Thus, SLS should not be considered an abnormality nor a syndrome if on average the herd litter size in parity 2 is larger than in parity 1.     

Plasticity in Vegetative Growth over Contrasted Growing Sites of an F1 Olive Tree Progeny during Its Juvenile Phase

Climatic changes impact fruit tree growth and severely limit their production. Investigating the tree ability to cope with environmental variations is thus necessary to adapt breeding and management strategies in order to ensure sustainable production. In this study, we assessed the genetic parameters and genotype by environment interaction (GxE) during the early tree growth. One hundred and twenty olive seedlings derived from the cross ‘Olivière’ x ‘Arbequina’ were examined across two sites with contrasted environments, accounting for ontogenetic trends over three years. Models including the year of growth, branching order, environment, genotype effects, and their interactions were built with variance function and covariance structure of residuals when necessary. After selection of a model, broad sense heritabilities were estimated. Despite strong environmental effect on most traits, no GxE was found. Moreover, the internal structure of traits co-variation was similar in both sites. Ontogenetic growth variation, related to (i) the overall tree form and (ii) the growth and branching habit at growth unit scale, was not altered by the environment. Finally, a moderate to strong genetic control was identified for traits at the whole tree scale and at internode scale. Among all studied traits, the maximal internode length exhibited the highest heritability (H² = 0.74). Considering the determinant role of this trait in tree architecture and its stability across environments, this study consolidates its relevance for breeding.