evaluation in American Angus Daniela Lourenco S. Tsuruta, B. - - PowerPoint PPT Presentation

evaluation in american angus
SMART_READER_LITE
LIVE PREVIEW

evaluation in American Angus Daniela Lourenco S. Tsuruta, B. - - PowerPoint PPT Presentation

Issues in commercial application of single-step genomic BLUP for genetic evaluation in American Angus Daniela Lourenco S. Tsuruta, B. Fragomeni, Y. Masuda, I. Pocrnic, I. Aguilar, J.K. Bertrand, D. Moser, I. Misztal University of Georgia


slide-1
SLIDE 1

Issues in commercial application of single-step genomic BLUP for genetic evaluation in American Angus

Daniela Lourenco

  • S. Tsuruta, B. Fragomeni, Y. Masuda, I. Pocrnic, I. Aguilar, J.K. Bertrand,
  • D. Moser, I. Misztal

University of Georgia

July, 2016

slide-2
SLIDE 2

ssGBLUP theory vs. practice

Practice

  • Full capability
  • Implementation for American Angus
  • Challenges and problems
  • Studies on several livestock species
  • Simple

Theory

2

slide-3
SLIDE 3

Angus Data

  • 9.7M pedigree
  • BW, WW, PWG
  • SC
  • CE
  • Docility
  • Heifer pregnancy
  • Yearling height
  • Mature weight, height
  • Carcass weight, marbling, ribeye area, fat thickness
  • Dry matter intake

82,000 112,000 132,000 152,000 220,000 50,000 100,000 150,000 200,000 250,000 Number of Genotypes # Genotyped Animals

American Angus

07 2014 01 2015 07 2015 10 2015 07 2016

  • 220,000 genotyped animals

3

slide-4
SLIDE 4

ssGBLUP - lots of genotyped animals

CORE animals randomly sampled from genotyped population

0.97 0.99 0.99 0.99 5k 10k 15k 20k

Correlation (GEBV,GEBV_APY)

Growth Traits – 82k CORE invert NON-CORE

GAPY −1

Misztal et al., 2014

4

slide-5
SLIDE 5

ssGBLUP - lots of genotyped animals

  • How to choose number of core animals?
  • Ne, Me, ESM, Eigen of G
  • Limited dimensionality

Pocrnic et al., 2016 Misztal, 2016

5000 10000 15000 90 95 98 99 3654 6166 10605 14555 Number of Eigenvalues % of Variance

AAA – 82k

3654, 0.96 6166, 0.98 10605, 0.99 14555, 0.99 0.90 0.92 0.94 0.96 0.98 1.00 2000 4000 6000 8000 10000 12000 14000 16000 COR (GEBV,GEBV_APY) NUMBER OF EIGENVALUES

AAA - 82K

5

slide-6
SLIDE 6

ssGBLUP – Set core animals for evaluation

  • 220,000 genotyped animals
  • 99% eigen
  • Core = 19,000 high accuracy
  • Keep core constant 1 year
  • Add extra core
  • Genomic set up = 30min

CORE invert NON-CORE

19,000 201,000

6

slide-7
SLIDE 7

Growth Traits with external EBV

  • 220,000 genotyped animals
  • 19,000 core
  • 201,000 non-core
  • 9.7M pedigree
  • 7.4M BW
  • 8.1M WW
  • 4M PWG
  • Maternal effect for BW & WW
  • E = external
  • I = internal
  • T = PEV for E

Adapted from Legarra et al., 2007

  • External EBV from ~10k Red Angus

Computing time BLUP = 8h ssGBLUP = 12h

7

slide-8
SLIDE 8

Calving Ease is categorical

  • 220,000 genotyped animals
  • 19,000 core
  • 201,000 non-core
  • 8.7M Pedigree
  • 1.4M CE
  • 91% easy
  • 9% difficult
  • 2-trait BW-CE linear-threshold
  • BLUP = 12h
  • ssGBLUP = 4.5 days

8

slide-9
SLIDE 9

Scenario Description of parameters rounds hours correlation with genomic pcg rounds alpha beta traditional 40

  • 60

12

  • genomic

40 0.9 0.1 488 108

  • 1

100 0.9 0.1 81 43 0.999 2 100 0.85 0.15 62 32 0.999 3 200 0.9 0.1 24 25 0.999 4 200 0.85 0.15 19 19 0.999

𝐈−𝟐 = 𝐁−𝟐 + 𝟏 𝟏 𝟏 (α𝐇 + β𝐁𝟑𝟑)−𝟐– 𝐁𝟑𝟑

−𝟐

Calving Ease is categorical

Working on OMP – 30% faster

  • Cblup90iod2: 2 nested rounds

9

slide-10
SLIDE 10

Interim GEBV

  • 220,000 genotyped animals
  • 19,000 core
  • 201,000 non-core
  • Weekly evaluation
  • New genotypes daily

di = SNP weight = I

  • SNP effect
  • GEBVI = 𝐚

𝑣

  • COR (GEBVI_CORE,GEBVI_50k) = 0.98

10

slide-11
SLIDE 11

Accuracy of GEBV

  • GEBV published with accuracy
  • Measurement of precision

Diag(CZZ+) = PEV Traditional 𝑀𝐼𝑇𝑣𝑣

𝑗𝑗 =

1 (λ + 𝑒𝑗

𝑠 + 𝑒𝑗 𝑞)

  • Large datasets
  • Impossible to invert
  • 𝑒𝑗

𝑠and 𝑒𝑗 𝑞are approximated

(Misztal and Wiggans, 1988)

  • Accuracy = 1 - 𝑀𝐼𝑇-1

11

slide-12
SLIDE 12

Accuracy of GEBV

Genomic 𝑀𝐼𝑇𝑣𝑣

𝑗𝑗 =

1 (λ + 𝑒𝑗

𝑠 + 𝑒𝑗 𝑞 + 𝑒𝑗 𝑕)

Z′Z+ λA−1+ λ G−1−A22 −1

𝑒𝑗

𝑠

𝑒𝑗

𝑞

𝑒𝑗

𝑕

How to approximate 𝑒𝑗

𝑕?

12

slide-13
SLIDE 13

Accuracy of GEBV

  • Approximation 1
  • h2
  • # effective SNP
  • Average contribution from G
  • Accuracy of PA

𝑒𝑗

𝑕 = 1 + ℎ2 ∗ 𝐹𝑇𝑁 ∗ 𝐻 − 𝐵22 ∗ 𝐵𝑑𝑑𝑄𝐵

  • Approximation 2
  • 𝑒𝑗

𝑠 and 𝑒𝑗 𝑞

  • G
  • Minimum relationship (δ)
  • Accuracy of PA

13

𝑒𝑗

𝑕 = λ

𝑗≠𝑘 𝐵𝑑𝑑𝑄𝐵 𝑑𝑝𝑣𝑜𝑢(𝜀) +≈ 𝑗≠𝑘 𝑕𝑗𝑘(𝑒𝑗

𝑠 + 𝑒𝑗 𝑞)

𝑑𝑝𝑣𝑜𝑢(𝜀)

slide-14
SLIDE 14

Accuracy of GEBV

Approximation 1 Approximation 2

3’24” 0’12”

14

slide-15
SLIDE 15

Considerations

Several problems and challenges ssGBLUP ready for national Beef cattle evaluation

  • Angus in 2016
  • multi-trait, categorical, maternal
  • external info
  • interim GEBV
  • accuracy of GEBV

All solved

15

slide-16
SLIDE 16

Acknowledgements

16