NZ579456A - Biological markers on Bovine Chromosome 14 for size and uses therefor - Google Patents

Biological markers on Bovine Chromosome 14 for size and uses therefor

Info

Publication number
NZ579456A
NZ579456A NZ57945610A NZ57945610A NZ579456A NZ 579456 A NZ579456 A NZ 579456A NZ 57945610 A NZ57945610 A NZ 57945610A NZ 57945610 A NZ57945610 A NZ 57945610A NZ 579456 A NZ579456 A NZ 579456A
Authority
NZ
New Zealand
Prior art keywords
animal
proteins
fragments
isoforms
precursors
Prior art date
Application number
NZ57945610A
Inventor
Richard John Spelman
Bevin Lyal Harris
Michael Dominic Keehan
Stephen Richard Davis
Matthew Douglas Littlejohn
Michel Alphonse Jumien Georges
Wouter Robert Herman Coppieters
Latifa Karim
Tom Jean-Marc Druet
Haruko Takeda
Li Lin
Juan Antonio Castro Arias
Original Assignee
Livestock Improvement Corp Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Livestock Improvement Corp Ltd filed Critical Livestock Improvement Corp Ltd
Priority to NZ57945610A priority Critical patent/NZ579456A/en
Publication of NZ579456A publication Critical patent/NZ579456A/en

Links

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

Provided is a method for inferring the size potential of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the size potential of the animal and/or its offspring.

Description

*10059240376* <5 Wlf INTELLECTUAL PROPERTY OFFICE OF N.2. 3 1 AUG 2010 RECEIVED PATENTS FORM NO. 5 Our ref: SHR507647NZPR NEW ZEALAND PATENTS ACT 1953 COMPLETE SPECIFICATION BIOLOGICAL MARKERS AND USES THEREFOR We, LIVESTOCK IMPROVEMENT CORPORATION LIMITED, a New Zealand company, of Cnr Ruakura and Morrinsville Roads, Hamilton 3216, New Zealand, hereby declare the invention, for which we pray that a patent may be granted to us and the method by which it is to be performed, to be particularly described in and by the following statement: (This page is followed page 1A) 301494760_1 .DOC:SHR:qakld SHR507647NZPR 301525122 - 1A- BIOLOGICAL MARKERS AND USES THEREFOR FIELD The present invention relates generally to methods for inferring the size potential of an 5 animal and/or its offspring, particularly but not exclusively, methods for identifying and selecting animals on the basis of their liveweight and/or growth rate potential. The invention also relates to biological markers suitable for use in such methods.
BACKGROUND Traditional breeding programmes for livestock have involved selecting animals based on the presence of certain phenotypic characteristics. Recent methods have incorporated selection for genotypic characteristics known to be associated with improved production traits. Selection of animals on the basis of genotypic characteristics allows for earlier and more specific selection of animals of interest than does selection on the basis of phenotypic 15 characteristics.
The size of an animal may have an effect on its value for a particular purpose. For example, for beef farming it is desirable for animals to have a relatively larger weight or growth rate, whereas for dairy farming efficiencies can be gained with smaller animals, as 20 less feed is required to support milk production.
A number of researchers have reported an association between certain genetic regions, Quantitative Trait Loci (QTL), and/or certain polymorphisms and weight in bovine species (Mizoshita et al. (2004), Takasuga et al. (2007), Buchanan et al. (2005), Kneeland et al. 25 (2004), Nkrumah et al. (2007), Mizoshita et al. (2005)). Gudbjartsson et al., (2008) and Lettre et al., (2008) have reported certain genetic regions being associated with height in humans. However, there is a need to identify specific markers and genotypes which will allow for improvements in the ability to predict the size (including weight and/or height) potential of livestock, for breeding, selection, and farm management purposes for example.
Bibliographic details of the publications referred to herein are collected at the end of the description.
Received at IPONZ on 24.02.2012 SHR507647NZPR ^ 302016060 ^ OBJECT It is an object of the present invention to provide improved methods for selecting animals, for predicting the size potential of an animal and/or its offspring, of screening a population of animals, breeding animals, managing animals, and/or estimating the worth of an animal 5 and/or its offspring, or to at least to provide the public with a useful choice.
STATEMENT OF INVENTION The inventors have identified a number of genetic markers whose sequence or genotype can be used to predict or infer the potential size of an animal and/or its offspring. Such 10 information may be used in methods for selecting, screening and breeding animals, farm management, and for estimating an animals worth to a particular industry, for example.
In one aspect, the invention provides a method for inferring the size potential of an animal and/or its offspring, the method comprising at least the step of analysing a nucleic acid 15 from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in 20 linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the size potential of the animal and/or its offspring.
In a particular embodiment, the invention provides a method for inferring the growth rate potential of an animal and/or its offspring, the method comprising at least the step of 25 analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648,23187238,23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in 30 linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the growth rate potential of the animal and/or its offspring.
In a particular embodiment, the invention provides a method for inferring the liveweight potential of an animal and/or its offspring, the method comprising at least the step of Received at IPONZ on 24.02.2012 SHR507647NZPR q 302016060 ^ analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260,23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734,23264810, 23264854, 23272034, 5 and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the liveweight potential of the animal and/or its offspring.
In one embodiment, the presence of one or more of 22986260T, 23186380T, 23186648A, 10 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a larger size for the animal and/or its offspring. In one embodiment, the presence of one or more of these markers infers a higher growth rate and/or higher liveweight for the animal and/or its offspring.
In one embodiment, the presence of one or more of 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C infers a smaller size for the animal and/or its offspring. In 20 one embodiment, the presence of one or more of these markers infers a lower growth rate and/or lower liveweight for the animal and/or its offspring.
In another aspect the invention provides a method for selecting or rejecting an animal, the method comprising at least the step of analysing a nucleic acid from said animal to 25 determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus\ 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more 30 thereof, wherein the nucleotide present at the one or more positions infers the size potential of the animal and/or its offspring.
In one embodiment, the method comprises at least the steps of: SHR50764 7NZPR 302086287 received 16 april 2012. 4 a) analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989,23215765,23218552,23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof; and, b) selecting or rejecting an animal based on the nucleotide present at the one or more positions.
In one embodiment, the presence of one or more of 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(delction), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a larger size for the animal and/or its 15 offspring. In one embodiment, the presence of one or more of these markers infers a higher growth rate and/or higher liveweight for the animal and/or its offspring.
In one embodiment, the presence of one or more of 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 20 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C infers a smaller size for the animal and/or its offspring. In one embodiment, the presence of one or more of these markers infers a lower growth rate and/or lower liveweight for the animal and/or its offspring.
In another aspect, the invention provides a method of identifying an animal which is more likely or less likely to produce one or more desirable trait, the method comprising at least the step of analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238,23205280,23211989, 23215765, 30 2321.8552,23228064, 23232324, 23235239, 23243588,23243734,23264810,23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the animal is more likely or less likely to produce the one or more desirable trait.
In one embodiment, the method comprises at least the steps of; Received at IPONZ on 24.02.2012 SHR507647NZPR c 302016060 ^ a) analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552,23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof; and, b) identifying the animals which are more likely or less likely to produce one or more desirable trait on the basis of the nucleotide present at the one or more positions.
In one embodiment, the one or more desirable trait is size. In one embodiment, the trait is growth rate. In one embodiment, the trait is liveweight.
In one embodiment, animals which are more likely to have a smaller size are identified on the basis of the presence of one or more of 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C. In one embodiment, animals which are more likely to have a 20 lower growth rate and/or lower liveweight are identified on the basis of the presence of one or more of these markers.
In one embodiment, animals which are more likely to have a larger size are identified on the basis of the presence of one or more 22986260T, 23186380T, 23186648A, 25 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T. In one embodiment, animals more likely to have a higher growth rate and/or higher liveweight are identified on the basis of the presence of one or more of these markers.
In one particular embodiment, the invention provides a method for identifying an animal having a higher growth rate potential and/or higher liveweight potential, the method comprising at least the step of analysing a nucleic acid from said animal to determine the SHR507647NZPR 302016060 Received at IPONZ on 24.02.2012 nucleotide sequence at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989,23215765,23218552,23228064, 23232324, 23235239, 23243588,23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the presence of one or more of 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a higher growth rate and/or higher liveweight potential.
In another particular embodiment, the invention provides a method for identifying an animal having a lower liveweight potential and/or a lower growth rate potential, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide sequence at one or more positions corresponding to the following 15 positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280,23211989,23215765,23218552,23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the presence of one or more of 22986260C, 23186380C, 23186648G, 20 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C infers a lower liveweight and/or lower growth rate potential.
In another aspect, the invention provides a method for estimating the worth of an animal 25 and/or its offspring, the method comprising at least the step of analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, Received at IPONZ on 24.02.2012 SHR5Q7647NZPR « 302016060 ' and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present infers of the size potential of the animal and/or its offspring.
In particular embodiment, the invention provides a method for estimating the worth of an animal and/or its offspring, the method comprising at least the step of analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 10 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present infers of the growth rate potential of the animal and/or its offspring.
In one particular embodiment, the invention provides a method for estimating the worth of an animal and/or its offspring, the method comprising at least the step of analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 20 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present infers of the liveweight potential of the animal and/or its offspring.
In one embodiment, the presence of one or more of 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a larger size potential of the animal and/or its offspring. In one embodiment, the presence of one or more of these markers infers a 30 higher growth rate and/or higher liveweight for the animal and/or its offspring.
In one embodiment, the present of one or more of 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, Received at IPONZ on 24.02.2012 SHR507647NZPR q 302016060 O 23272034G, and 23273471C infers a smaller size potential of the animal and/or its offspring. In one embodiment, the presence of one or more of these markers infers a lower growth rate and/or lower liveweight for the animal and/or its offspring.
In another aspect, the invention provides a method for breeding animals to produce offspring, which comprises selecting at least a first animal that has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 2318723 8A} 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG) 11, 23264854G, 23272034A, and 10 23273471T, wherein the position of the genetic marker is read in relation to chromosome 14 of Bos taurus, and mating said first animal with a second animal to produce offspring.
In one embodiment, the second animal is selected on the basis it has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 15 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 2327347IT.
In another broad aspect, the invention provides a method for breeding animals to produce 20 offspring, which comprises selecting at least a first animal that has one or more of the following genotypes 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C, wherein the position of the genetic marker is read in relation to chromosome 14 of Bos 25 taurus, and mating said parent with a second animal to produce offspring.
In one embodiment, the second animal is selected on the basis it has one or more of the following genotypes 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 30 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C.
Received at IPONZ on 24.02.2012 SHR507647NZPR q 302036060 " In one embodiment, the methods of the invention further involve taking a sample from the animal.
In one embodiment, the methods further involve analysing one or more additional 5 biological markers.
Preferably, analysis of the nucleotide sequence of the one or more genetic markers occurs using one or more of: polymerase chain reaction (PCR); gel electrophoresis; Southern blotting; nucleic acid sequencing; restriction fragment length polymorphism (RFLP); 10 single-strand confirmation polymphism (SSCP); LCR (ligase chain reaction); denaturing gradient gel electrophoresis (DGGE); allele-specific oligonucleotides (ASOs); proteins which recognize nucleic acid mismatches; RNAse protection; oligonucleotide array hybridisation; denaturing HPLC (dHPLC); high resolution melting (HRM); and, matrix-assisted laser desorption/ionisation time-of-flight mass spectroscopy (MALDI-TOF MS).
In another aspect the invention provides a method for inferring the size potential of an animal and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more 20 proteins, precursors, isoforms, or fragments in said animal, wherein the level infers the size potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In another aspect the invention provides a method for inferring the growth rate potential of 25 an animal and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal, wherein the level infers the growth rate potential of the animal and/or its offspring, and wherein the one or 30 more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK. received 16 april 2012.
SHR5G7647NZPR i 302086124 A U In another aspect the invention provides a method for inferring the liveweight potential of an animal and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one 5 or more proteins, precursors, isoforms, or fragments in said animal., wherein the level infers the liveweight potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CIICIID7, RDHE2, SDR16C6, and PENK, In another aspect, the invention provides a method of identifying an animal which is more or less likely to produce one or more desirable trait, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, 15 wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK, wherein the level infers the animal is more likely or less likely to produce the one or more desirable trait. In one embodiment, the one or more desirable trait is size. In one embodiment, the trait is growth rate. In. one embodiment, the trait is liveweight.
In one embodiment, the methods involve comparing the level of expression of one or more of said proteins, precursor thereof, isoforms thereof, fragments thereof, and/or nucleic acids encoding said proteins against a standard, wherein a difference in the level of expression in the sample compared to the standard infers the size potential, growth rate 25 potential and/or liveweight potential of the animal and/or its offspring. In one embodiment, an increase in the level of expression infers a higher growth rate and/or higher liveweight potential for the animal and/or its offspring.
In one particular embodiment, the invention provides a method for identifying an animal 30 having a higher growth rate potential and/or a higher liveweight potential, the method comprising at least the step of observing the level of expression of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal Received at IPONZ on 24.02.2012 SHR507647NZPR i i 302016060 A A , wherein an increase in the level of expression compared to a standard infers a higher growth rate potential and/or higher liveweight potential, and wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In another particular embodiment, the invention provides a method for identifying an animal having a lower growth rate potential and/or a lower liveweight potential, the method comprising at least the step of observing the level of expression of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more 10 fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal, wherein a decrease or substantially no increase in the level of expression compared to a standard infers a lower growth rate potential and/or lower liveweight potential, and wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In another broad aspect the invention provides a method for selecting or rejecting an animal the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, 20 precursors, isoforms, or fragments in an animal and selecting or rejecting an animal based on the level of one or more of said proteins, precursors, isoforms or fragments thereof, and/or nucleic acids encoding same, wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In one embodiment, the method involves comparing the level of expression of one or more of said proteins, precursors thereof, isoforms thereof, fragments thereof, and/or nucleic acids encoding said proteins against a standard, wherein a difference in the level of expression in the sample compared to the standard infers the size potential, growth rate potential and/or liveweight potential of the animal and/or its offspring. In one 30 embodiment, an increase in the level of expression infers a higher growth rate and/or higher liveweight potential for the animal and/or its offspring. In one embodiment, a decrease or substantially no increase in the level of expression compared to the standard Received at IPONZ on 24.02.2012 SHR507647NZPR ^ â– *) 302016060 infers a lower growth rate potential and/or lower liveweight potential for the animal and/or its offspring.
In one embodiment, the method is performed for the purpose of selecting or rejecting an 5 animal for milking purposes. In one embodiment, the method is performed for the purpose of selecting or rejecting an animal for beef farming. In another embodiment, the method is performed for the purpose of selecting or rejecting an animal for breeding purposes.
In one embodiment, where the method is performed for the purpose of selecting or 10 rejecting an animal for milking purposes, the animal is selected if there is a decrease or no substantial increase in the level of expression compared to the standard or the animal is rejected where there is an increase in the level of expression compared to the standard.
In one embodiment, where the method is performed for the purpose of selecting an animal 15 for beef farming, the animal is selected if there is an increase in the level of expression compared to the standard or the animal is rejected where there is a decrease or no substantial increase in the level of expression compared to the standard.
In another aspect the invention provides a method for estimating the worth of an animal 20 and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal, wherein the level infers the size potential of the animal and/or its offspring, and wherein the one or more proteins are selected from 25 RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In one embodiment, the level infers the growth rate potential and/or the liveweight potential of the animal and/or its offspring.
In one embodiment, the method involves comparing the level of expression of one or more of said proteins, precursor thereof, isoforms thereof, fragments thereof, and/or nucleic acids encoding said proteins against a standard, wherein a difference in the level of expression in the sample compared to the standard infers the size potential, growth rate Received at IPONZ on 24.02.2012 SHR50764TNZPR -| o 302016060 A potential and/or liveweight potential of the animal and/or its offspring. In one embodiment, an increase in the level of expression infers a higher growth rate potential and/or higher liveweight potential for the animal and/or its offspring. In one embodiment, a decrease or substantially no increase in the level of expression compared to the standard 5 infers a lower growth rate potential and/or lower liveweight potential for the animal and/or its offspring.
In another aspect, the invention provides a method for breeding animals to produce offspring, the method comprising observing the level of one or more proteins, one or more 10 precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding one or more proteins, precursors, isoforms, or fragments in at least a first animal, selecting the first animal where there is an increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, and mating said first animal with a 15 second animal to produce offspring, wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In one embodiment, the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more 20 fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in the second animal, and selecting the second animal where there is an increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, 25 SDR16C6, and PENK.
In another aspect, the invention provides a method for breeding animals to produce offspring, the method comprising observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or 30 one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in at least a first animal, selecting the first animal where there is a decrease or Received at IPONZ on 24.02.2012 SHR507647NZPR 1 'I A 302016060 AJ/i substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, and mating said first animal with a second animal to produce offspring, wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, 5 and PENK.
In one embodiment, the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, 10 precursors, isoforms, or fragments in the second animal, and selecting the second animal where there is a decrease or substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more proteins are selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK.
In one embodiment the methods comprise at least the steps of: a) taking a sample from an animal; b) detecting one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding any one or more thereof in the sample, wherein the one or more protein is selected from RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6, and PENK; and, c) comparing the level of the one or more proteins, precursors, fragments, isoforms and/or nucleic acids against a standard.
Preferably at least an approximately 1.2 fold increase in the level of one or more of the proteins, a precursor thereof, an isoform thereof, a fragment thereof, and/or a nucleic acid encoding same, is indicative of a higher growth rate potential and/or a higher liveweight potential of the animal and/or its offspring.
In one embodiment the methods further involve analysis of one or more additional biological markers.
Received at IPONZ on 24.02.2012 SHR5Q7647NZPR 1 I'D 302016060 LJD In certain embodiments of the invention the level of the one or more proteins, precursors, fragments, isoforms and nucleic acids encoding same is determined using an immunoassay, separation based on characteristics such as molecular weight and isoelectric point, gel electrophoresis, Western Blotting or mass spectroscopy. Preferably the 5 immunoassay is an ELISA. Preferably the gel electrophoresis is 2D gel electrophoresis or gel-free systems based on microfluidics technologies.
SHR507647NZPR 301525122 In one embodiment of the methods of the invention the animal is bovine. In a particular embodiment the bovine animal is Bos taurus or Bos indicus. In a particular embodiment, the animal is chosen from the group consisting Jersey, Holstein, Friesian or crossbred dairy 5 cattle. In another embodiment, the animal is Simmental.
The invention may also be said broadly to consist in the parts, elements and features referred to or indicated in the specification of the application, individually or collectively, in any or all combinations of two or more of said parts, elements or features, and where 10 specific integers are mentioned herein which have known equivalents in the art to which the invention relates, such known equivalents are deemed to be incorporated herein as if individually set forth.
FIGURES These and other aspects of the present invention, which should be considered in all its novel aspects, will become apparent from the following description, which is given by way of example only, with reference to the accompanying figures, in which: Figure 1: Structure of the Holstein-Friesian (HF) x Jersey (J) intercross generated for QTL 20 mapping. Following conventional genetic formalism matings are described mentioning females first and males second. Hence "HF x J" corresponds to a mating between HF dams and J sires. Males are labeled with squares/rectangles, females with circles/ovals. Number of individuals of a given genotype are given in parentheses.
Figure 2: (A) Location scores (F-values) obtained across the entire bovine genome when 25 analyzing weight at birth, 6 (6M), 12 (12 M), 18 (18 M) and 24 months (24M) and height at 18 months in the HFxJ inter-cross with a line-cross model (Haley et al., 1994). (B) Location scores (F-values) obtained for BTA14 when analyzing weight and height in the HFxJ inter-cross with a paternal half-sib pedigree model implemented with HSQM (Coppieters et al., 1998). The lines correspond to weight at birth, 6,12,18 and > 24 (=live 30 weight) months and "height" as labelled. The black dotted line corresponds to the 5% genome-wide significance threshold determined with a permutation test (Churchill & Doerge, 1994). The lower horizontal bar (CI) corresponds to the 95% CI for the QTL location determined by bootstrapping (Visscher et al., 1996). The higher horizontal bar (SIAG) shows the position of the 1.1 MB critical interval defined by Mizoshita et al.
SHR507647NZPR 301525122 (2005). (C) Highest chromosome-wide log(l/p) values obtained within each of the six sire families for height (white bars) and live weight (black bars). The corresponding map positions are given above the bars. Note the very distinct map positions for sire 4. The black dotted line marks the 5% chromosome-wide significance threshold. The black arrows point towards the four sire families segregating for the QTL. The corresponding sires (1, 2, 3 and 5) are therefore predicted to be of "Qq" genotype. (D) Sire-specific allele substitution effects on weight (gray bars of increasing darkness: birth, 6,12,18 months and live weight) and height (white bars) expressed in Kgs and mm, respectively. Family-specific allele substitution effects were determined at the most significant QTL position as determined in the across-family analysis.
Figure 3: (A) Black line (labeled "HSQM"): Linkage-based QTL analysis of live weight in the HFxJ intercross using HSQM (Coppieters et al., 1998) and a 56 microsatellite marker map (cfr. Fig. 2B). Results are expressed as F-values (right axis). Dark dots (labeled "SP-T"): Single-point linkage + LD analysis of live weight in the HFxJ intercross 15 using Dualphase (Druet & Georges, 2009) and a 925 SNP + 56 microsatellite marker map. Results are expressed as a likelihood ratio (LR) test (21nLR)(left axis). (B) Black line (labeled "MP-T"): Haplotype-based linkage + LD analysis of live weight in the HFxJ intercross using Dualphase and a 925 SNP + 56 microsatellite marker map. Results are expressed as a likelihood ratio (LR) test (21nLR)(left axis). (C) Black line (labeled 20 "OUTBRED"): Haplotype-based linkage + LD analysis of breeding value for live weight in 2,700 progeny-tested sires from the New Zealand outbred dairy cattle population using Dualphase and 293 SNP markers from the USDA 50K SNP chip (Van tassel et al., 2008) spanning a 15 Mb BTA14 segment spanning the QTL location. Results are expressed as a likelihood ratio (LR) test (21nLR)(left axis). (D) Dark dots (labeled "TAQ"): Single-point 25 linkage + LD analysis of live weight in the HFxJ intercross using Dualphase and 11 candidate QTN identified by resequencing the 750 Kb critical interval in the six F1 sires. (A-D) The X-axis corresponds to the chromosomal position in base pairs. The lower horizontal line corresponds to the QTL critical defined by Mizoshita et al. (2005), the higher horizontal line to the 750 Kb segment sequenced in the present study. Every graph 30 shows the results of all analysis in gray watermarks to facilitate cross comparison.
Figure 4: (A) Dark gray line: Single QTL haplotype-based combined L+LD analysis of live weight in the HFxJ intercross (cfr. Fig. 3B). Black line: Two QTL haplotype-based combined L+LD analysis of live weight conducted with the same mixed model augmented SHR507647NZPR 301525122 with a random effect corresponding to the hidden state at the most likely QTL position obtained in the single QTL analysis (black arrow). The black line measures the increase in LRT (above that of the single-QTL analysis at the position of the black arrow) by adding a second QTL at the corresponding position. (B) Dark gray line: Single QTL haplotype-5 based combined L+LD analysis of breeding value for live weight in the outbred New Zealand dairy cattle population (cfr. Fig. 3C). Black line: Two QTL haplotype-based combined L+LD analysis of breeding value for live weight in the outbred New Zealand dairy cattle population with the same mixed model augmented with a random effect corresponding to the hidden state at the most likely QTL position obtained in the single 10 QTL analysis (black arrow). The black line measures the increase in LRT (above that of the single-QTL analysis at the position of the arrow) by adding a second QTL at the corresponding position. Modest evidence (LRT: 20 = lod score: 4.3) for an additional QTL within the interval defined by Mizoshita et al. (2005) is observed (gray arrow).
Figure 5: Effect on live weight (in Kgs; X-axis) and frequency (number of observations in 15 the studied population; Y-axis) for the 20 hidden haplotype states modeled with Dualphase (Druet and Georges, 2009) in the New Zealand outbred dairy cattle population. Shades of gray distinguish the breed origin of the corresponding animals: Holstein-Friesian (black), Jersey (gray), crossbred (white). The number of haplotype states in each class is given above the corresponding bars.
Figure 6: Upper track (labeled "Genes"): A. Organization of the eight genes mapping to the 750Kb critical region (LYN, RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6 and PENK). Middle track (labeled "SNP"): Positions of the 14 candidate Quantitative Trait Nucleotides (QTN), plus splice site variant detected in CHCHD7 (half height). Brown track (labeled "Phastcons conserved elements, 5-way Vertebrate Multiz Alignment"): 25 Location of "Phastcons" multispecies conserved elements. (B) Orthologous locus in zebrafish (D. rerio) Figure 7: Combined linkage + LD analysis of weight at birth (A), weight at 6 months (B), weight at 12 months (C), weight at 18 months (D) and height at 18 months (E) using 925 SNPs + 56 microsatellites in single point analysis (gray dots), 925 SNPs + 56 30 microsatellites in haplotype-based analysis (grey line), and 10 of the 14 candidate QTN (black circled dots, labeled). The higher and lower horizontal lines respectively mark (i) SHR507647NZPR 301525122 the 750 Kb critical region sequenced in this work, and (ii) the interval defined by Mizoshita et al. (2005).
Figure 8: (A) Identification of the Q (grey), q (darker grey), and "recombinant" (grey + darker grey) haplotypes in a bovine diversity panel. Breeds and corresponding numbers of 5 genotyped animals are given. Breed-specific haplotype frequencies are given and highlighted with a black box on white text when > 0.03. (B) Haplotype effect on birth weight in Simmental, assuming that the "R " haplotype (cfr. A) is q (left panel) or Q (right panel). Animals are sorted by predicted QTL genotype as shown. Smaller dots (with a shadow) correspond to individual observation. Larger dots correspond to genotype means.
Arrows correspond to the 95% confidence intervals of the genotype mean (mean ± 1.96*SEM). The regression of birth weight on QTL genotype is near significant (p=0.06) when assuming that R=q (left panel) and non-significant under the alternative hypothesis (R=Q). The difference in birth weight between animals predicted to be of qq and QQ genotype assuming R=q is significant (t-test; p=0.02).
Figure 9: Exploring the integrity of the PLAG1 (A-D), CHCHD7 (E-F) and MOS (G-H), ORFs by RT-PCR performed with cDNA prepared from bone (Bo), muscle (M), brain (Br) and liver (L) mRNA extracted from fetuses that were homozygous QQ or qq for the pQTN. Animals that carried the FJX_303486_1:1 CHCHD7 splice site variant (cfr. SUppl. Fig. 11) are indicated with "sp". (1) Negative controls.
Figure 10: Effect of pQTN genotype on the expression level of the eight genes mapping to 750-kb CI for the QTL in fetal liver, bone, muscle and brain, estimated by QRTPCR (mid-darkness bars), an allelic imbalance test using 3'UTR SNPs (darkest bars) or an allelic imbalance test using an intronic SNP (lightest bars). The X-axis measures the slope of the regression (QRTPCR) or the ratio of Q over q allele (allelic imbalance tests) on a log2-scale. The vertical black lines correspond to the absence of an effect of pQTN genotype on expression. #: p<0.10; *:0.01 Figure 11: (A) QRT-PCR results for LYN, RPS20, MOS, PLAG1, CHCHD7, RDHE2, SDR16C6 and PENK and MOS in fetal bone, muscle, brain and liver samples. Small grey circles correspond to the expression level in a given animal relative to the mean expression level across all animals. Relative expression levels are expressed on a log2-scale, and SHR507647NZPR 301525122 correspond to the mean relative expression level estimated with one or more amplicons (the number of tested amplicons is given in the upper left corner of each graph). Animals were sorted by pQTN genotype as indicated at the top of the figure (qq, Qq and QQ). Large darker circles correspond to the mean relative expression of the corresponding 5 pQTN genotype, "a" corresponds to the slope of the regression line (in log2-units). p corresponds to the statistical significance of the regression. (B) Results of the allelic imbalance test in heterozygous Qq animals for PLAGl and CHCHD7. The X-axis corresponds to the ratio of peak heights corresponding to the Q and q allele as measured by direct sequencing and Peakpicker (Ge et al., 2005). Small gray circles correspond to the 10 allelic ratio for individual animals, large darker circles correspond to the overall mean. Numbers on the right of the graphs correspond to the mean value +_ 1.96 x SEM, and (for cDNAs) to the p-value of the difference with genomic DNA (gDNA) evaluated by a t-test. Fetuses that were heterozygous for the CHCHD7 FJX_303486_1:1 splice site variant were excluded from both QRT-PCR and allelic imbalance analyses of PLAGl and CHCHD7.
Figure 12: Testing the effect of pQTN FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l on bidirectional promoter strength using a luciferase reporter assay. (A) Schematic representation of the supposedly bidirectional promoter driving expression of the head-to-head oriented PLAGl and CHCHD7 genes (dark line at top of figure), with corresponding Phastcons conservation scores (lighter bars), and multispecies 20 sequence alignment of a segment encompassing the FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l candidate pQTN. The black arrows mark the positions of the "long" and "short" fragments cloned in the pGL4 luciferase reporter vector in the "forward" (towards PLAGl) and "reverse" (towards CHCHD7) orientation. (B) Ratios of firefly to renila luminescence obtained after transfection of Cosl cells with (i) a 25 promoterless pGL4 vector, (ii) two distinct, sequence-verified preparations of the pGL4 vector endowed with the TK promoter, (iii) pairs of sequence-verified preparations of the vector endowed with the q or Q version, of the long or short fragment, cloned either (as denoted under figure) in forward or reverse orientation. Error bars correspond to standard errors of the mean computed from replicates. (C) Schematic representation of the 30 "recombinant" "Q-q" and "q-Q" promoter fragments that were generated by swapping the Q and q residues at the FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l sites as shown. (D) Firefly/Renilla luminescence obtained with the non-recombinant « Q-Q » promoter as well as recombinant « Q-q » ane «q-Q» promoters cloned in forward and reverse SHR507647NZPR 301525122 orientation, relative to the cognate non-recombinant « q-q » promoter. Error bars correspond to standard errors of the mean.
Figure 13: Representative results of EMS A experiments conducted with radiolabeled 29-5 mer (sQ and sq) and 74/80-mer probes (LQ and Lq) spanning the FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l pQTN. "sCtrl" corresponds to an unrelated 25-mer control duplex used as control competitor. Results shown were obtained with nuclear extracts derived from fetal bone (A) and C2C12 cells (B&C). Complexes with differential affinity for the Q and q allele are marked by the arrows on the left of the gels (short probe; Qq). An abundant complex with equal affinity for the Q and q allele detected with the long probe is marked by the green arrow. Free probes are labeled "SQ/q", "LQ" and "Lq". The bar graphs at the bottom of the figure quantify the abundance of two corresponding complexes (longer probe -darker bars, 15 shorter probe - lighter bars) relative to the Q probe (in the absence of cold competitor) determined by densitometry.
Figure 14: (A) Mapping eQTL influencing CHCHD7 expression levels in liver (lighter and upper line) and adipose (darker and lower line line) using a 56 microsatellite BTA14 20 map and HSQM (Coppieters et al., 1998). Location scores correspond to F-values. The black dotted line corresponds to the 5% chromosome-wide significance threshold determined with a permutation test (Churchill & Doerge, 1994). Inlet: Results of the analysis conducted within each of the six sire families separately. The p-values obtained on BTA14 for each of the six families is indicated by the height of the first (liver) or 25 second (adipose) bars, and indicate segregation of the eQTL in pedigree "1" only (arrow). (B) Haplotype cluster-based association analysis conducted as described in M&M for CHCHD7 expression in liver (lighter and upper line) and adipose (darker and lower line). Markers considered are the 925 BTA14 SNPs + 56 BTA14 microsatellites. Location scores correspond to log(l/p), where p corresponds to the p-value of the corresponding F-30 test. (C) Single SNP association analysis conducted as described in M&M for CHCHD7 expression in liver (lighter dots) and adipose (darker dots). Markers considered were the 925 BTA14 SNPs + 56 BTA14 microsatellites. The arrows point towards the p-values obtained for the CHCHD7 FJX_303486_1:1 splice site variant. (D) Two SNP SNP association analysis conducted as described in M&M for CHCHD7 expression in liver SHR507647NZPR 301525122 (lighter dots) and adipose (darker dots). Markers considered were the 925 BTA14 SNPs + 56 BTA14 microsatellites. Log(l/p) values correspond to the p-value of adding the corresponding SNP in a model already including the CHCHD7 FJX_303486_1:1 splice site variant. (A-D) The position of the sequenced 750 Kb interval containing CHCHD7 is 5 marked by the black horizontal line.
Figure 15: (A) Mapping eQTL influencing RPS20 expression levels in liver (lighter and upper line) and adipose (darker and lower line) using a 56 microsatellite BTA14 map and HSQM (Coppieters et al., 1998b). Location scores correspond to F-values. The black dotted line corresponds to the 5% chromosome-wide significance threshold determined 10 with a permutation test (Churchill & Doerge, 1994). Inlet: Results of the analysis conducted within each of the six sire families separately. The p-values obtained on BTA14 for each of the six families is indicated by the height of the first (liver) or second (adipose) bars, and indicate segregation of the eQTL in pedigree "1", "5" and "6" (arrows). (B) Haplotype cluster-based association analysis conducted as described in M&M for RPS20 15 expression in liver (upper and darker line) and adipose (lower and lighter line). Markers considered are the 925 BTA14 SNPs + 56 BTA14 microsatellites. Location scores correspond to log(l/p), where p corresponds to the p-value of the corresponding F-test. The position of the sequenced 750 Kb interval containing CHCHD7 is marked by the black horizontal line.
Figure 16: (A) Agarose gel of RT-PCR products obtained from AA, AT and TT individuals for the CHCHD7 FJX_303486_1:1 splice site variant. The homozygous mutant (AA) animal shows three bands shown (by sequencing) to correspond to mRNAs with structures D and E as depicted, i.e. having lost exon 3. The largest band correspond to a D/E heteroduplex. The homozygous wild-type (TT) animals show three bands corresponding 25 to mRNAs with structures A, B and C as depicted. The largest band is a mixture of B/C heteroduplex plus band A. Heterozygous (AT) animals show all fragments (A-E) plus additional heteroduplexes. The darker portions of the different mRNA species correspond to predicted ORF, while the lighter segments correspond to 5' and 3' UTRs. (B) Relative CHCHD7 transcript levels determined by QRT-PCR in muscle, brain and liver samples of 30 heterozygous AT and homozygous wild-type TT fetuses. Smaller dots correspond to individual expression levels, while the larger dots correspond to the genotypic average. Relative expression levels are expressed as log2 = (xy/x^ , where xtJ /3cy is the expression SKR507647NZPR 301525122 level of individual i measured with amplicon j relative to the other individuals, "a" is the slope of the regression (in log2 units), and p the significance of its deviation from 0.
PREFERRED EMBODIMENT(S) The following is a description of the present invention, including preferred embodiments 5 thereof, given in general terms. The invention is further elucidated from the disclosure given under the heading "Examples" herein below, which provides experimental data supporting the invention.
The inventors have identified that particular alleles of a number of genetic markers in a 10 region of chromosome 14 in Bos taurus correlate with the weight and/or height of an animal. The inventors note a particular correlation to weight at birth, at 6 months of age, 12 months of age, 18 months of age and 24 months of age. The genotype of one or more of these particular genetic markers may therefore be used to infer or predict the size potential of an animal and/or its offspring. The size of an animal may contribute to its 15 value to a particular industry. For example, typically a larger animal or an animal with a high growth rate is preferable for beef farming and efficiencies can be gained in dairy farming by using smaller animals which require less feed.
The inventors also note that there is a correlation between increases in the expression of 20 certain genes and the size of an animal. Thus, observing the level of expression of one or more of these genes may also be used to infer or predict the size potential of an animal and/or its offspring.
Thus, analysis of the genetic markers of the invention, or the level of expression of 25 particular genes, may assist in: predicting phentotypic performance; identifying animals more or less likely to have a desired trait; the selection or rejection of animals for breeding purposes; managing animals in order to maximise their individual potential performance and value; estimating the worth or economic value of an animal; improving profits related to selling animals and/or products produced from the animals; improving the genetics of a 30 population of animals by selecting and breeding desirable animals; cloning animals likely to have a specific trait; predicting the suitability of an animal and/or its offspring to use in different industries.
SHR507647NZPR 301525122 The term "genetic marker" as used herein refers to nucleic acids or specific genetic loci (including specific nucleotide positions) that are polymorphic or contain sequence variations within a population, the alleles of which can be detected and distinguished by one or more analytic methods. The term "genetic marker" further includes within its scope 5 a plurality of genetic markers co-segregating, in the form of a "haplotype". In this context, the term "haplotype" refers to a plurality of genetic markers that are generally inherited together. Typically, genetic markers within a haplotype are in linkage disequilibrium. Reference herein to analysing a nucleic acid to determine the nucleotide sequence of a "genetic marker" (for example the genetic markers of Table 1) or at a particular genetic 10 position should be taken to include analysing and determining the nucleotide sequence on either strand of the nucleic acid.
The term "single nucleotide polymorphism" (SNP) refers to nucleic acid sequence variations that occur when a single nucleotide in the genome sequence is altered. A single 15 nucleotide polymorphism may also be a single nucleotide insertion or deletion. The different nucleotides within a SNP are referred to as an allele.
"VNTR", or "variable number tandem repeat", refers to a tandem repeat of a nucleic acid sequence at a genetic locus in which the number of repeated DNA segments varies 20 between individuals in a population. The different repeat numbers within a VNTR at a particular locus may be referred to as an allele.
The term "genotype" as used herein means the genetic constitution or nucleotide sequence at one or more genetic locus, in particular the nucleotide sequence of an allele of a genetic 25 locus.
As used herein, "isoform(s)" should be taken broadly and includes those generated by posttranslational modification (for example, glycosylation), and proteolytic processing. The term "precursor(s)" should also be interpreted broadly and includes for example 30 preproteins, prepeptides, preproproteins and prepropeptides including an N-terminal leader sequence.
Reference herein to the term "size" should be taken broadly, and includes total size, weight (including liveweight), height, and/or growth rate.
SHR507647NZPR 301525122 "Growth rate" as used herein is intended to mean the average daily weight gain of an animal.
The phrase "infer size potential", and like phrases, should be taken broadly to mean that the presence of a particular allele of a genetic marker indicates an increased likelihood that an animal and/or its offspring has the capacity to have a certain size characteristic (at one or more time points in its life; for example, at birth, or at an age at which it will be put to productive use in a particular industry such as the beef or dairy industry) compared to if it 10 had a different allele of the genetic marker. The phrase should not be taken to imply that the actual size of an animal is predicted. In a preferred embodiment, the presence of a particular allele of a genetic marker of the invention (or a particular level of expression of a specified protein) provides an indication that the animal and/or its offspring is more likely than not to have a particular size characteristic relative to if it had a different allele. 15 In one embodiment of the invention, an animal with one copy of a Q allele of a genetic marker (as outlined herein) will have from an approximately 60% to approximately 70% probability of having a larger size than if it had a different allele of that genetic marker. In one particular embodiment, an animal with one copy of a Q allele of a genetic marker (as outlined herein) will have an approximately 65% probability of having a larger size than if 20 it had a different allele of that genetic marker. In one embodiment of the invention, an animal with two copies of a Q allele of a genetic marker (as outlined herein) will have from an approximately 60% to approximately 70% probability of having a larger size than if it only had a single copy of the Q allele of that genetic marker. In one particular embodiment, an animal with two copies of a Q allele of a genetic marker (as outlined 25 herein) will have an approximately 65% probability of having a larger size than if it only had a single copy of the Q allele of that genetic marker. In one embodiment, an animal with two copies of a Q allele of a genetic marker (as outlined herein) will have from an approximately 73% to approximately 83% probability of having a larger size than if it had two copies of a different allele of that genetic marker. In one particular embodiment, an 30 animal having two copies of a Q allele of a genetic marker (as outline herein) will have an approximately 78% probability of having a larger size than if it had two copies of a different allele of that genetic marker.
SHR507647NZPR 301525122 The terms "low(er)", "high(er)", "small(er)", and "large(er)" when used herein in relation to size, including weight, height, and/or growth rate, are intended to mean that an animal will have a size that is "low(er)", "high(er)", "small(er)", "large(er)" than if it had a different allele of a genetic marker (or a different level of expression of a specified 5 protein). In one embodiment, the terms indicate that the animal may have a size that is "low(er)", "high(er)", "small(er)", "large(er)" than another animal or the average in a population. In one particular embodiment, the following values are used as reference points: a height average at 24 months of approximately 1.25 metres with a standard deviation (SD) of approximately 3.7cm; weight average at 24 months of approximately 10 437 kg with a SD of approximately 40kg; average growth rate from birth to 24 months of approximately 0.58 kg/day with a SD of approximately 0.05 kg/day; wherein the population is a population of Holstein Fresian/Jersey crossbreed cattle.
The term "animal" is used herein primarily in reference to mammals within the Bovidae 15 family. In particular embodiments, the animal is a bovine animal. More particularly the animal is Bos taurus or Bos indicus. In one particular embodiment the animal is a beef or dairy breed. By way of further example, the animal may be chosen from the group of animals including, but not limited to, Jersey, Holstein, Friesian, Ayrshire, crossbred dairy cattle, Angus, Hereford, Simmental and crossbred beef cattle.
As used herein the "worth" of an animal refers to an index used to evaluate the value of an animal, for breeding purposes or herd management, for example. The "worth" is the sum of the estimated value of one or more characteristics which may be associated with the animal, typically weighted by an economic value. Exemplary characteristics include milk 25 fat, protein, milk volume, liveweight, fertility, and milk somatic cells. The term "worth" should be taken to encompass "breeding worth" and other known indexes used to assess the value of an animal. Persons skilled in the art to which the invention relates will readily appreciate methods and formulae suitable for estimating breeding worth on the basis of any number of different characteristics. Further details of parameters and methods for 30 calculating the worth of an animal are provided in the "Examples" section herein after.
In certain embodiments, the methods of the invention allow for "selection" or "rejection" of an animal based on the genetic marker present or the level of a specified protein (including precursors, fragments and/or isoforms thereof) and/or nucleic acids encoding SHR507647NZPR 301525122 said protein, precursor, fragment and/or isoform. Typically an animal is selected for purposes relating to the beef industry or rejected for purposes relating to the dairy industry where the genetic marker or level of protein (including precursors, fragments and/or isoforms thereof) and/or nucleic acid encoding same infers a larger size potential for the 5 animal. Similarly, an animal is selected for purposes relating to the dairy industry or rejected for purposes relating to the beef industry where the genetic marker or level of protein (including precursors, fragments and/or isoforms thereof) and/or nucleic acid encoding same infers a smaller size potential for the animal.
The methods of the invention relate generally to inferring the size potential of an animal and/or its offspring. In one particular embodiment, the methods relate to inferring the liveweight potential of an animal and/or its offspring. In another particular embodiment, the methods relate to inferring the growth rate of an animal and/or its offspring. As mentioned hereinbefore, such information can be used, for example, to select animals for 15 breeding purposes and/or estimate an animal's worth. The invention may also be used to predict bull or cow phenotype performance and as such can be used in production management systems known as Marker Assisted Selection. Animals more or less suitable for a particular production system can be screened early at birth or as embryo's to predict life time performance and segregated or managed to suit their genotype and therefore 20 predicted phenotype.
In one embodiment of the invention, the methods involve the analysis of a nucleic acid from the animal to determine the nucleotide sequence of at least one genetic marker as described herein, wherein the sequence infers the size potential of the animal and/or its 25 offspring. In another embodiment, the methods involve the analysis of the level of a specified protein or nucleic acid encoding it, wherein the level infers the size potential of the animal and/or its offspring.
Typically, the methods involve taking a sample from an animal to be tested. The sample 30 may be any appropriate tissue or body fluid sample. In one embodiment, the sample is blood, muscle, bone, somatic cell(s), saliva, or semen. In other embodiments, the sample is taken from the liver or brain. Such samples can be taken from the animal using standard techniques known in the art. It should be appreciated that a sample may be taken from an animal at any stage of life, including prior to birth; for example, an embryo, feotus, SHR507647NZPR 301525122 blastocyst. Individual gametes could also be tested using the methods of the invention.
This may assist in breeding and/or cloning programmes. A sample may also be taken after the death of an animal. The samples are analysed using techniques which allow for the observation or analysis of levels of nucleic acid and/or the sequence of a particular nucleic 5 acid, as will be described further herein after.
The methods of the invention may be combined with one or more other methods of use in assessing genotype, predicting phenotype, selecting an animal based on certain characteristics, estimating breeding values or estimating worth and the like. Accordingly, 10 the methods of the invention may include, in addition to analysis of the markers identified herein, analysis of additional genetic markers, and/or the level of expression of certain genes/proteins, and/or one or more phenotypic traits, for example.
Genetic Markers In one embodiment, methods of the invention involve analysing one or more genetic markers to identify the sequence or genotype of the genetic marker(s).
Genetic markers of use in the invention are listed in Table 1 below, along with the genotypes which the inventors have identified to be indicative of a smaller size (q allele), 20 or a larger size (Q allele). The sequence and position provided for each genetic marker is based on the genomic sequence of chromosome 14 in bovine build Btau 4.0 (gi|158137890|gb|CM000190.3| in the GenBank database http://www.ncbi.nlm.nih.gov/) The position of each of the genetic markers should be read in accordance with base position being the start site of the polymorphism, given that the first nucleotide in the 25 sequence (gi| 158137890|gb|CM000190.31) is denoted as position one. Further sequence information in the region of the genetic markers of Table 1 is provided in Table 10 herein after.
SHR507647NZPR 301525122 Table 1 Marker Chromosome Position Q All q All FJX 21374 1:1 14 22986260 T C FJX 221494 14 23186380 T C FJX 221762 14 23186648 A G FJX 222352 1:1 14 23187238 A G FJX_240394 14 23205280 A G FJX_247103 14 23211989 - G FJX 250879 1:1 14 23215765 C ' G FJX 253666 14 23218552 T C FJX 263178 1:1 14 23228064 G T FJX 267438 2:2 14 23232324 A G FJX 270353 1:1 14 23235239 C G FJX_278702_1:1 14 23243588 «.
TT FJX-278848 14 23243734 A G 23264810 (CCG)11 (CCG)9 23264854 G A 23272034 A G 23273471 T C FJX_PLAPROTRI_l:l 14 FJX_PLAPROSNP_l:l 14 FJX_307148_1:1 14 FJX_308585_1:1 14 Whilst the nucleotide variations listed in Table 1 provide specific examples of markers of use in the methods of the invention, the inventors contemplate that any variation in a 5 nucleotide or the nucleotide sequence at any one or more the genetic loci referred to may be indicative of size potential.
In one particular embodiment, the marker(s) are chosen from the group outlined in Table 8 herein after. In another particular embodiment, the marker(s) are chosen from 10 FJXPLAPROTRIl: 1 and FJX_PLAPROSNP_l: 1.
It should also be appreciated that one could analyse the nucleic acid sequence of either strand of the nucleic acid to identify the sequence at a particular genetic loci or position; ie instead of analysing the strand associated with the sequence variants listed in Table 1, the 15 nucleotide sequence of the opposite or complementary strand of DNA could be analysed. Persons of skill in the art will readily appreciate nucleic acid sequence variations on such SHR507647NZPR 301525122 opposite strand which correlate with the genotypes of Table 1, having regard to the information contained herein and nucleic acid base pairing principles (ie, A pairs with T and C pairs with G).
It should be appreciated that the invention also encompasses use of one or more genetic markers which are in linkage disequilibrium with one or more of the markers mentioned herein before. Such markers may be analysed instead of or in addition to the genetic markers mentioned herein before.
"Linkage disequilibrium" should be taken broadly to refer to the tendency of the presence of an allele at one genetic loci (for example a genetic marker of Table 1) to predict the presence of an allele at one or more other genetic loci (for example a distinct genetic marker). The genetic loci need not necessarily be on the same chromosome. However, in a preferred embodiment, the genetic loci are located on the same chromosome.
One measure of linkage disequilibrium is DELTA2, which is calculated using the formula described by Devlin et al (Genomics 29 (2):311-22 (1995)), and is a measure of how well an allele X at a first genetic loci predicts the occurrence of an allele Y at a second genetic loci. A DELTA2 value of 1.0 indicates the prediction is perfect (for example, if Y is present then X is present). It should be appreciated that reference to linkage disequilibrium herein should not be taken to imply a DELTA2 value of 1.0. In particular embodiments, the linkage disequilibrium between an allele at one genetic locus and an allele at a second genetic locus, has a DELTA2 value of at least 0.75, at least 0.80, at least 0.85, at least 0.90, at least 0.95, and most preferably 1.0.
Skilled persons will readily appreciate methods for determining whether any two alleles are in linkage disequilibrium. However, by way of example, see Genetic Data Analysis II, Weir, Sinauer Associates, Inc. Publishers, Sunderland, Mass., 1996.
Nucleic acids can be analysed to determine the genotype/sequence of the genetic markers described herein according to any appropriate technique. Such techniques include for example polymerase chain reaction (PCR), including allele-specific PCR, gel electrophoresis, the use of oligonucleotide probe hybridisation, Southern blotting, direct sequencing, restriction digestion, restriction fragment length polymorphism (RFLP), SHR507647NZPR 301525122 single-strand confirmation polymorphism (SSCP), LCR (ligase chain reaction), denaturing gradient gel electrophoresis (DGGE), the use of allele-specific oligonucleotides (ASOs), the use of proteins which recognize nucleic acid mismatches, such as E.coli mutS protein, RNAse protection assays, oligonucleotide array hybridisation (for example microarray), 5 denaturing HPLC (dHPLC), fluorescence quenching PCR (TaqManâ„¢ Applied Biosystems, CA 94404, USA), High Resolution Melting (HRM), and matrix-assisted laser desorption/ionisation time-of-flight mass spectroscopy (MALDI-TOF MS). Combinations of two or more of such techniques may be used. Such combination may increase the sensitivity of the analysis being conducted.
The technique(s) used will depend on the nature of the marker to be detected as will be appreciated by skilled persons. For example, single nucleotide polymorphisms (SNPs), may be analysed using those techniques capable of resolving a single nucleotide difference between sequences; for example, direct sequencing or LCR, allele-specific PCR, RFLP, 15 SSCP, DGGE, using allele-specific oligonucleotides (ASOs), or proteins which recognize nucleic acid mismatches, oligonucleotide array hybridisation, dHPLC, fluorescence quenching PCR and matrix MALDI-TOF MS.
Any one or more of the techniques mentioned hereinbefore (including for example, SSCP, 20 RFLP, DGGE, dHPLC and direct sequencing) may be used to analyse the genetic markers which may include insertion or deletion of one or more nucleotide.
It should be appreciated that certain of the techniques of use in analysing the genetic markers of the invention will utilise one or more oligonucleotides which hybridise to a 25 genetic region encompassing the marker, adjacent to the marker, or flanking the marker. Such oligonucleotides may be DNA, RNA or derivatised forms thereof and include nucleic acid primers, such as PCR and LCR primers, and nucleic acid probes.
Persons of ordinary skill in the art to which the invention relates will readily appreciate 30 appropriate oligonucleotides of use in the invention having regard to the nucleic acid sequence of chromosome 14 (as detailed herein before), particularly in the genetic regions proximal to a particular marker, the nature of the genetic markers to be analysed, and the general principles of nucleic acid hybridisation. The nucleic acids will be capable of hybridising in a specific manner to a target nucleic acid and in the case of primers they will SHR507647NZPR 301525122 be capable of priming a PCR or like reaction. While such nucleic acids will preferably have 100% complementarity to their target region of the mRNA or cDNA of the protein of interest, they may contain one or more non-complementary nucleotides at a particular position while still substantially retaining specificity for the target nucleic acid to which 5 they are designed to bind. By way of example, the nucleic acids may have approximately 80%, approximately 90%, approximately 95%, or approximately 99% complementarity or homology to its target. By way of further example, in certain cases, the oligonucleotides may be designed such that a mismatch at a particular nucleotide position is indicative of the nature of the genetic marker being analysed (for example, a SNP). By way of example, 10 a mismatch in the nucleotide present at the 3' end of an LCR primer will inhibit the reaction providing an indication of the nature of the nucleotide at that position. Mismatches may similarly be utilised in techniques including RNAse protection assays and allele-specific PCR, as well as in fluorescence quenching PCR, for example. Typically, the nucleic acids will hybridise to their target nucleic acid under stringent hybridisation 15 conditions (see for example, Sambrook and Russell, Molecular Cloning: A Laboratory Manual, 2001, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York).
The oligonucleotide probes or primers may be of any length as is appropriate for a particular application, having regard to the sequence of the genetic region to which they 20 are designed to bind. A probe or primer will typically be capable of forming a stable hybrid with the complementary sequence to which it is designed to hybridise. Accordingly, the length is dependent on the nucleic acid composition and percent homology between the oligonucleotide and its complementary sequence, as well as the hybridisation conditions which are utilised (for example, temperature and salt 25 concentrations). Such hybridisation factors are well known in the art to which the invention relates. By way of example, oligonucleotides of use in the present invention may be from 2 to 500 nucleotides in length. In one embodiment, particularly where they are used as primers, the oligonucleotides may be of approximately 15 nucleotides to 30 nucleotides in length.
Oligonucleotide probes and primers of use in the invention may be prepared by any number of conventional DNA synthesis methods including recombinant techniques and chemical synthesis, or they may be purchase commercially. It will be appreciated that the usefulness of any probe or primer may be evaluated, at least notionally, using appropriate SHR507647NZPR 301525122 software and sequence information for the nucleic acid encoding the protein of interest. For example, software packages such as Primer3 (http://primer3.sourceforge.net/), PC 01igo5 (National Bioscience Inc), Amplify (University of Wisconsin), and the PrimerSelect program (DNAStar Inc) may be used to design and evaluate primers.
Where amplification techniques (for example PCR) are used in methods of the invention amplification may be conducted according to conventional procedures in the art to which this invention relates, such as described in US Patent No 4,683,202. By way of example PCR reactions will generally include 0.1pM-l|iM of each primer, 200(oM each dNTP, 3-10 7mM MgCb, and 1U Taq DNA polymerase. Further, exemplary PCR cycling conditions include: denaturation at a temperature of approximately 94°C for 30 to 60 seconds, annealing at a temperature calculated on the basis of the sequence and length of the primer (as herein after discussed) for 30 to 60 seconds, and extension at a temperature of approximately 70°C to 72°C for 30 to 60 seconds. By way of example, between 25 and 45 15 cycles are run.
It will be appreciated by those of ordinary skill in the art that the amplification conditions provided herein are merely exemplary and may be varied so as to optimise conditions where, for example, alternative PCR cyclers or DNA polymerases are used, where the 20 quality of the template DNA differs, or where variations of the primers not specifically exemplified herein are used, without departing from the scope of the present invention. The PCR conditions may be altered or optimised by changing the concentration of the various constituents within the reaction and/or changing the constituents of the reaction, altering the number of amplification cycles, the denaturation, annealing or extension times 25 or temperatures, or the quantity of template DNA, for example. Those of skill in the art will appreciate there are a number of other ways in which PCR conditions may be optimised to overcome variability between reactions.
It will be understood that whilst not specifically exemplified herein, appropriate annealing 30 temperatures for any primer within the scope of the present invention may be derived from the calculated melting temperature of that primer. Such melting temperatures may be calculated using standard formulas, such as that described in Sambrook and Russell, 2001. As will be understood by those of ordinary skill in the art to which this invention relates annealing temperatures may be above or below the melting temperature but generally an SHR507647NZPR 301525122 annealing temperature of approximately 5°C below the calculated melting temperature of the primer is suitable.
Oligonucleotides used for detection and/or analysis of the genetic markers of the invention 5 may be modified to facilitate such detection. Similarly, nucleic acid products obtained using techniques such as PCR may be modified to facilitate detection and/or analysis. For example, the nucleic acid molecules may be labelled to facilitate visual identification using techniques standard in the art. By way of example nucleic acids may be radio-labelled using P32 as may be described in Sambrook and Russell, 2001. Further, nucleic acids may 10 be appropriately labelled for use in colorigenic, fluorogenic or chemiluminescence procedures. Specific examples are provided herein after, in the "Examples" section.
It will be appreciated that the methods of this embodiment of the invention may employ one or more control samples. Such control samples may be positive or negative controls 15 for a particular genetic marker. The type of control samples used may vary depending on such factors as the nature of the genetic marker being analysed and the specific technique being used for such detection and analysis. Positive controls may include samples having known nucleic acid sequences. Negative controls may include samples having no nucleic acid present. By way of general example, in analysing a SNP positive control samples 20 could include nucleic acids known to have a particular nucleotide at the relevant position.
In order to facilitate detection of a genetic marker in accordance with the invention, a sample may be processed prior to analysis. For example, the sample may be processed to isolate nucleic acid from the sample to be analysed or to amplify a specific genetic region 25 to be analysed.
In one embodiment, nucleic acid is isolated or extracted from the sample prior to analysis. In one embodiment, genomic DNA is isolated or extracted from the sample. In an alternative embodiment, mRNA may be isolated or extracted from the sample. In such a 30 case, the mRNA may be converted to cDNA using reverse transcription techniques known in the art. Techniques for isolating nucleic acids from samples will be readily appreciated by skilled persons. By way of Example, methods of use in isolating nucleic acids are described in Sambrook and Russell, 2001.
SHR507647NZPR 301525122 In an alternative form of this embodiment of the invention analysis of the nucleic acid may occur in situ obviating the need to extract nucleic acid from the sample. This may be done using PCR for example. Skilled persons will readily appreciate appropriate techniques and methodology to this end (see for example, Sambrook and Russell, 2001).
Expression levels In another embodiment, the methods of the invention involve observing the level of one or more of the proteins RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK, including any one or more isoforms, precursors and fragments thereof, and/or any nucleic 10 acids encoding one or more of the foregoing.
The proteins RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK are described further in the "Examples" section herein after. Exemplary amino acid and gene sequences for each protein are described in http ://www.ensembl.org/. By way of example 15 see: RPS20 - ENSBTAT00000025484; MOS - ENSBTAT00000025483; PLAGl -ENSBTAT00000005251; CHCHD7 - ENSBTAT00000047298; RDHE2 -ENSBTAT00000052250; SDR16C6 - ENSBTAT00000052250; and, PENK -ENSBTAT00000006478.
The methods of this embodiment of the invention will typically involve taking a sample from an animal, observing the level of expression of one or more of the proteins (including isoforms, precursors or fragments thereof) or nucleic acids encoding same, and comparing the level detected against one or more standard. Any difference in the level of expression between the sample and the standard infers the size potential of an animal and/or its 25 offspring, including growth rate potential and/or liveweight potential.
In one particular embodiment, the sample is a muscle, bone, brain or liver sample.
In one embodiment, at least an approximately 1.2-fold difference in the level of one or 30 more of the proteins, a precursor, fragment, or isoform thereof and/or nucleic acid encoding them compared to the standard is indicative of size potential. In one embodiment, there is at least an approximately 1.2-fold increase in the level of one or more of the proteins, a precursor, fragment, or isoform thereof and/or nucleic acid encoding any one thereof.
SHR507647NZPR 301525122 In accordance with this embodiment of the invention, a standard is a control sample having a known level of one or more of the proteins (including isoforms, precursors and fragments thereof) or a nucleic acid encoding same, which is tested concurrently with the sample 5 from an animal to be tested. However, in another embodiment, the standard could be a printed chart or electronic information or the like containing previously generated data considered to provide an appropriate standard and which the test samples could be compared to on the basis of colour, fluorescence levels, or numerical values, for example.
In accordance with this embodiment of the invention, the standard is preferably a level of one or more of the proteins (including isoforms, precursors or fragments thereof) or a nucleic acid encoding same, which is associated with an expression level in an animal or animals having a particular size (including growth rate and liveweight).
As mentioned herein before, the inventors have identified that an increase in expression of one or more of the above proteins correlates with a larger size. Whilst any increase in the level of expression one or more of these proteins (including isoforms, precursors and fragments thereof or nucleic acids encoding same) may be considered to infer size potential of an animal, in a one embodiment, an increase in the level of expression of at 20 least approximately 1.2 fold infers an increase in the size potential of an animal and/or its offspring.
A decrease or substantially no increase in the level of expression compared to a standard may also infer a lower size (including growth rate and/or liveweight) potential for the 25 animal and/or its offspring.
The one or more proteins (including precursors, fragments and isoforms thereof) and nucleic acids encoding same may be detected and the levels thereof compared to a standard using any one or a combination of techniques which are of use in identifying, quantifying 30 and/or highlighting differential levels or expression of one or more proteins. Such techniques will be readily appreciated by persons of ordinary skill in the art to which the invention relates. However, by way of example protein purification methods, immunological techniques, separation of proteins based on characteristics such as molecular weight and isoelectric point including gel electrophoresis and microfluidics- SHR507647NZPR 301525122 based technologies as for example in gel-free protein separation techniques, and mass spectroscopy (MS) utilizing isobaric label based MS such as iTRAQ or label-free approaches such as multiple reaction monitoring (MRM) may be employed.
Appropriate immunological techniques include enzyme linked immunosorbent assay (ELISA) (sandwich ELISA, double sandwich ELISA, direct ELISA, microparticle ELISA), radioimmunoassay (RIA), immunoprecipitation, Western blotting, immunohistochemical staining, antibody arrays, or agglutination assays. Protocols for carrying out such techniques are readily available; for example, see "Antibodies a 10 Laboratory Manual", Cold Spring Harbor Laboratory Press (1988).
Antibodies of use in such immunological techniques may be purchased commercially or produced according to standard methodology in the art having regard to the nature of the proteins to be tested. For example, polyclonal antibodies and monoclonal antibodies may 15 be produced in accordance with the procedures described in the text "Antibodies a Laboratory Manual" (Cold Spring Harbor Laboratory Press, 1988) using one or more of the proteins or a fragment thereof as antigen. Preferably monoclonal antibodies are used.
Nucleic acid-based techniques of use in the invention include differential display 20 procedures, Northern Blotting, and competitive PCR. Persons skilled in the art to which the invention relates will readily appreciate methodology for performing these techniques.
Nucleic acids, such as oligonucleotide probes and primers, of use in detecting expression levels of proteins in accordance with the invention (for example using Northern blotting or 25 competitive PCR) will be readily appreciated by skilled persons having regard to the information contained herein and any published amino acid and/or nucleic acid sequence information for RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK. The nucleic acids will be capable of hybridising in a specific manner to a mRNA or cDNA associated with RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK and in 30 the case of primers they will be capable of priming a PCR or like reaction.
Mass spectroscopy techniques of use in the invention are described for example in "Proteins and proteomics-A laboratory manual" (RJ Simpson, Cold Spring Harbour Laboratory Press (2002).
SHR507647NZPR 301525122 The difference in the levels RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK or nucleic acids encoding RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK in a sample versus a standard may be compared using standard technology 5 having regard to the method employed to detect the protein or nucleic acid. For example, colorimetric and fluorometric techniques may be used in which a detection molecule (such as an antibody or nucleic acid probe or primer) is labelled with a molecule which can be visualised by the naked eye or otherwise detected using a spectrophotometer, or fluorometer for example. Alternatively, detection molecules could be labelled with radio-10 isotopes. Incorporating labels into nucleic acids during PCR amplification where it is employed (as opposed to labelling a detection molecule such as a probe or primer), is also contemplated.
Methods for labelling molecules and subsequently measuring the intensity of signals 15 generated will be known to those of skill in the art to which the invention relates.
It should be appreciated that in addition to analysing samples and standards, the methods of the invention may include the testing of one or more positive or negative control samples to ensure the integrity of the results. For example, one could include a sample containing no protein/nucleic acid and one or more samples containing a known level of 20 protein/nucleic acid so that results can be calibrated across different runs of the method.
The sample may be processed prior to analysing one or more of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK and/or nucleic acids encoding same to facilitate analysis of the proteins or nucleic acids. Skilled persons will readily appreciate 25 appropriate processing steps and techniques suitable for performing them.
In one embodiment, high abundance proteins which have the potential to make it difficult to analyse, such as detect and/or measure the level of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and/or PENK, may be removed from the sample. For example, Top6 30 or Top7 depletion may be used (as described herein after under the heading "Multi affinity removal system (MARS) for major abundance proteins"). The sample may also be subject to proteolytic digestion. As such detection of a protein or isoform in accordance with the invention should be taken to include detection of any one or more fragments thereof. Fragments should be of a length sufficient to ensure specificity to one or more of RPS20, SHR507647NZPR 301525122 MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK. Such fragments will for example be at least 8 amino acids in length, more preferably at least 10,15 or 20 amino acids in length.
Processing steps for preparing the sample for analysis of nucleic acids encoding one or more of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK, may include lysing cells, isolating mRNA, and generating cDNA using standard procedures such as reverse transcription-PCR as will be known in the art to which the invention relates. In one embodiment, mRNA may be observed in situ.
Skilled persons may readily appreciate other means by which the sample may be processed for use in the invention.
Kits The invention also relates to kits which are of use in a method of the invention.
In one embodiment, the kit comprises at least one or more reagents suitable for analysis of the sequence of one or more of the genetic markers referred to herein. Reagents suitable for analysis of one or more of the variants include one or more nucleic acid probes and/or 20 primers as herein before described.
In another embodiment, the kit comprises at least one or more reagents suitable for detection of one or more of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK (including precursors, isoforms and fragments thereof), and/or nucleic acids 25 encoding same.
By way of example, where an immunological procedure is used the kit may comprise one or more antibody specific to one or more of the proteins (including precursors, isoforms and fragments thereof). In a particular embodiment, ELISA is used and the kit comprises 30 one or more capture and/or detection antibody for one or more of the proteins (including precursors, isoforms and fragments thereof).
By way of further example, where a method of the invention involves detection of the level of nucleic acids encoding one or more of the proteins (including precursors, isoforms and SHR507647NZPR 301525122 fragments thereof), it may comprise one or more nucleic acid probes and/or primers which have specificity for the target nucleic acids.
In another embodiment, the kit comprises at least one reagent suitable for analysis of the 5 sequence of the one or more genetic markers and at least one reagent suitable for detection of one or more of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK (including precursors, isoforms and fragments thereof), and/or nucleic acids encoding same.
Reagents of use in processing samples for analysis may also be contained in the kits of the invention. The kits may also comprise one or more standard and/or other controls including nucleic acids whose sequence or genotype at a particular position is known, or containing known levels of one or more of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK (including precursors, isoforms and fragments thereof) or nucleic 15 acids encoding same. Further, kits of the invention can also comprise instructions for the use the components of the kit as well as printed charts or the like that could be used as standards against which results obtained from test samples could be compared. Reagents may be held in any suitable container.
Breeding The invention also relates to methods of breeding animals to produce offspring.
In one embodiment the method involves selecting at least a first animal that has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 25 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11,23264854G, 23272034A, and 2327347IT and mating said first animal with a second animal to produce offspring, wherein the offspring has a smaller size potential. In one embodiment, the second animal is selected on the basis that it also has one or more of the following 30 genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11,23264854G, 23272034A, and 23273471T.
Received at IPONZ on 24.02.2012 SHR507647NZPR 302016060 39 In another embodiment, the method involves selecting at least a first animal that has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 5 23272034A, and 23273471T, and mating said first animal with a second animal to produce offspring, wherein the offspring has a larger size potential. In one embodiment, the second animal is selected on the basis that it also has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 10 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T.
In one embodiment, the method involves observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, 15 and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in at least a first animal, selecting the first animal where there is an increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, and mating said first animal with a second animal to produce offspring, wherein the one or more proteins are selected from 20 RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK. In one embodiment, the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in the second animal, and selecting the second animal where there is an increase 25 in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK. Preferably the first and/or second animals are selected for breeding for beef farming where the level of expression is increased compared to a standard.
In another embodiment the method involves observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or Received at IPONZ on 24.02.2012 SHR507647NZPR 302016060 40 more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in at least a first animal, selecting the first animal where there is substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to 5 a standard, and mating said first animal with a second animal to produce offspring, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK. In one embodiment, the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids 10 encoding the one or more proteins, precursors, isoforms or fragments in the second animal, and selecting the second animal where there is substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK. Preferably the first 15 and/or second animals are selected for breeding for dairy farming where there is no substantial increase in the level of expression compared to a standard.
In one embodiment, an animal is rejected for breeding for the purposes of dairy farming where there is an increase in the level of expression compared to the standard. In one 20 embodiment, the animal is rejected for breeding for the purposes of beef farming wherein there is substantially no increase in the level of expression compared to the standard.
It should be appreciated that any appropriate breeding methods may be utilised in this embodiment of the invention including natural insemination, artificial insemination and in 25 vitro fertilisation. Accordingly, the word "mating" should be construed broadly and not limited to the physical pairing of two animals.
Preferably, the methods herein before described may be used to determine the genotype of one or both potential parents and/or the level of expression of one or more specified 30 protein.
The invention is now further elucidated by the following non-limiting examples.
SHR507647NZPR 302016060 Received at IPONZ on 24.02.2012 40A EXAMPLES Example 1 A QTL with major effect on body size maps to bovine chromosome 14, To identify QTL affecting traits of importance for the dairy sector, a Holstein-Friesian (HF) x Jersey (J) line-cross comprising 864 F2 cows (Fig. 1) was generated. To that end, 36 HF bulls were mated with 430 J dams, and 24 J bulls with 366 HF dams, yielding six F1 sires (three HFxJ and three JxHF) and 796 F1 dams (430 HFxJ and 366 JxHF). F1 animals were mated in order to obtain approximately equal proportions of the four possible 10 genotypes ((HfxJ)x(HFxJ); (HfxJ)x(JxHF); (JxHF)x(HFxJ); (JxHF)x(JxHF)). This experimental design primarily targets QTL segregating between HF and J, but - if analyzed appropriately - also allows the identification of QTL segregating within breeds.
SHR507647NZPR 301525122 More than 500-hundred traits were measured on the F2 animals (Spelman et al, 2007). More than twenty of these pertained to body size, including weight at birth, 6, 8, 12, 18 and > 24 (-'live weight") months, as well as height at withers at 18 months. All animals 5 of the pedigree were genotyped for 294 microsatellite markers spanning the bovine genome. When applying a line-cross model (i.e. assuming fixation of alternate QTL alleles in HF and J)(Haley et al., 1994) to weight and height, the genome-wide location scores shown in Fig. 2A were obtained. The strongest signal, exceeding genome-wide significance of P-value <0.00001, was on the proximal half of chromosome BTA14.
To increase the information content of BTA14, increased the microsatellite density was increased from eight to 56. The data was re-analyzed using the linear regression option of the HSQM software (Coppieters et al., 1998), treating the data as six independent paternal half-sib pedigrees. Across-family analysis yielded the location scores shown in Fig. 2B, 15 confirming the presence of a genome-wide significant QTL affecting body size on proximal BTA14. The effect on weight was clearly detectable at birth, but its significance increased with age. It was accompanied by a highly significant, co-localized effect on height. Bootstrap analysis (Visscher et al., 1996) of live weight (the most significant trait) defined a 18.4 Mb confidence interval (CI) spanning positions 9.0 Mb to 27.4 Mb on 20 BTA14. Within-family analyses yielded significant evidence for QTL segregation in four of the six sire families, suggesting (under a biallelic QTL model) that sires 1, 2, 3 and 5 were of "Qq" QTL genotype, while the two remaining sires (4 and 6) were either of "QQ" or "qq" genotype (Fig. 2C). QTL position and allele substitution effects were consistent across the four "Qq" sires (Fig. 2D): ~20 Kgs for live weight and ~2 cms for height, 25 corresponding to predicted differences between full-grown QQ and qq animals of -40 Kgs and ~4 cms.
Analysis of the scientific literature indicates that the same QTL is likely to be segregating in other populations, thus increasing the interest of identifying the corresponding causative 30 mutation(s). Mizoshita et al. (2004) and Takasuga et al. (2007) reported QTL affecting weight of Japanese Black Cattle on BTA14. Buchanan et al. (2005) reported association between polymorphisms in the CRH gene (located on chromosome 14) and growth in Charolais-cross cattle. Kneeland et al. (2004) reported QTL influencing birth weight as well as average daily gain pre- and post-weaning in commercial beef cattle. Nkrumah et SHR507647NZPR 301525122 al. (2007) reported QTL affecting average daily gain in offspring of Angus, Charolais and Alberta Hybrid bulls. Using haplotype-sharing methods, Mizoshita et al. (2005) fine-mapped one of their growth QTL to a 1.1 Mb interval (flanked by microsatellites BMS1914 and INRA094) encompassed by the CI defined in this study (Fig. 2B).
Combined linkage and linkage disequilibrium (LD) analysis in the HFxJ intercross and outbred dairy cattle population positions the QTL within a ~750 Kb chromosome segment To refine the map position of the QTL the BTA14 marker density was increased focusing 10 on the region spanning the QTL peak. The entire HFxJ cross was genotyped for 925 SNP markers corresponding to (i) 629 public domain SNPs, and (ii) 296 novel SNPs. The latter were identified by resequencing coding exons as well as highly conserved sequence elements in the six F1 sires. Combined linkage and linkage disequilibrium analyses (e.g. Blott et al., 2003; Druet & Georges, 2009) of live weight (i.e. the most significant trait) 15 were performed. Highly significant location scores were obtained, maximizing at position 23,716,140 (single marker analysis; lodscore: 31.4; Fig. 3A) or 23,221,130 (haplotype-based analyses; lodscore: 34.1; Fig. 3B), i.e. ~1 Mb proximal from the 1.1 Mb interval defined by Mizoshita et al. (2005). Fitting the most likely QTL position in the model and rescanning the chromosome for residual QTL effects didn't reveal any evidence for a 20 second QTL in the vicinity (Fig. 4). It is noteworthy that at this stage not a single SNP nor haplotype matched the QTL genotypes of the six F1 sires (1-2-3-5: Qq; 4-6: QQ/qq) (data not shown).
The effect of chromosome 14 on body size was then evaluated in the New-Zealand outbred dairy cattle population. A dataset corresponding to 2,700 progeny-tested bulls genotyped for the USDA 50K Illumina SNP chip (Van Tassell et al., 2008) was used. Phenotypes were breeding values for live weight. The data were analyzed using the same mixed model as above (extracting both linkage and linkage disequilibrium information) with addition of a regression on percentage Jersey blood given available pedigree information. The obtained location scores are shown in Fig. 3C. A maximum LRT of 156 (corresponding to a lod score of 34.0) was obtained at position 23.425.845, i.e. exactly the same location as in the HFxJ cross despite the use of a distinct SNP panel. These results considerably strengthened the confidence in the location of the QTL. Figure 5 shows the effect on breeding value for live weight of the different haplotype clusters as well as their respective SHR507647NZPR 301525122 frequencies in the New Zealand dairy cattle population. The bimodal distribution of haplotype effects supports a bi-allelic QTL. As expected, the q allele associated with lower weight/height is virtually fixed in the Jersey population, while the Q allele associated with higher weight/height predominates in the Holstein Friesian population. It 5 is noteworthy that, in the outbred population, fitting the most likely QTL position in the model and scanning BTA14 for additional QTL effects did reveal modest evidence (lod score 4.3) for the existence of a second QTL in the interval defined by Mizoshita et al. (2005)(Fig. 4).
From there on, further analyses focused on a ~750 Kb segment spanning all most likely QTL positions obtained in the different analyses (22,967,680 - 23,718,980).
Massive parallel resequencing of the -750Kb critical region identifies a cluster of 16 candidate QTN.
A 103 long-range PCR products spanning the entire 750 Kb interval was generated. PCR products of the predicted size were obtained from HF and J genomic DNA thereby (i) confirming the accuracy of the local sequence assembly, and (ii) indicating that both HF and J alleles could be amplified with comparable efficacy. The same 103 long-range PCR products were then amplified from genomic DNA of the six F1 sires, pooled, size-20 fractionated by nebulization, and appended with adaptors allowing for massive parallel resequencing on a Roche FLX instrument. The adaptors included individual specific multiplex identifiers (MIDs) allowing pooling and simultaneous sequencing of the eight libraries (HF, J and six F1 sires). The obtained sequence traces were analyzed using the GS Reference Mapper software from Roche. The average sequence depth of non repetitive 25 sequences was~20-fold per individual. Mining of DNA sequence polymorphisms using GS Reference Mapper revealed a total of 8,851 putative DNA sequence polymorphisms (DSP) corresponding to an average nucleotide diversity of 0.3%. Assuming that the QTL is indeed bi-allelic (as suggested by the bimodal distribution of haplotype effects), the causative "Quantitative Trait Nucleotide" (hereafter referred to as pQTN for phenotypic 30 QTN) has to be heterozygous in the four F1 sires that segregate for the QTL, homozygous in the two non-segregating F1 sires, and homozygous for alternate alleles in the HF and J samples. Applying this filter to the 8,851 DSP yielded only 17 candidate pQTN: 16 SNPs and one VNTR marker (Table 2). Interestingly, all but one of these clustered in a 87 Kb region spanning MOS, PLAGl and CHCHD7 (Fig. 6).
SHR507647NZPR 301525122 Genotyping assays could be developed for 9 of the 16 SNPs and the VNTR and were used to genotype the HFxJ intercross. All markers in the 87 Kb window proved to be in perfect LD with each other (D'=l; r2 > 0.985), but not with the more distant one (0.859 An orthologous region was shown in two independent genome wide association studies (GWAS) to control height in humans (Gudbjartsson et al., 2008; Lettre et al., 2008). 15 Table 2 Marker Chromosome Position Q All q All Phastcons FJX_21374_1:1 14 22986260 T C 0.625 FJX 221494 14 23186380 T c 0.571 FJX_221762 14 23186648 A G 0.500 FJX_222352_1:1 14 23187238 A G 0.002 FJX 240394 14 23205280 A G 0.002 FJX_247103 14 23211989 - G 0.002 FJX_250879_1:1 14 23215765 C G 0.905 FJX 253666 14 23218552 T C 0.102 FJX_263178_1:1 14 23228064 G T 0.015 FJX 267438 2:2 14 23232324 A G 0.239 FJX 270353 1:1 14 23235239 C G 0.002 FJX 278702 1:1 14 23243588 _ TT 0.004 FJX_278848 14 23243734 A G 0.004 FJX PLAPROTRI 1:1 14 23264810 (CCG) 11 (CCG) 9 0.999 FJX PLAPROSNP 1:1 14 23264854 G ! A 0.996 FJX_307148_1:1 14 23272034 A G 0.033 FJX 308585 1:1 14 23273471 T 1 c NA SHR507647NZPR 301525122 Exploiting haplotype diversity in other breeds excludes 5 of 16 candidate pQTN as causative pQTN.
The 16 QTN being in complete LD within the HFxJ F2 population precluded further genetic differentiation of causative from "passenger" variants in this pedigree material. To 5 overcome this limitation, a panel of 159 animals representing 12 breeds for 12 of the 16 candidate pQTN were genotyped. While the two haplotypes ("g" and "q") observed in the HfxJ cross accounted for 93.8% of the chromosomes in this diversity panel 10 additional "recombinant" haplotypes were observed. One of these ("i?"), carrying the Q allele at 5/12 positions and the q allele at 7/12 positions, reached a frequency of 37% in the Simmental 10 population (Fig. 8A). Therefore DNA samples from 44 unrelated Simmental animals with birth weight information were collected and genotyped for the 12 polymorphisms. The phase of the 44 animals was manually determined. The frequency of the Q, q and R haplotypes in this sample was 22%, 35% and 34%, respectively. The remaining 8% of the chromosomes corresponded to four minor haplotypes. The birth weight of the 44 animals 15 on QTL genotype was regressed assuming (i) that the R haplotype was Q, or (ii) that it was q. While no effect on birth weight was observed under the first hypothesis, a near significant (p=0.06) substitution effect on birth weight with the predicted sign was observed under the second. The difference in birth weight between predicted qq and Qq animals (representing 93% of the sample) was significant (p = 0.02) (Fig. 8B). These 20 results strongly suggest that the 5 Q positions of the R haplotype can be excluded as causative pQTN . However, this does not preclude these 5 Q positions being useful as markers in a method of the present invention.
SHR507647NZPR 301525122 Table 8: Top 8 SNPs Marker Chromosome Position Q All q All Phastcons FJX_21374_1:1 14 22986260 T C 0.625 FJX_221494 14 23186380 T C 0.571 FJX_221762 14 23186648 A G 0.500 FJX_222352_1:1 14 23187238 A G 0.002 FJX_250879_1:1 14 23215765 C G 0.905 FJX_253666 14 23218552 T C 0.102 FJX_263178_1:1 14 23228064 G T 0.015 (CCG) (CCG) FJX_P LAP ROTR I_1:1 14 23264810 11 9 0.999 FJX_PLAPROSNP_1:1 14 23264854 G A 0.996 Unaffected MOS, PLAGl and CHCHD7 open reading frames support regulatory nature ofpQTN.
The candidate pQTN span 87 Kb encompassing three genes (MOS, PLAGl and CHCHD) (Fig. 6A). None of the 16 candidate pQTN resides in the open reading frame (ORF) of either MOS, PLAGl or CHCHD7. Nevertheless, the pQTN could affect the splicing reaction thereby altering coding capacity. To examine this possibility, tissue samples (muscle, bone, brain and liver) were collected from 79 outbred fetuses (Table 9). Fetuses 10 were genotyped for the candidate pQTN, two homozygous "QQ" and two homozygous "qq" fetus were selected, and RT-PCR reactions using overlapping amplicons that would jointly span the complete MOS, PLAGl and CHCHD7 were performed. No evidence for genotype-specific RT-PCR products (Fig. 9) was obtained. These findings indicate that the pQTN affect the expression profile of the causative gene(s) rather than their structure, i.e. 15 they are regulatory pQTN.
SHR507647NZPR 301525122 Table 9: Genotypes of 79 fetuses used for expression analysis N° PLAGl 3'UTR VNTR FJX 250879 1:1 FJX PLAPROTRI 1:1 CHCHD7 SPL.VAR. 3'UTR FJX 303486 1:1 FJX 307148 1:1 14 C/C 341/341 T/T A/A 32 C/G 341/335 T/T A/G G/G 335/335 T/T G/G 9 C/G 341/335 A/T A/G 8 G/G 335/335 A/T G/G 1 G/G 335/335 A/A G/G Alleles associated with increased weight/height are in bold.
The pQTN affect expression of a regulon including RPS20, MOS, PLAGl, CHCHD7, 5 RDHE2, SDR16C6 and PENK.
Being regulatory, the pQTN could not only affect the expression of the spanned MOS, PLAGl and CHCHD7, but also of more distant genes. Thus the effect of pQTN genotype on the expression level of all eight genes in the 750 Kb interval in the fetal tissue samples described above was examined. First expression by QRT-PCR, using up to three 10 amplicons per gene was assayed. Data were normalized using from two to five housekeeping control genes selected from eight candidates using geNorm (Vandesompele et al. 2002). Normalized expression levels were expressed on a log2-scale, relative to the mean of all animals. For a given gene, the individual's average relative expression across amplicons was computed. The average relative expression on the number of Q alleles in 15 pQTN genotype was regressed. In addition to QRTPCR, allelic imbalance assays for PLAGl and CHCHD7 were performed, taking advantage of pQTN located in the 3'UTRs (FJX_250879_1:1 and FLX_307148_1 :1). Comparing the allelic ratio in transcripts of heterozygous individuals is potentially more sensitive (as the comparison of the allelic output is performed within the same sample thereby minimizing the effect of confounding 20 factors), and will also detect differences in promoter strength for genes whose transcription levels are under negative feedback regulation (i.e. rendering steady state transcriptional levels independent of promoter strength). Amplicons encompassing the respective SNPs were amplified from genomic and cDNA and allelic ratios compared by direct sequencing and Peakpicker. (Ge B et al. 2005).
Fig. 10 and Fig. 11 summarize the obtained results. Except for LYN, a significant effect of pQTN genotype on the transcript level of each gene in at least one tissue was observed. The Q allele was associated with increased expression (average 1.2-fold) in 26 of the 29 SHR507647NZPR 301525122 assays performed (where 14.5 are expected if pQTN genotype does not affect expression; p = 2E-12). Significant effects were observed for all four examined tissues. The pQTN effect possibly extends over a longer chromosome segment in brain (from RPS20 to PENK) than in the other tissues (from PLAGl to PENK).
These findings suggest that the pQTN affect the function of (a) long-range cis-acting element(s) that regulate(s) the expression of RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6 and PENK, but apparently not LYN, in multiple tissues. It is noteworthy in this regard that synteny is conserved in zebrafish (chromosome 7) for RPS20, PLAGl, 10 CHCHD7, RDHE2 (=SDRC16C5) and PENK but not for LYN which lies at the proximal end of the cluster in mammals (Fig. 6). The coregulation of multiple genes implies a transcriptional effect. In agreement with this hypothesis, the inventor(s) observed the same degree of pQTN-dependent allelic imbalance when assaying PLAGl pre-mRNAs by means of the FJX_267438_2:2 intronic pQTN candidate (Fig. 10 and Fig. 1 IB).
Reporter and electrophoretic mobility shifty assays support the causality of two pQTN in the PLAG1-CHCHD7 bidirectional promoter Which one(s) of the eight candidate pQTN causes the observed effect on transcription rate? It was noted that two of the eight candidate pQTN (FJX_PLAPROTRI_l:l and 20 FJX_PLAPROSNP_l:l) affect an extremely conserved sequence segment located in the intergenic region separating PLAGl and CHCHD7 (Fig. 5A). FJX_PLAPROTRI_l:l (phastcons score: 0.999) is a (CCG)n trinucleotide repeat with either 9 (q) or 11 (Q) copies located immediately upstream of the presumed human PLAGl transcriptional start site, while FJX_PLAPROSNP_l:l (phastcons score: 0.996) is a A (q) to G (Q) SNP 12-bp 25 upstream from FJX_PLAPROTRI_l:l. Of note, PLAGl and CHCHD7 are positioned head-to-head, separated by only ~500-bp supposed to encompass a bidirectional promoter (Trinklein ND et al. 2004). Of the six remaining candidate pQTN, only FJX_250879_1:1 also affects a conserved element (phastcons score: 0.905) located in the 3'UTR of PLAGl (Table 8).
The effect of the FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l DSP on the presumed bidirectional promoter activity of the PLAG1-CHCHD7 intergenic region was tested. Both allelic forms (Q and q) of a 378-bp and 659-bp fragment centered around the two pQTN were cloned in both orientations in the pGL4 luciferase vector (Fig. 12A).
SHR507647NZPR 301525122 Cosl cells were transfected and measured luciferase activity after 24 hours. When compared to the promoter-less reference vector, both fragments indeed increased luciferase activity: ~9-fold (short) and ~20-fold (long) in the PLAGl direction, and ~90-fold (short) and 44-fold (long) in the CHCHD7 direction (Fig. 12B). Most importantly, the level of 5 luciferase activity was systematically higher with the Q constructs than with the q constructs, the difference being ~1.5-fold, i.e. a magnitude comparable to that observed in vivo (Fig. 12B).
To determine the relative contribution of the FJX_PLAPROTRI_l:l and 10 FJX_PLAPROSNP_l:l variants to the observed effect, recombinant constructs (i.e. being Q for FJX_PLAPROTRI_l:l and q for FJX_PLAPROSNP_l:l or vice versa) were generated, for both the long and short fragments (Fig. 12C). While both variants appeared to be necessary to recapitulate the full eQTL effect in the PLAGl orientation, the FJX_PLAPROSNP_l:l SNP appeared to play a predominant effect in the CHCHD7 15 orientation (Fig. 12D).
To examine whether the effect of the FJX PLAPROTRI l :1 and FJX_PLAPROSNP_l:l variants might reflect differential binding of trans-acting factors, electrophoretic mobility shift assays (EMSA) were performed. First, radiolabeled double-stranded 29-mers 20 centered around the FJX_PLAPROSNP_l:l and nuclear extracts from Cosl, C2C12 and ATDC5 cells, as well as from fetal bone, brain, muscle and liver were used. With all extracts, a number of complexes were observed, some with the same mobility across tissues and some specific, that were systematically 1.2 to 3 x more abundant with the q than with the Q duplex. Moreover, cold q duplex tended to be more proficient than cold Q 25 duplex in displacing either of the radiolabeled probes (Fig. 13). Then 74/80-mer probes encompassing both the FJX_PLAPROTRI_l:l and FJX_PLAPROSNP_l:l variants (Cosl and C2C12 cells) were used. As with the 29-mer oligonucleotide, large complexes that were ~2.5 x more abundant with the q than with the Q probe were observed. As these would be displaced by cold 29-mers (q more effective than Q), they were assumed to be 30 related to the complexes detected with the 29-mer sequence alone. In addition, with C2C12 extracts we detected a jg-specific complex that was not displaced by the 29-mer (Fig. 13).
SHR507647NZPR 301525122 Taken together, these findings suggest that distinct regulatory complexes assemble on the interrogated Q and q sequences, and that this is likely to explain their different promoter strength in reporter assays, and to contribute to the observed differential transcription rate of the Q and q alleles in vivo. Complex assembly in the PLAG1-CHCHD7 intergenic 5 segment is likely to be mediated by broadly expressed and abundant trans-acting factors, as the complexes were readily detected in all examined tissues/cell types.
Exploiting a naturally occurring null allele to exclude the causality of CHCHD7.
Which one(s) of the seven genes affected by the pQTN is (are) controlling stature? RPS20 10 encodes ribosomal protein S20. RPS20 mutations cause p53-mediated skin darkening with pleiotropic effects including reduced body size (McGowan et al. 2008). MOS encodes a protein kinase that is specifically expressed in oocytes where it is involved in the control of meiotic maturation. However, ectopic expression of MOS in somatic cells induces oncogenic transformation (Sagata 1997). PLAGl (Pleomorphic Adenoma Gene 1) was 15 discovered as an oncogene whose ectopic expression due to translocation-induced promoter swapping with ubiquitously expressed genes (such as P-catenin), causes pleomorphic adenomas especially of the salivary gland (Van Dyck et al. 2007). PLAGl encodes a transcription factor that is strongly and broadly expressed during foetal development (particularly in the anterior pituitary), yet severely downregulated after birth. 20 PLAGl has been shown to regulate transcription of several growth factors including IGF2, a key regulator of body size (Van Dyck et al. (2007) Voz et al. (2004)). Interestingly, the most striking phenotype of PLAGl knock-out mice is growth retardation (Hensen et al. 2004). CHCHD7 encodes a widely expressed protein of unknown function, named for the coiled-coil-helix-coiled-coil-helix (CHCH) domain it contains. It was identified as a 25 PLAGl fusion partner in tumors of the salivary gland (Van Dyck et al. 2007). RDHE2 (alias SDR16C5) and SDR16C6 encode members 5 and 6 of the 16C family of short-chain alcohol dehydogenases/reductases. RDHE2 catalyzes the first and rate-limiting step generating retinal from retinol. Imbalances in endogenous retinoids profoundly perturb development and affect growth (Ross et al. 2000). PENK encodes the precursor of met-30 and leu-enkephalins, which play a role in pain perception and response to stress (Kiefer et al. 2002). While the candidacy of PLAGl seem strongest, supported by GWAS signals maximizing on top of the same gene in human, available evidence neither proves PLAGFs causality, nor disproves a contribution of one or more of the other genes.
SHR507647NZPR 301525122 To gain additional insights in the causality of the positional candidate genes, genome-wide expression data generated for 429 of the HFxJ F2 animals was used by hybridizing liver and adipose cDNA (sampled at 60 or 70 months of age) on Affymetrix Bovine 24K expression arrays. The array included probes interrogating LYN, RPS20, MOS, CHCHD7 5 and RDHE2 but not PLAGl and SDR16C6. QTL mapping using HSQM (Coppieters, et al. 1998) and the 56 microsatellite BTA14 map revealed chromosome-wide significant eQTL effects for CHCHD7 in liver and adipose, and for RPS20 in liver. Fine-mapping using the 925 SNP map positioned both eQTL in the resequenced 750-Kb segment, strongly suggesting cis-acting eQTL (Fig. 14 & 15).
Within-family analysis of the CHCHD7 data showed that only the first F1 sire ("1") was segregating for the eQTL, i.e. a segregation pattern in F1 sires distinct from that of the QTL on height/weight. The resequencing data was examined and a SNP matching the eQTL segregation pattern of the six F1 sires that was predicted to affect the donor splice 15 site of CHCHD7 exon 3 (FJX_303486_1:1) was identified. The HFxJ pedigree was genotyped for the corresponding SNP and single-marker association analysis was performed. FJX_303486_1:1 had a highly significant effect on CHCHD7 expression in both liver and fat. The effect was superior to that of all other markers, thus strongly suggesting that FJX_303486_1:1 is the causative "eQTN" (Fig. 14C). Of note, including 20 FJX_303486_1:1 in the association model caused marked drops in the significance of all other markers. Nevertheless, several of them remained significant, including the pQTN (Fig. 14D).
To confirm that FJX_303486_1:1 affects splicing, the 79 fetuses were genotyped: 61 were 25 homozygous TT, 17 were heterozygous AT and one was homozygous AA (without obvious phenotype) (Table 9). RT-PCR experiments were conducted on RNA samples from fetuses of the three genotypes using primers spanning the entire CHCHD7 transcript. RT-PCR products from AA individuals were shortened by 70 bp when compared to TT individuals, which was shown by direct sequencing to result from the skipping of exon 3. The ensuing 30 mutant transcripts have truncated ORF closing within exon four (of five) (Fig. 16A), and are therefore predicted to undergo nonsense mediated RNA decay which may at least in part account for the observed eQTN effect. To confirm the effect of the FJX_303486_1:1 splice variant on CHCHD7 transcript levels, the available QRT-PCR data was exploited and compared CHCHD7 expression levels between qq (pQTN) fetuses that were SHR507647NZPR 301525122 homozygous wild-type versus heterozygous for the FJX_303486_1:1 splice site variant. CHCHD7 transcript levels were indeed found to be reduced ~1.3-fold in AT vs TT animals in all examined tissues, corresponding to mutant transcripts levels ~ lA those of wild-type transcripts (Fig. 16B). Of note, the effect of pQTN genotype on the transcript levels of 5 the eight positional candidate genes is independent of the FJX_303486_1:1 splice site variant (Fig. 10 and Fig. 11).
Assuming that variation in CHCHD7 expression levels influences body size, the allele substitution effects on height/growth should be larger in the offspring of sire 1 10 (heterozygous AT at FJX_303486_1:1) than in any of the other sire-families (sires being homozygous TT). Note that as FJX_303486_1:1 perturbs splicing, it must affect CHCHD7 levels in all tissues. As can be seen from Fig. 2D the QTL allele substitution effects in sire-family 1 did not differ significantly from those in the other sire-families (2, 3 and 5) segregating for the body size QTL. Moreover, and under the same assumption, the 15 FJX_303486_1:1 should have an independent effect on body size. This hypothesis was tested by estimating the effect of the FJX_303486_1:1 SNP on height/weight in the HFxJ line-cross using a mixed model including the pQTN and an individual animal effect. There was no evidence for the slightest residual effect of FJX_303486_1:1 on body size when including the pQTN in the model (data not shown). Finally, assuming that 20 FJX_303486_1:1 creates a CHCHD7 null allele and if CHCHD7 is causally involved in the stature QTL, the [qT/qT] vs [qT/QT] phenotypic contrast should differ from the [qA/qT] vs [qA/QT] contrast. This "quantitative complementation assay" (QCA) (Georges 2007), was performed using a mixed linear model including an individual animal effect. There was no evidence for a significant difference between the two contrasts. Taken together, 25 these results allow formal exclusion of CHCHD7 as being causally involved in the determinism of the stature QTL.
Within family analysis of the RPS20 liver data indicated that sires 1, 5 and 6 were segregating for this eQTL, i.e. a segregation pattern distinct from both the stature QTL and 30 CHCHD7 eQTL (Suppl. Fig. 10). Contrary to CHCHD7 however, the apparent tissue-specificity of the eQTL does not allow exclusion RPS20 on that basis. Analysis of the sequence traces did not reveal obvious candidate eQTN in the immediate vicinity of the RPS20 gene. to Oi Table 10 — Nucleic acid sequence in region of genetic markers SEQID No.
Name Sequence 1 FJX_21374_1:1 Q ACAACTTGAGATCCCTAAACCAGTGTAACCAGGAATGCTTCCCCATAGAACAGTTOAACKjAGACCTTCAAGAAGGTGOAOAAGAAAACATAACAGGAGAG T CTGGCAGAACTGCAAAAGAACTCAGTTCAGTTCAGTTCAGTTCAGTCGCTCAGTCGTGTCCGACTCTTTGTGACCCCATGAATCGCAGCACGCCAGGC C 2 FJX_221494 Q ATCACTGTCTCCTCTTGCTTCTTTCACAGATAAACTCACTTCTTTTATCCTAAAAGTAAAGTCAACCAAACCAGCAATCGCAAAACACTTCTGTGCATCT T ACCCCCTCTTAAGGTTTCTGCTTCTCCCTTCTGCCTCAACAATCCAAGCATAGCTTCTATATCAAGCTGTTGCTTAAAGTGGTGTTGTTAAAGGCAAGT C 3 FJX_221762 Q AGAGACATCACTTTGCCGACAAAGGTCCATATAGTTACAGCTATGGTTTTCCCATTAGTCATGTATGGGTGTAAGAGTTGGACCATAAAGAACACTAAGC A CCAAAGAATTGATGCTTTTGAACTGTGGTGTTGGAGAAGACTCTTGAGAGTCCTTTGGATTGCAAGGAGCTCAAACCAGGCAATCCTAAAGGAAATC AAC 4 FJX_222352_1:1 Q CCTTGATTTCTAGGATGGGAATCTGGTTCTCCTTGTCTGCTTTCCTCCATCCCAATGTTAAATTTTACTATTTCTCAAGGTTATTATTTCTTCTCACATC A TGTAATCTTTCCTTGTTGTCTCATGCAGATAACTGGAAAAGGGACCCAGTTCAGAAGTATTTTAGCAAAATTATACAAGTTCAAAATATCACAAGAAA TC FJX 240394 Q TGGATTTTTACTGATTGTATACTTACAGGAATCCTAAAGGAAATCAATCCTGAATATTCATTGGAAGGACTGATGCTAAAGCTGAAGCTCCAATACTTTG A CCACCTGATGCAAAGAGCTGACTCATTGGAAAAGACCCTGATGCTGGGAAAGATTGAGAGCAGGAGGAGAAGAGGGCAACAGAGGATGAGATGGT TGGAT 6 FJX_247103 Q AGGAGAGCCCAAGGACTCCTTTATAGGTACAATCTCCATTAGGCACTGTGCTGTTTATCTCTTTTCTAAAGCAGTCAGAACTTGGGAACTTTTTTTAACC- 7 FJX_250879_1:1 Q TCAGAGGTAACAGTCTGCCATTGCTAATGAAGAATACATGACGTCTAAAAGATCTCCTCCAATGTCGCCTGGATCTCAGTAGGCTGGAAGCAACTTGGGT C AATATTTTTTTTTGTATTTTAACTTACATACTGGGAAGACAGTGTGTTTTGAGACACGTAAACCTTGTAAAAAAAATAAACCATTGTAAAACCTTATGA A 8 FJX_253666 Q CCAAAAACAGCAGGGAAAAATGAATTAACACTGAATGGAAAAAAAAAAAAAAAACACTGAATGAGAAAGATACCACAGTGCGGGGGGTGGTGGTGGTTGG T GAGCCGGGTTGTTCTGTTTTCTGCTTTTGGTGCTCAGGGAGAGAAAACTGAATCAAAGGGAGCAAGAGTGATGTAAGTGATATGCTGTGACTTGACC CCT 9 FJX_263178_1:1 Q CTTACCTGCATCGATAACTCTTTCATATTTAAGCTAGAAGTCAAGCATATACTTTGCTAATTAAGGGCAAAAGCTATCATAATCAAAGTTTCAGTCTCCTG TATGCTAATCATAAAACTCATTAGTAAACAATTAAGAAAATGTCAACATTTGCTTAAAACATGTTTTTTATGTATAGTAAACAGGTGGGTTGGAAAAG GA FJX_267438_1:1 Q TGCCCCTACAGGCTGCAGACTCTAGCAGCCTGACTGTGTCCTCTTTCCAAGGTCAGCGAACTGTTGCTGTTCCACCAGGATCCTTATTCTAGCATCTGCC A CATCCCACCTGCCAGGGCCAGCCCTCTCCCCACCTGCTCGGAGGTACCCGGACCCTTCCCAGTGGATAGAAAGAGATACAGTGGAAACTCTGGGTAA GAC 11 FJX_270353_1:1 Q TTTAAAAAATTAAGTTTCTTAGCTATTGCTGTTGGGAGTGTAAATTGTTGAAGCTTTTCTGGAATATGTTGTAGAAATAGGTACTCAGAGCCTTTTAAGT C TCCATTTATGCTGATCTGATAAATAATTTCTAAAAATGTAAGGAGAAAAACCAAGAAATAAAGATTTAAATATAAGCATGTTCGAAGCATATTTTAA GGA 12 FJX_278702_1:1 Q ACAGGACAAAATCAAACCCCAAAAAACCCCAAAAAATCAATTCTTCAAGAAAATTTTAAGGGGAAAAAAAAGAAGAGGTAACCTACAGATCAAAATAGAC AGAGATATGTCAACTGACTATATCAGTCAGTCAGTCAGTTCAGTGGCTCAGTCATGTCCAACTCTTTGCAACCCCATGAATCGCAGCACGCCGGGCCT CC l.TmTnT.OV3XXX9WV3XVVXV3XVX3DVVVV3OD0VVXXVVX3DXXX3VXVXV30VV3X0VV0VX30VVXXXVXV3XXX3X3VVXV93XV30X33VXX3 b i:i %L\£9Z Xfd 9Z X33 33VOXX3VOXOX3DXVXVDXOVVX9XVOXOVOVV3DVOODVVV3XVVOX3VVVVOVDVODOV3X30XOOXXXX39X3XXXXOX3XXOXXODD33DVO 3 OOXXOOXOOXOOXOOOOOOOOXOYOVOOVXYOYWOVOXWOXOVOVVVVVVVVVVVVVVWOOXWOXaVDWXXWOXVVVWOOOYaOVOYVVWOD b999£SZ~Xld sz V o XODOXX3VV3OVVODX3DOVXDV3X3XVOOX33D3XOXVV33X33X3XVOVVVVX3XO3VDXV3VXVVOVV9XVVX3DXXV33DX3X0V3VVXOOV0V3X b I:I_6i80£Z~XH n 333X3333X0X3XVXVXVX0VVX0X3003XV3D00VW30V3VX0X30X3XXX0VXX0X0XX0XW3VXX330XXVVX0X9VXX3XXX3XXXXXXXXXXX O33WXXXXXXX3WOOOXX3VVOV3XDV3DVVVX3XXXX3X3XVXXXOX3DXDX3V3OOVXXV33X3XVV3VXOOVXVXXX33X3VOOVV333OVOV0DV b eouw~xrj £Z XY£)OX XOOXVOVOXYOOVOVOWaOOOVOWOVOOYOOVOOVOVOXXVOVWOOOXOOXVOXDDDYOWWOOXXYOXDYOXaOVOVWaOXVOXODVOO o OXXX3VXW33X3DWOX3OWVX3OXV0X3VODVVOOXXV3XXVXVVDX33XVV3XVVVOOVVVX33XVVOOV3VXX3VXVXOXXVOX3VXXXXXVOOX b to£QVZ~Xid zz ox VWOW3V3XVXVVW3IJ£>VV3VXVXXVVVV3DVXXXXVXOWOV3XXOV333VOOOVVVVOOX3VVXVOV30XV3X3XOXXOXX33XXX3XVVXOX O 3XV3V3X3XX3XXXVXLVXXODVV3X3XXXVX3VXXXXVWXXDXW333XV33X33XXX3QX3XDXX3313XXQQX3XVVDDOXVDDVX3XXXVOXX33 b Vl~ZS£ZZZ~Xtd iz 3W 3XWV0DVVVX33XW300V33VW3X30V00VV30XXV00XXX33X0V0V0XX3X3V0VV0VDD1XDX0DX0I3VV0XXXX30XV0XXW0WV33 0 30WX3V3VV0VWXV33V00XX0V0VVX0X090XVX0XV3X0VXXV333XXXX00XVX3DV3VXX0VXVXV33X00VVV3V0330XXX3V3XV3V0V0V b Z9L\ZZXM oz 0 X0W3DQWVXX0XXDXQDXDVVVXX30XX0X3DVV3XVXVX3XX3DVXV3DW33XW3W3X33DX3XX333X3XX30X3XXXD0VVXX3X33333V 3 X3XV30X0X3XX3V3WW303XW30Y33WV33W3X0VWX0WWX33XVXXXX3XX3V3X3VWXV0V3V3XXX3XX30XX3XD3X3X0X3Y3XV b t6PJZZ~Xtd 61 X3 300V33D3V30V3D3XW0XV3333V0X0XXX3X3V033X0X03X0V3X303X0V3XX0V3XX0V3XX0V3XX0V3X3W0VWV30X3W0V3D0X3 3 OY9YOOV3WXV3VVWOWOVf)OXOOWOW3XX33V9VOOWOXXOY3WOVXV3333XX30XYVOOV33WXOXOV33VWX333XVOVOXX3VV3V b VCVL£\Z~XH 81 XV 330XXOOOXOVO0X3V3WOW300V33X3XXVOO0X333XV3303XX300V33V3330V3D0XVOVXV3333V03OX3X3VO33XOXO3XOV3XX30 X XOWX30X3DX30X3XX3XXXOXOX3XV3WXXOOVW33XOWWW3VX3XOXV3XOXOOXXXX3DVXOXXOOXOOXOOXOXXXXV33XOVXOXVOXX I:rS8S80£_Xf;I LI XV X3XXVOXOXVW3XXXXVXXX33W30XVX3VWXXXXXOOVOV3XOOV333WW33XOVX33XOOXX03X30V33XOX3VW3V3VOX30X3003V V OVXXVXOVOOV30XOOX33XOOX3XVOXX3X3XX333XOXOXXXOVXVOV3WWOXOVWOVXVXOXOXV3XXXVW3XVOXXXVXXV3V3VOVWVXO 6 l:I_81?U0£_XfJ 91 O 3333XX333Q3333X3330333D3DX3D330333D33333DX333330333DD330333X33033333033V33D03X3D333DD33333X33D33033 0 3330303X03V033033033033033D33033003033033D330X30V30XX3W3330V30X33330X0000XVX0X0VX0XXV330V3XV33XXXXX Dl :I dNSOtfdVld XfJ £1 33 3X3330333030X303303330333330X33333033300330333X33033333033V33003X303330033333X3303303303330303X03V 11(033)0X30V30XX3W3330V30X33330X303DXVX0X0VX0XXV330V3XV33XXXXX3303VXX33XXX0XXVW00X330003003D3V00X30003XXX 6 I :i raxo^dvid xrj n X3X3WXX0V0XW33XXXX3X0V0V3X330V333X333XW333339X33X33X3XX3333XV3X0X3X33XV3X3XV330V33XV33VXV0X0V3X0V V 3XV33XO3V3X3VO3333V3XX9VOV333X3W33V3XV33XOX333X330O03303V30V303XWOXV3333W30XXX3X3W33XOXV3XOV3X30 b 8fr88£Z~XfJ ei TATGCTAATCATAAAACTCATTAGTAAACAATTAAGAAAATGTCAACATTTGCTTAAAACATGTTTTTTATGTATAGTAAACAGGTGGGTTOGAAAAG GA 27 FJX_267438_1:1 q TGCCCCTACAGGCTGCAGACTCTAGCAGCCTGACTGTGTCCTCTTTCCAAGGTCAGCGAACTGTTGCTGTTCCACCAGGATCCTTATTCTAGCATCTGCCG CATCCCACCTGCCAGGGCCAGCCCTCTCCCCACCTGCTCGGAGGTACCCGGACCCTTCCCAGTGGATAGAAAGAGATACAGTGGAAACTCTGGGTAA GAC 28 FJX_270353_1:1 q TTTAAAAAATTAAGTTTCTTAGCTATTGCTGTTGGGAGTGTAAATTGTTGAAGCTTTTCTGGAATATGTTGTAGAAATAGGTACTCAGAGCCTTTTAAGT G TCCATTTATGCTGATCTGATAAATAATTTCTAAAAATGTAAGGAGAAAAACCAAGAAATAAAGATTTAAATATAAGCATGTTCGAAGCATATTTTAA GGA 29 FJX_278702_1:1 q ACAGGACAAAATCAAACCCCAAAAAACCCCAAAAAATCAATTCTTCAAGAAAATTTTAAGGGGAAAAAAAAGAAGAGGTAACCTACAGATCAAAATAGAC TT AGAGATATGTCAACTGACTATATCAGTCAGTCAGTCAGTTCAGTGGCTCAGTCATGTCCAACTCTTTGCAACCCCATGAATCGCAGCACGCCGGGCCT CC FJX 278848 q GCTCAGTCATGTCCAACTCTTTGCAACCCCATGAATCGCAGCACGCCGGGCCTCCCTGTCCATCACCAACTCCCAGAGTTCACCCCGACTCACGTCCATC G AGTCAGTGATACCATCCAGCCATCTCATCCTCTGTCATCCCCTTCTCCTCCTGCCCCCAATCCCTCCCAGCCTCAGAGTCTTTTCCAATGAGTTAACTCT 31 FJX_PLAPROTRI_l: i q TTTCGGGCTGGACGCGGCGGGCCTGGAAATTGTTTCCTTACGCCTTTTTCCATCAGCCATTGTAGTGTATGCGCTGCCCCTGCAGCCCAACTTGCAGCTG(CCG)9 ACGTCGCGCCCGCCGCCGCCTCCCCCGGCCCGCTCGGCCACCGCCCCCGCCTCCCGCCGGCCCGCCCCCTGCCCCCGCCCGCCGCTGCGCCCGCCCTC CC 32 FJX_PLAPROSNP_l: i q TTTTTCCATCAGCCATTGTAGTGTATGCGCTGCCCCTGCAGCCCAACTTGCAGCTGCCGCCGCCGCCGCCGCCGCCGCCGCCGCCGCCGACGTCGCGCCC A CCGCCGCCTCCCCCGGCCCGCTCGGCCACCGCCCCCGCCTCCCGCCGGCCCGCCCCCTGCCCCCGCCCGCCGCTGCGCCCGCCCTCCCCGCCCTTCCCC G 33 FJX_307148_1:1 GTAAAAGACACATTATTTGATCAAATTTCATGTGTATAGAAAGTGAAAACAGATAGTTTGTGTCCCTTCTCTTCATCTGGTCCTGGTGCAGGAGTATTAC G ACGGCTGCTGACACAAACTGTCCAGCTCGTTGGTCCTAGTCCAAAACCCAGGTCAGAGGTTTTTAAACTATGCAACCTTTATTTTCAAATGTGATTCT AT 34 FJX_308585_1:1 q TTGATGTAGTCCATTTTGTGGTGGTGGTTGTAGCTTTTGGTGTCATGTCTACAAAAAAGTCCAAAGGTTAACATCTGTGTTTCTTCTGCTGCTGCTAAGT C GCTTCAGTCGTGTCCGACTCTGCGACCCCATAGATGGCAGCCCACCAGGCTTCGCCATCCCTGGGATTCTCCAGGCAAGAACACTGGAGTGGGTTGCC AT Q - nucleic acid sequence Q allele q - nucleic acid sequence q allele SHR507647NZPR 301525122 Materials and methods Genotyping. Genotyping the 56 microsatellite markers on the HFxJ F2 pedigree was performed using standard procedures labeling the PCR products with fluorescent primers and size fractionating them on a ABI36730 capillary sequencer (e.g. Coppieters et al., 5 1998a). Genotyping the 925 BTA14 SNP markers on the HFxJ F2 pedigree was performed with a custom-made Golden Gate assay on an Illumina Beadstation500 using standard procedures. Genotyping the pQTN was done using Taqman assays-by-design using standard procedures (ABI, Foster City, CA). The primers used to analyse the genetic markers are listed in Table 7 below.
QTL mapping and fine-mapping. Weight/Height: (i) QTL mapping in the HFxJ F2 population under a line-cross model was conducted with a regression based statistical approach (Haley et al 1995). (ii) QTL mapping in the HFxJ population under a half-sib pedigree model was conducted with HSQM (Coppieters et al., 1998b). Significance 15 thresholds were determined by permutation testing following Churchill & Doerge (1994). Confidence intervals for the QTL were determined by bootstrapping following Visscher et al. (1996). (iii) QTL fine-mapping was achieved by simultaneously exploiting linkage and LD using a haplotype-based approach described recently by Druet & Georges (2009). The mixed model includes a mean (fixed), a haplotype effect (random), an animal effect 20 (random) and an error term (random). The haplotype effect corresponds to haplotype clusters defined using a Hidden Markov Model (HMM). The covariance between distinct clusters is assumed to be zero. Variance components and individual effects are estimated using AI-REML. Results are reported as LRT=21nLR, where LR corresponds to the ratio between the likelihood of the data under HI (model with haplotype effect) and likelihood 25 of the data under HO (model without haplotype effect). LRT is asymptotically distributed as a chi-squared variable with ~ one degree of freedom. In some instances, two QTL positions were fitted simultaneously. Evidence in favor of a QTL at position 2 (conditional on the presence of a QTL at position 1) was then evaluated from the increase in the LRT achieved by adding this second QTL in the model, (iv) Association analysis was conducted 30 using a mixed model including a mean, a regression on the number of alleles "1" for the considered SNPs (fixed), an animal effect (random) and an error (random). Variance components and effects were estimated by REML. The significance of the SNP effect was estimated using an F-test. In some instances, association analysis only considered the SHR507647NZPR 301525122 marker allele inherited from the mother. Transcriptomic data: (i) eQTL mapping in the HFxJ population under a half-sib pedigree model was conducted with HSQM (Coppieters et al., 1998b) as for the phenotypic data, (ii) Association analysis was conducted using a linear model including a mean, an experiment (1 or 2) effect (fixed), a regression on the 5 number of alleles "1" for the considered SNPs (fixed), and an error term. The significance of the SNP effect was estimated using an F-test. In some instances two distinct SNP effects were fitted simultaneously. The significance of the second SNP effect conditional on the first one was estimated by an F-test. Haplotype-based association analysis was conducted by multiple regression against the number of copies of each of the 20 HMM-defined 10 haplotype clusters (Druet & Georges, 2009) Sequencing. Resequencing the ~750 Kb critical region was achieved by amplifying the entire interval as 103 long range PCR products, equimolar pooling of amplification products, size-fractionation by nebulization, ligation of adapters including multiplex 15 identifiers, and sequencing on a Roche FLX system using standard procedures. The flow grams were analyzed using the GS Sequencer and Reference Mapper software modules from Roche.
Table 7A Assay ID FJX 21374 1:1 FJX 221494 FJX 221762 FJX 222352 1:1 FJX_240394 FJX 247103 FJX_250879_1:1 FJX_253666 FJX_263178_1:1 FJX 267438 2:2 FJX 270353 FJX 278702 1:1 FJX 278848 FJX_PLAPRO TRI 1:1 FJX_PLAPROS NP 1:1 Genotyping Forward Primer method Name (SEQ ID No) Taqman FJX_21374-161F (35) Taqman FJX_221494-FJX2F (37) Taqman FJX_221762-2195F (39) Taqman FJX_222352-FJX3F (41) RFLP-PCR FJX_240394-1 F (Haelll) (87) Re- FJX_247103-1F sequencing (89) Taqman FJX_250879-2531F (43) Taqman FJX_253666-FJX9F (45) Taqman FJX_263178-FJX4F (47) Taqman FJX_267438-RPTF (49) Taqman FJX_270353-FJX5F (51) Taqman FJX_278702-FJX6F (53) RFLP-PCR (Taq I) Forward Primer Seq.
CCTTCAAGAAGGTGGAGAAGAAAACA AAACCAGCAATCGCAAAACACTT AGTCATGTATGGGTGTAAGAGTTGGA GCTTTCCTCCATCCCAATGTTAAATTTTAC TA tattcattggaaggactgatgcta CATTAGGCACTGTGCTGTTTATCT GATCTCAGTAGGCTGGAAGCA AAAACACTGAATGAGAAAGATACCACAGT GCTAATTAAGGGCAAAAGCTATCATAATC AA CTGTTCCACCAGGATCCTTATTCT TGTTGTAGAAATAGGTACTCAGAGCCT CAAAAAACCCCAAAAAATCAATTCTTCAAG AAAA fjx_278848 (91) atccaacagagatgcagttaacaa FJX 307148 1:1 FJX 308585 1:1 Microsat Fragment length analysis 5UTR_SNP3F (55) Re sequencing (Sanger method ABI3730) Taqman Taqman 5UTR_SNP3F (85) FJX_307148-3031F (57) FJX_308585-FJX7F (59) agaactcaccgcg666ctttaacat AGAACTCACCGCGGGGCTTTAACAT CAGATAGTTTGTGTCCCTTCTCTTCAT AGTCCAAAGGTTAACATCTGTGTTTCT Reverse Primer Name (SEQ ID No) FJX_21374-161R (36) FJX_221494-FJX2R (38) FJX_221762-2195R (40) FJX_222352-FJX3R (42) FJX_240394-1 R (88) FJX_247103-1R (90) FJX_250879-2531 R (44) FJX_253666-FJX9R (46) FJX_263178-FJX4R (48) FJX_267438-RPTR (50) FJX_270353-FJX5R (52) FJX_278702-FJX6R (54) FJX_278848 (92) 5utr_snp2r (56) Reverse Primer Seq.
GACTGAACTGAACTGAACTGAACTG CAGAAGGGAGAAGCAGAAACCT CAACACCACAGTTCAAAAGCATCAA GTCCCTTTTCCAGTTATCTGCATGA AACAAGGAAACATGCTGAAATACA cgctatataccacatggagctaga CACACTGTCTTCCCAGTATGTAAGT CACCAAAAGCAGAAAACAGAACAAC ATGTTG ACATTTTCTTAATTGTTTACTAAT GAGT TTTATGA GGAAGGGTCCGGGTACCT TTAAATCTTTATTTCTTGG11 I IICTCCTTACAT TTTTAGAA CACTGAACTGACTGACTGACTGATA Taattgatgcttttgaactgtggt ggaggagcgcgcggggaagg 5UTR_SNP2R (86) FJX_307148-3031R (58) FJX_308585-FJX7R (60) GGAGGAGCGCGCGGGGAAGG ACGAGCTGGACAGTTTGTGT TGGTGGGCTGCCATCTATG Table 7B Assay ID Reporter 1 Name (SEQ ID No) Reporter 1 Dye Reporter 1 Sequence FJX 21374- FJX_21374_1:1 161V2 (61) VIC TCTGCCAGACTCTCC FJX 221494- FJX_221494 FJX2V1 (63) VIC TGTGCATCTCACCCCCT FJX 221762- FJX_221762 2195V1 (65) VIC AAGAACACTAAGCACCAAAG FJX 222352- FJX_222352_1:1 FJX3V1 (67) VIC CTTCTCACATCATGTAATC FJX_240394 FJX_247103 FJX 250879- FJX_250879_1:1 2531V1 (69) VIC ACTTGGGTCAATATTT FJX 253666- FJX_253666 FJX9V2 (71) VIC CCGGCTCACCAACCA FJX 263178- FJX_263178_1:1 FJX4V1 (73) VIC TTTCAGTCTCCTGTATGCTA FJX 267438- FJX_267438_2:2 RPTV1 (75) VIC CATCTGCCACATCCCA FJX 270353- FJX_270353_1:1 FJX5V2 (77) VIC CAGCATAAATGGAGACTTAA FJX 278702- FJX_278702_1:1 FJX6V2 (79) VIC TGACATATCTCTCAAGTCTATT FJX_278848 FJX_PLAPROTRI_1:1 - - - FJX_PLAPROSNP_1:1 _ - .
FJX 307148- FJX_307148_1:1 3031V2 (81) VIC CAGCAGCCGTTGTAAT FJX 308585- FJX_308585_1:1 FJX7V2 (83) VIC ACTGAAGCAACTTAGC Reporter 2 Name (SEQ ID No) FJX_21374-161M2 (62) FJX 221494-FJX2M1 (64) FJX_221762-2195M1 (66) F JX_222352-F JX3M1 (68) FJX_250879-2531M1 (70) FJX_253666-FJX9M2 (72) FJX 263178-FJX4M1 (74) FJX_267438-RPTM1 (76) FJX 270353-FJX5M2 (78)" FJX_278702-FJX6M2 (80) FJX_307148-3031M2 (82) FJX_308585-FJX7M2 (84) Reporter 2 Dye Reporter 2 Sequence Design Strand FAM TGCCAGGCTCTCC Reverse FAM CTGTGCATCTTACCCCCT Forward FAM AACACTAAGCGCCAAAG Forward FAM CTCACATCGTGTAATC Forward FAM C11GGGTGAATA111 Forward FAM CGGCTCGCCAACCA Reverse FAM TTTCAGTCTCCTTTATGCTA Forward FAM TCTGCCGCATCCCA Forward FAM AGCATAAATGGACACTTAA Reverse ^ FAM TGACATATCTCTCGTCTATT Reverse FAM CAGCAGCCGTCGTAAT Reverse FAM CTGAAGCGACTTAGC Reverse SHR507647NZPR 301525122 60 Example 2 SNP Effect and Breeding Value Calculations The effect of the C allele for the FJX_250879_1:1 SNP has been demonstrated to increase 5 live weight in mature dairy cattle by approximately 10-20 kg in comparison to the G allele. The inventors have demonstrated this effect in 3 breeds of cattle (Holstein-Friesian, Jersey and Ayrshire). Live weight has a negative relative economic value in the New Zealand dairy breeding index; Breeding Worth (BW). This reflects that a larger cow requires more energy for maintenance. Animals that are smaller but produce the same 10 output of milk solids are desirable.
Analysis of the genotype of the FJX_250879_1:1 SNP (and other genetic markers identified herein) may assist in calculations of the BW of an animal. The effect of the genotype of the FJX_250879_1:1 SNP for 34 dairy traits that are of economic or farmer 15 importance was studied.
Two thousand one hundred and fifty Holstein-Friesian (HF) progeny tested sires were genotyped for FJX_250879_1:1 and FJX_307148_1:1. There was total concordance between the genotypes for the two SNPs, i.e. all animals that had a C allele for 20 FJX_250879_1:1 had a corresponding A allele for FJX_307148_1:1. Therefore, the results that are presented here are applicable to both SNPs.
Genotype results of the 2150 HF sires identified that the A allele frequency for FJX_307148_1:1 is 0.87. There were only 38 sires that have the GG genotype for the 25 FJX_307148_1:1 marker. For the statistical analysis the allelic effects were estimated between animals with the AA genotype (1678 animals) and AC genotype (479 animals).
The model fitted included the genotype as a fixed effect and also fitted year of birth and percent of North American genetics as covariates.
The effect of the genotype of the FJX_307148_1:1 SNP for 34 dairy traits is shown in Table 3 below. There are a number of traits that are significant at the 5% threshold level. For example, live weight has a P-value of <0.0001 with the effect of the G allele being 61 .9 kg less than the A allele. Other traits that are significant at the 5% level are; body condition score, calving difficulty, capacity, dairy conformation, fat percent, fat persistency, fore udder, gestation length, legs, milk volume, overall opinion, protein yield, protein persistency, rear teat, rump angle and width, somatic cell, stature and temperament.
Table 3: Allelic effect of the G allele compared to the A allele for marker FJX_307148_1:1 for 34 dairy traits.
Trait Trait Effect p-value 1 Adaptability -0.63 0.37 2 Body Condition Score -0.03 <0.0001 3 Breeding Worth 2.11 0.5 4 Calving Difficulty -1.60 <0.0001 Capacity -0.05 0.004 6 Dairy Conformation -0.03 0.03 7 Fat % 0.04 0.009 8 Fat yield -0.13 0.83 9 Fat Persistency 0.36 0.02 Fertility -0.003 0.99 11 Fore Udder -0.04 0.02 12 Front Teat 0.02 0.27 13 Gestation length -0.51 0.02 14 Legs 0.01 0.02 Live weight -10.9 <0.0001 16 Milk volume -45 0.003 17 Milk Persistency 2.9 0.47 18 Milking Speed 0.02 0.18 19 Overall Opinion -0.04 0.004 PM21 -0.82 0.27 21 Protein % 0.0002 0.98 22 Protein yield -1.57 0.0003 23 Protein Persistency 0.31 0.004 24 Rear Teat 0.05 0.02 Rear Udder 0.03 0.17 26 Residual Survival -0.88 0.13 27 Rump Angle -0.06 0.0001 28 Rump Width -0.15 <0.0001 29 Somatic Cell 0.06 0.001 -0.20 <0.0001 Stature SHR507647NZPR 301525122 62 31 Temperament -0.05 0.0007 32 Total Longevity -3.56 0.72 33 Udder Overall -0.004 0.83 34 Udder Support 0.001 0.95 Materials and Methods Genotyping As mentioned above, two thousand and fifty Holstein-Friesian progeny tested sires were 5 genotyped for the SNPs FJX_250879_1:1 and FJX_307148_1:1. Genotyping was conducted using Taqman assays-by-design using standard procedures (ABI, Foster City, CA). The assays were performed using the relevant primers herein before described in Example 1 (Table 7).
The statistical model fitted for each of the phenotypes (breeding values) was a linear model with genotype as a fixed effect and year of birth and percent of North American genetics fitted as covariates.
Breeding Values Estimated breeding values were calculated for a number of traits in progeny tested sires, in accordance with the following information.
Calculation of Breeding Values Estimated Breeding Values (EBVs) are an estimate of a bull or cow's genetic merit for any 20 given trait. Such genetic merit may also reflect somatic cell and/or germ cell characteristics, depending on the trait(s) being examined. EBVs can therefore not only predict phenotype performance and herd lifetime productivity (that is, provide a genetic diagnostic of phenotype) but also provide predictive insight into the ability of an animal to pass on superior genetic material (such as increased protein yield) to its offspring.
Breeding values are calculated/estimated for production traits (expressed in units of measurements, that is, litres of milk and kilograms of fat and protein), fertility traits, conformation traits, as well as other health, management and survival traits; for example, individual type traits, liveweight, somatic cell count, daughter fertility, calving difficulty 30 and gestation length.
SHR507647NZPR 301525122 63 In reference to BLUP: 'Best' is a reference to the method which gives the best estimates for the breeding value, or to put it more precisely, this minimizes the variance of difference between the estimates and the true breeding values. 'Linear' means that there is a linear 5 relation between parameters in the statistical model. 'Unbiased' means that the breeding value estimates are central, that they are expected to be normally distributed with the true breeding value as mean value. 'Prediction' normally refers to the future. But prediction also refers to the estimation of realized values of a random variable drawn from a population with known variance- and co-variance-structure. A simpler and less precise formulation: 'prediction' is estimation of a given value of a random variable. It is better to use the expression 'prediction of breeding values', but 'breeding value estimates' are commonly used. With BLUP the breeding values are calculated for all animals simultaneously, and at the same time environmental effects can be corrected for.
Animal Evaluation Animals were evaluated according to New Zealand's current national genetic evaluation system to dairy cattle, known as "Animal Evaluation" (www.aeu.org.nz). The genetic evaluation system is conducted with a common base for all breeds and crosses. The genetic base for all evaluations apart from Calving Difficulty is a group of 1985-born cows which had all traits recorded in 1987. These cows were of all breeds and crosses. The New Zealand national breeding objective is to identify animals whose progeny will be the most efficient converters of feed into farmer profit. The evaluation system is designed to identify the most profitable and efficient sires regardless of breed.
Trait Evaluation Production Traits - The current models for the national genetic evaluation of milk volume, fat yield and protein yield are multiple lactation (ML) random regression (RR) test-day (TD) models where lactations 1,2 and 3 are modelled as individual lactations and lactations 4, 5 and 6 together are modelled as one lactation under the same genetic control.
Each lactation (1-6) is considered a different trait for modelling the permanent environmental effects. These statistical models include effects for: - herd-year-season-test-day contemporary group, with season referring to spring or autumn calving period; SHR507647NZPR 301525122 64 - age at calving in days (nested within breed); - induced lactation; - heterosis class; - stage of lactation fixed effects of days in milk at test-day (standardised 5 Legendre polynomials of order 2 to 5); - random additive genetic effect of the animal at each lactation day (four lactations and standardised Legendre polynomials of order 1 to 3); - random non-additive genetic and permanent environment effect of the animal at each lactation day (six lactations and standardised Legendre polynomials of order lto 3).
The BVs reported from these models are calculated as the simple average of the four 270-day lactation BVs. For the national evaluation of the production traits, genetic groups were assigned by breed, sex of missing parent, birth year and country of origin. Four breed 15 classes were assigned genetic grouping, namely, Holstein-Friesian, Jersey, Ayrshire-Red, and other breeds. Genetic groups were assigned in five year intervals from 1960 to 1980 then yearly, with the first birth year group being prior to 1960. Country of origin was defined as New Zealand, North American and Other. Missing male parent of New Zealand origin had two categories: (i) unknown male but known from mating records to be an 20 artificial insemination proven sire, or (ii) completely unknown male. If genetic groups had less than 200 animals per group, birth years were clustered. No clustering occurred across breed, origin or sex genetic groups. Foreign information for sires for is blended with New Zealand information for calculating the breeding values for the dairy production traits (and for the female fertility, udder health, longevity and conformation traits). The blending 25 procedure uses data from the International Genetic Evaluation Service provided by Interbull (http://www.interbull .slu.se).
Liveweisht - The statistical model for liveweight is a repeated record, single trait, additive effects repeatability model. It includes effects for herd-year-season-age contemporary 30 group, age at calving in months (nested within breed), stage of lactation when weighed (nested within age), heterosis, genetic merit of the animal, and random non-additive genetic and permanent environment effects.
SHR507647NZPR 301525122 65 Cow Fertility - The objective for herd reproductive performance in most New Zealand herds is to achieve high pregnancy rates in a short time period following the planned start of mating, and to maintain calving intervals very close to 365 days. In this system successful reproduction depends on two factors which display genetic variation. The first 5 factor is the ability of the cow to resume cycling soon after calving, and to be mated early in the herd's mating period. The second factor is the cow's ability to conceive, sustain a pregnancy and calve early in the herd's subsequent calving period. Animal Evaluation has developed a genetic evaluation of cow fertility that incorporates both these aspects of successful reproductive performance in seasonal dairying. Mating records - Being 10 presented for mating in the first 21 days of the herd's mating period (PM21) is scored 1 for the cows that are mated in this period, and 0 for cows which fail to be mated. Calvins records - In the New Zealand circumstances, successfully bearing a calf in the first 42 days of the herd's calving period is the best indicator of successful reproduction (CR42, standing for Calving Rate in the first 42 days of the herd's calving period). CR42 is scored 15 1 for cows that successfully re-calve in the first 42 days, and 0 for cows that fail to re-calve in this period. CR42 is coded "missing" for cows which leave the herd prior to the re-calving period for reasons other than low fertility, or which have not yet had the opportunity to re-calve to the recorded mating. The analysis includes eight traits: PM21 in the first three parities, and the associated CR42 in the second, third and fourth lactations — 20 together with Body Condition Score at day 60 of the first lactation, and first lactation milk yield relative to contemporaries. These records are evaluated using a multiple trait animal model (www.aeu.org.nz/page.cfm?id=59&nid=35). The reported BV relates to the comparative percentage likelihood of a bull's daughters to re-calve for their second lactation in the first 42 days of the herd's calving period.
Somatic Cells in milk- Test day records for somatic cell counts are transformed into somatic cell scores (SCS) by taking the log (base 2) of test day SCC/1000. SCS is analysed in a multiple trait random regression animal model. The two traits are first lactation SCS — and second and third lactation SCS analysed as repeated observations of a single trait. 30 The statistical model for analysis of a cow with SCS records includes effects for herdyear-season-age-testday contemporary groups, induced lactation, heterosis, age at calving (in months nested within breed class), stage of lactation, genetic group, genetic merit of the SHR507647NZPR 301525122 66 animal, random non-additive genetic and permanent environment effects, and random residual effects. Results for breeding values averaged over the two traits are reported.
Herd-life - Herd-life is defined as the interval from the date a cow has her first calf to the 5 date when she has her last herd test, and is recorded in days. It is evaluated using a single-trait animal model. The multiple-trait (MT) animal model used for the national genetic evaluation of survival contains 4 survival traits (SV12, SV13, SV14, SV15) and 9 predictor traits. The model can cope with missing data on any combination of traits. SY15 is the trait reported by the AE system (converted into days of herd life). Each individual 10 trait record was modelled as: yijkl = hysij + htiswnshs=16'Z+ aik + eijkl where yijkl is the record for ith trait, hySij is the jth herd-year-season fixed effect for a cow's first lactation for trait i, with season referring to spring or autumn calving period, htis is the linear regression coefficient for the sth heterosis effect for trait i, Whns is the sth heterosis covariate for animal n, aik is the random additive genetic effect of animal k for trait i, eyki is the random residual associated with record yijkl.
Genetic correlations between two survival traits and 9 predictor traits (all significantly different from zero) are summarised in Table 4.
Table 4: SV12 SV15 Protein yield 0.424 0.418 see -0.156 -0.040 CR42 0.841 0.478 BCS 0.370 0.240 Milking speed 0.072 0.220 Overall opinion 0.330 0.379 SHR507647NZPR 301525122 67 Legs -0.119 -0.074 Udder overall 0.113 0.113 Dairy conformation 0.235 0.183 Herd-life breeding values are reported as Total Longevity.
Residual Survival - For inclusion in the total economic merit index called Breeding Worth 5 (BW), Residual Survival has been defined in a way that ensures that herd-life is counted only once in the index. Residual Survival is defined as "Herd-life after accounting for the genetic effects of production, liveweight, fertility and milk somatic cells on herd life." The Residual Survival trait included in BW is calculated from the following equation, where EBV stands for Estimated Breeding Value.
EBV(Total Long'ty) = 5.434*EBV(Milkfat) + 4.408*EBV(Protein) + 0.03815*EBV(Milk) - 0.489*EBV(Liveweight) + 27.847*EBV(Fertility) - 65.13 l*EBV(Somatic Cells) + Residual Survival The estimated Residual Survival Breeding Values calculated by this method are uncorrelated with the production, liveweight and fertility traits included in the BW index. This property is important in order to ensure that effects on herd life are counted only once in the BW index.
Management and Conformation Traits - The models for the linear management and conformation traits (apart from Udder overall) are single record, multiple trait, additive genetic effects models. The statistical model for analysis of a cow with linear type scores includes effects for herd-year-season contemporary group, stage of lactation class when scored and age at first calving class (in months nested within breed), heterosis, genetic group, animal genetic merit and the random residual. The four farmer-scored management traits (Adaptability to milking, Shed temperament, Milking speed, Overall opinion) are evaluated together in a multiple trait evaluation that takes the genetic correlations amongst the four traits into account. The inspector-scored traits associated with body conformation (Stature, Capacity, Rump angle, Rump width, Legs, Dairy conformation) are evaluated together in a multiple trait evaluation that takes the genetic correlations amongst these six SHR507647NZPR 301525122 68 traits into account. The inspector-scored traits associated with udder conformation (but excluding Udder overall) are analysed together in a multiple trait evaluation. These traits are Udder support, Front udder, Rear udder, Front teat placement, and Rear teat placement. They are evaluated together in a multiple trait evaluation that takes the genetic correlations 5 amongst the five traits into account. Udder overall is evaluated in a single trait model. Breeding Values for individual traits are expressed in the units in which the trait is measured.
Calvinz Difficulty - Calving Difficulty Breeding Values are supplied to help Artificial 10 Breeding organisations assess the suitability of bulls for mating with yearling heifers; and to give farmers knowledge about bulls which cause higher than usual rates of calving assistance when mated to their cows and heifers. A sire's Calving Difficulty Breeding Value (BV) predicts the percentage of assisted calvings expected when he is mated to yearling heifers. The higher the BV, the higher the expected percentage of assisted 15 calvings. The average BV of sires born in 1985 is set to zero. All breeds have been evaluated together for the calculation of Calving Difficulty BVs. Although the Calving Difficulty BV is expressed in terms of assisted births in first-calving heifers, the BV can also be used to identify bulls that are expected to increase rates of calving assistance for cows carrying the bulls' calves. The records used for the analysis comprise over 5 million 20 records of calving assistance collected from New Zealand herds since 1994. The records were extracted in January 2006. The BV for Calving Difficulty was estimated with a multiple-trait sire model using BLUP methodology (best linear unbiased prediction). The two traits were: - assistance for calves born to heifers, and - assistance for calves born to cows. The statistical model used in the estimation was: Yjjkim = HERD-YEAR-AGEi + CALF-SEXj + CALF-BREEDmk+ HETEROSISm + SIREi + Residualjjkim in which Yijkim is the observation of calving assistance (unassisted = 0, assisted = 1) of CALFm, in 30 HERD-YEAR-AGE GROUPi, of SEXj, with CALFm's BREED GROUPSk, BREED HETEROZYGOSITY of CALFm, and with SIREi Residual,jnm is the variation that is not explained by the model.
SHR507647NZPR 301525122 69 Body Condition Score Breeding Values - Body condition score (BCS) is commonly used as a method to assess body energy reserves. In recent years there has been considerable interest in assessing the potential use of BCS in dairy cattle breeding programmes. In particular, breeders are interested in knowing about sires that transmit lower or higher BCS 5 for their daughters in early lactation when cows' body reserves are being mobilised for lactation at the same time as the cows are expected to get back in calf. Consequently Estimated Breeding Values for day 60 of first lactation heifers have been calculated for all AE Enrolled sires.
* The Dexcel BCS scale of 1 to 10 is used in NZ 10 * A cow with a BCS of less than 3.0 is considered emaciated * A cow with a BCS greater than 7.0 is considered obese (Kevin Macdonald and John Roche, Condition Scoring Made Easy, Dexcel, 2004.) Persistency Breeding Values - Persistency Breeding Values have been calculated for 15 milkfat, protein, and volume. The Persistency BVs are estimates of the genetic ability of cows to maintain production subsequent to peak of lactation. They are reported on a scale where higher values correspond to greater persistency of lactation after the peak. Describing the mathematical technique used is beyond the scope of this Introduction. Two explanatory points are relevant. Cows with persistent lactation curves can have high total 20 lactation yields, average total lactation yields or low total lactation yields. The persistency measure is designed to be independent of total yield - so farmers ought not to use high Persistency Breeding Values in an attempt to select for higher yields. Persistency Breeding Values are reported in the same units as the dairy production traits - kilograms of milkfat, kilograms of protein, and litres of volume. To a close approximation, a Persistency BV of 25 +5 kg can be interpreted as the genetic predisposition to yield 5 kg more in the lactation period after the 115th day of lactation compared to the earlier lactation period. A Persistency BV of -5 kg can be interpreted as the genetic pre-disposition to yield 5 kg less in the lactation period after the 115th day of lactation compared to the earlier lactation period. The total lactation period used for the evaluation is 270 days.
Gestation length - Gestation length is the number of days from date of insemination to the date of parturition. The average number of days for Gestation length in dairy bovine is 282 SHR507647NZPR 301525122 70 days. A breeding value of -10 can be interpreted as an animal that will leave progeny on average 5 days shorter than the dairy bovine average.
Economic Evaluation To compare individual animals on net farm profitability for breeding replacements, the Breeding Worth (BW) has been developed. The BW is the sum of the Breeding Values for milkfat, protein, milk volume, liveweight, fertility, milk somatic cells and residual survival each weighted by an economic value. The BW economic values for each trait represent the expected net income per unit of genetic change (per unit of feed) from breeding 10 replacements. The economic values are calculated using a bio-economic devised to value technological change. The model includes income streams from milk production, cull cows and bobby calf sales and cost streams associated with maintaining and growing cows and replacements, the feed required for production and dairy cash expenses. Predictions of future milk component prices are taken into account. The prices and costs used in the farm 15 model are reviewed annually. The unit of feed adopted for reporting economic values is 4.5 tonnes of dry matter of feed containing 10.5 megajoules of metabolisable energy per kilogram. The economic values used in the BW are given in the following table.
Milkfat BV (kg) Protein BV (kg) Milk BV (1) Liveweight BV (kg) Fertility BV (%) Somatic Cell BV (score) Residual Survival BV (days) 1.646 7.846 -0.083 -1.171 2.790 -28.924 0.042 cBreedinz animals for profit and efficiency Selecting animals on Breeding Worth identifies the animals expected to be the most profitable and efficient for breeding replacements.
Breeding Worth BW ($) - The expected ability of an animal to breed replacements which are efficient converters of feed into profit. A Breeding Worth of 206 indicates the bull is expected to generate an extra $206 profit per year per unit of feed, through breeding replacements, compared with using a bull with a BW of zero.
Reliability - Associated with the BW, the reliability figure is the amount of confidence we can place in the figure. The more information included in the evaluation, the greater the SHR507647NZPR 301525122 71 reliability and less likely it is to change with additional records. Reliability ranges from 0%, meaning we know nothing about the animal or any of its ancestors, to 99%. If a sire has a Milkfat Breeding Value of +25 with 80% reliability it is expected that his future Milkfat Breeding Value would be in the interval of+15 to +35 in the great majority of 5 cases. If a sire has a Milkfat Breeding Value of +25 with 95% reliability it is expected that his future Milkfat Breeding Value would be in the interval of +20 to +30 in the great majority of cases. The reliability reflects the amount of information used in the calculation of the Breeding Value.
Breeding Value - The genetic merit of an animal for individual traits relative to a base zero (1985 born cow). A Breeding Value of+10 kg protein indicates the bull will produce daughters, which on average are genetically superior by 5kg protein per lactation above the base cow (a bull can pass 1/2 of its genetic merit on average onto offspring).
Fat %, Protein % - The Fat % and Protein % are Breeding Value estimates for the sire. The Breeding Values for Fat % and Protein % are not expressed as deviations but include the base averages (1985 cows) The base averages are: Fat % = 4.79 Protein % = 3.69 The Protein % and Fat % Breeding Values are across breed values.
Traits Other than Production The Traits Other than Production (TOP) Evaluation System is a national scheme for assessment of non-productive characteristics of dairy cattle (bulls and cows) based on linear assessment. Its main objective is to provide bull and herd owners with accurate and easy-to-use information for decision making.
Daughters are evaluated for 16 traits using a linear assessment on a scale from 1 to 9 where 1 and 9 represent the biological extremes. The traits are scored across breed and are defined as follows: SHR507647NZPR 301525122 72 Traits scored bv the herd owner: Adaptability to milking - describes how soon the animal settled into the milking routine after calving.
Slowly 1 5 9 Quickly Shed Temperament - describes the temperament of the animal in the shed while being handled and milked. It is a different trait to adaptability to milking and should be assessed once animals have settled into the milking routine.
Vicious 1 5 9 Placid Milking Speed - describes the milking speed of the animal, i.e. the time from putting cups on to the time milk flow stops or the cups are taken off.
Slow 1 5 9 Fast Overall Opinion - describes the farmer's overall acceptance of the animal as a herd member.
Undesirable 1 5 9 Desirable Traits scored bv the Inspector: Stature - describes the height at the shoulders of the animal. Each score represents 5 cm height at the withers.
Under 105cm 1 5 9 Over 140cm Capacity - describes the capacity of the animal as a combination of strength and depth of 25 chest and body as viewed from side, rear and front in relation to the physical size of the animal.
Frail 1 5 9 Capacious Rump Angle - describes the angle of a line between the centre of the hips and the top of the pins.
Pins high 1 5 9 Pins low/Sloping SHR507647NZPR 301525122 73 Rump Width - describes the width of the pins, hips and thurls relative to the size of the animal.
Narrow 1 5 9 Wide Legs - describes the straightness or curvature of the back legs from an imaginary line between the thurls and the mid hoof while the animal is walking.
Straight 1 5 9 Sickled/Curved Udder Support - describes the strength of the suspensory ligament as viewed from the rear. 10 It also includes udder depth relative to the hocks.
Weak 1 5 9 Strong Front Udder - describes how well the front udder is attached to the body wall.
Loose 1 5 9 Strong Rear Udder - describes the height and width of the rear udder attachment, as distinct from udder support.
Low 1 5 9 High Front Teat Placement - describes the placement of the front teats (at the point of attachment to the udder) relative to the centre of the quarter as viewed from the rear.
Wide 1 —- 5 -— 9 Close Rear Teat Placement - describes the placement of the rear teats (at the point of attachment 25 to the udder) relative to the centre of the quarter as viewed from the rear.
Wide 1 5 9 Close Udder Overall - All traits pertaining to the udder including those traits that have been linearly scored.
Undesirable 1 5 9 Desirable Dairy Conformation - All traits pertaining to dairy conformation including those body traits that have been linearly scored, but excluding all the udder traits.
SHR507647NZPR 301525122 74 Undesirable 1 5 9 Desirable All two-year-old daughters in the herd are inspected to avoid bias through selection of daughters and to allow valid comparisons.
Example 3 Further analysis of the FJX_250879_1:1 SNP in the Ayrshire and Jersey populations was undertaken for live weight using breeding values. One hundred and eighteen Ayrshire sires were genotyped with an allele frequency of 71% for the G allele. In this population, 10 significant differences (P <0.05) were found between the 3 genotypic groups for live weight. The average allelic effect for the C allele in the Ayrshire population is 17.3 kg compared to the G allele.
Table 5: Genotyic effects on live weight in the Ayrshire population for marker 15 FJX_307148_1:1 relative to the CC genotype.
Genotypic state Live weight (kg) Number genotyped CC 0.0 12 CG -23.3 44 GG -34.7 62 In the Jersey population 1309 sires were genotyped for the FJX_250879_1:1 SNP with an allele frequency of 98.5% for the G allele. Only 2 sires were CC genotype, 35 CG and 1272 GG. The live weight effect between the GG and GC genotypic groups was 20 statistically significant for live weight (P <0.05). The allelic effect for the C allele in the Jersey population is 17.9 kg compared to the G allele.
Table 6: Genotyic effects on live weight in the Jersey population for marker FJX_307148_1:1 relative to the GC genotype.
Genotypic state Live weight (kg) Number genotyped CG 0.0 GG -17.9 1272 SHR507647NZPR 301525122 75 The results from these two breeds demonstrate the SNP effects are consistent across breeds.
Materials and methods Genotyping was conducted using Taqman assays-by-design using standard procedures (ABI, Foster City, CA). The assays were performed using the relevant primers referred to herein before in Example 1 (Table 7).
The statistical model fitted for each of the phenotypes (breeding values) was a linear 10 model with genotype as a fixed effect and year of birth fitted as a covariate.
The invention has been described herein, with reference to certain preferred embodiments, in order to enable the reader to practice the invention without undue experimentation. However, a person having ordinary skill in the art will readily recognise that many of the 15 components and parameters may be varied or modified to a certain extent or substituted for known equivalents without departing from the scope of the invention. Accordingly, as noted herein before, it should be appreciated that such modifications and equivalents are herein incorporated as if individually set forth. Similarly, the invention also includes all of the steps, features, compositions and compounds referred to or indicated in this 20 specification, individually or collectively, and any and all combinations of any two or more of said steps or features.
Furthermore, titles, headings, or the like are provided to enhance the reader's comprehension of this document, and should not be read as limiting the scope of the 25 present invention.
The entire disclosures of all applications, patents and publications, cited above and below, if any, are hereby incoiporated by reference. However, the reference to any applications, patents and publications in this specification is not, and should not be taken as, an 30 acknowledgment or any form of suggestion that they constitute valid prior art or form part of the common general knowledge in any country in the world. 76 Throughout this specification and any claims which follow, unless the context requires otherwise, the words "comprise", "comprising" and the like, are to be construed in an inclusive sense as opposed to an exclusive sense, that is to say, in the sense of "including, but not limited to".
SHR507647NZPR 301525122 77 REFERENCES Blott S, Kim JJ, Moisio S, Schmidt-Kiintzel A, Cornet A, Berzi P, Cambisano N, Ford C, Grisart B, Johnson D, Karim L, Simon P, Snell R, Spelman R, Wong J, Vilkki J, Georges M, Farnir F, Coppieters W. (2003) Molecular dissection of a quantitative trait locus: a 5 phenylalanine-to-tyrosine substitution in the transmembrane domain of the bovine growth hormone receptor is associated with a major effect on milk yield and composition. Genetics 163:253-266.
Buchanan FC, Thue TD, Yu P, Winkelman-Sim DC. (2005) Single nucleotide 10 polymorphisms in the corticotrophin-releasing hormone and pro-opiomelancortin genes are associated with growth and carcass yield in beef cattle. Anim Genet. 36:127-131.
Churchill GA, Doerge RW. (1994) Empirical threshold values for quantitative trait mapping. Genetics 138:963-971.
Coppieters W, Riquet J, Arranz JJ, Berzi P, Cambisano N, Grisart B, Karim L, Marcq F, Moreau L, Nezer C, Simon P, Vanmanshoven P, Wagenaar D, Georges M. (1998a) A QTL with major effect on milk yield and composition maps to bovine chromosome 14. Mamm Genome 9:540-544.
Coppieters, W.; Kvasz, A.; Arranz, J.J.; Grisart, B.; Farnir, F.; Mackinnon, M. & Georges, M. (1998b). A rank-based non parametric method to map QTL in outbred half-sib pedigrees: application to milk production in a grand-daughter design. Genetics 149:1547-1555.
Druet T, Georges M (2010) An efficient haplotype-based approach for combined linkage and linkage disequilibrium QTL mapping using Hidden Markov Models. Genetics 184:789-798.
Ge B et al. (2005) Survey of allelic expression using EST mining. Genome Res. 15:1584-1591.
SHR507647NZPR 301525122 78 Georges M (2007) Mapping, fine-mapping and molecular dissection of Quantitatve Trait Loci in domestic animals. Annu. Rev. Genomics Hum. Genet. 8:131-62.
Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, 5 Zusmanovich P, Sulem P, Thorlacius S, Gylfason A, Steinberg S, Helgadottir A, Ingason A, Steinthorsdottir V, Olafsdottir EJ, Olafsdottir GH, Jonsson T, Borch-Johnsen K, Hansen T, Andersen G, Jorgensen T, Pedersen O, Aben KK, Witjes JA, Swinkels DW, den Heijer M, Franke B, Verbeek AL, Becker DM, Yanek LR, Becker LC, Tryggvadottir L, Rafiiar T, Gulcher J, Kiemeney LA, Kong A, Thorsteinsdottir U, Stefansson K. (2008) Many 10 sequence variants affecting diversity of adult human height. Nat Genet. 40:609-615.
Haley CS, Knott SA, Elsen JM. (1994) Mapping quantitative trait loci in crosses between outbred lines using least squares. Genetics 136:1195-1207.
Hensen K, Braem C, Declercq J, Van Dyck F, Dewerchin M, Fiette L, Denef C, Van de Ven WJ.(2004) Targeted disruption of the murine Plagl proto-oncogene causes growth retardation and reduced fertility. Dev Growth Differ.46:459-470.
Kiefer BL et al. (2002) Exploring the opioid system by gene knock-out. Prog Neurobiol. 20 66:285-306.
Kneeland J, Li C, Basarab J, Snelling WM, Benkel B, Murdoch B, Hansen C, Moore SS. (2004) Identification and fine mapping of quantitative trait loci for growth traits on bovine chromosomes 2, 6, 14, 19, 21, and 23 within one commercial line of Bos taurus. J Anim 25 Sci. 82:3405-3414.
Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, Eyheramendy S, Voight BF, Butler JL, Guiducci C, Illig T, Hackett R, Heid IM, Jacobs KB, Lyssenko V, Uda M; Diabetes Genetics Initiative; FUSION; KORA; Prostate, Lung Colorectal and 30 Ovarian Cancer Screening Trial; Nurses' Health Study; SardiNIA, Boehnke M, Chanock SJ, Groop LC, Hu FB, Isomaa B, Kraft P, Peltonen L, Salomaa V, Schlessinger D, Hunter DJ, Hayes RB, Abecasis GR, Wichmann HE, Mohlke KL, Hirschhorn JN.(2008) SHR507647NZPR 301525122 79 Identification of ten loci associated with height highlights new biological pathways in human growth._Nat Genet. 40:584-591.
Lynch M, Walsh B (1997) Genetics and analysis of quantitative traits. Sunderland 5 Massachusetts: Sinauer associates, Inc.
McGowan KA et al. (2008) Ribosomal mutations cause p53-mediated dark skin and pleiotropic effects. Nature Genetics 40: 963-970.
Mizoshita K, Watanabe T, Hayashi H, Kubota C, Yamakuchi H, Todoroki J, Sugimoto Y. (2004) Quantitative trait loci analysis for growth and carcass traits in a half-sib family of purebred Japanese Black (Wagyu) cattle. J Anim Sci. 82:3415-3420.
Mizoshita K, Takano A, Watanabe T, Takasuga A, Sugimoto Y. (2005) Identification of a 15 1.1-Mb region for a carcass weight QTL on bovine Chromosome 14. Mamm Genome 16:532-537.
Nkrumah JD, Sherman EL, Li C, Marques E, Crews DH Jr, Bartusiak R, Murdoch B, Wang Z, Basarab JA, Moore SS. (2007) Primary genome scan to identify putative 20 quantitative trait loci for feedlot growth rate, feed intake, and feed efficiency of beef cattle. J Anim Sci. 85:3170-3181.
Ross SA et al. (2000). Retinoids in embryonal development. Physiological Reviews 80: 1021-1054.
Sagata, N. (1997) What does Mos do in occytes and somatic cells? Bioessaysl9: 13-21. Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, et al. (2005) Evolutionary conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 15:1034-1050.
Spelman RJ; Miller FM, Hooper JD, Thielen M, Garrick DJ: Experimental design for QTL trial involving New Zealand Friesian and Jersey breeds, in: Proc Assoc Assoc Advmt Breed Genet 2001,14:393-396.
SHR507647NZPR 301525122 80 Takasuga A, Watanabe T, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujikawa A, Chiba K, Kobayashi N, Tatsuda K, Oe T, Furukawa-Kuroiwa M, Nishimura-Abe A, Fujita T, Inoue K, Mizoshita K, Ogino A, Sugimoto Y. (2007) Identification of 5 bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping. Mamm Genome 18:125-136.
Trinklein ND, Aldred SF, Hartman SJ, Schroeder DI, Otillar RP, Myers RM. (2004) An abundance of bidirectional promoters in the human genome. Genome Res. 14:62-66. 10 Vandesompele et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging multiple internal control genes. Genome Biol. 3:research0034.
Van Dyck F, Declercq J, Braem CV, Van de Ven WJ. (2007) PLAGl, the prototype of the PLAG gene family: versatility in tumour development. Int J Oncol. 30:765-774.
Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, Lawley CT, Haudenschild CD, Moore SS, Warren WC, Sonstegard TS. (2008) SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat Methods 5:247-252.
Visscher PM, Thompson R, Haley CS (1996) Confidence intervals in QTL mapping by bootstrapping. Genetics 143:1013-1020.
Voz ML, Mathys J, Hensen K, Pendeville H, Van valckenborgh I, Van Huffel C, Chavez 25 M, Van damme B, De Moor B, Moreau Y, Van de Ven WJM (2004) Microarray screeing for target genes of the proto-oncogene PLAGl. Oncogene 23: 179-191.
Weedon, M.N.; Frayling, T.M. (2008). Reaching new heights: insights into the genetics of human stature. Trends in Genetics 24: 595-603.
Wibowo et al (2008). Genome assembly anchored QTL map of bovine chromosome 14. Int. J. Biol. Sci. 4(6): 406-414. 507647NZPR 302016061 Received at IPONZ on 24.02.2012 81

Claims (41)

CLAIMS:
1. A method for inferring the size potential of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the size potential of the animal and/or its offspring.
2. A method for inferring the growth rate potential and/or liveweight potential of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the growth rate potential and/or the liveweight potential of the animal and/or its offspring.
3. A method for selecting or rejecting an animal from the Bovidae family, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos Taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the size potential of the animal and/or its offspring. 507647NZPR 302016061 Received at IPONZ on 24.02.2012 82
4. A method as claimed in claim 3 wherein the method comprises at least the steps of: a) analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos Taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof; and, b) selecting or rejecting an animal based on the nucleotide present at the one or more positions.
5. A method for estimating the worth of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present infers of the size potential of the animal and/or its offspring.
6. A method for estimating the worth of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064,23232324, 23235239, 23243588, 23243734,23264810,23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present infers of the growth rate potential and/or the liveweight potential of the animal and/or its offspring. 5G7647NZPR 302086289 83 received 16 april 2012.
7. A method as claimed in any one of claims 1 to 6 wherein the presence of one or more of 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 232I5765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a larger size, a higher growth rate and/or a higher liveweight for the animal and/or its offspring.
8. A method as claimed in any one of claims 1 to 6 wherein the presence of one or more of22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C infers a smaller size, a lower growth rate and/or a lower liveweight for the animal and/or its offspring.
9. A method of identifying an animal from the Bovidae family which is more likely or less likely to produce one or more desirable trait selected from size, growth rate and liveweight, the method comprising at least the step of analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380,23186648, 23187238, 23205280, 23211989,23215765, 23218552, 23228064,23232324, 23235239, 23243588, 23243734,23264810,23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the nucleotide present at the one or more positions infers the animal is more likely or less likely to produce the one or more desirable trait.
10. A method as claimed in claim 9 wherein the method comprises at least the steps of: a) analysing a nucleic acid from the animal to determine the nucleotide present at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238,23205280, 23211989, 23215765,23218552,23228064,23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471 or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof; and, 507647NZPR 302016061 Received at IPONZ on 24.02.2012 84 b) identifying whether the animal is more likely or less likely to produce one or more desirable trait on the basis of the nucleotide present at the one or more positions.
11. A method as claimed in claim 9 or 10 wherein animals which are more likely to have a smaller size, lower growth rate and/or lower liveweight are identified on the basis of the presence of one or more of 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 232 6 4854A, 23272034G, and 23273471C
12. A method as claimed in claim 9 or 10 wherein animals which are more likely to have a larger size, a higher growth rate and/or a higher liveweight are identified on the basis of the presence of one or more 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(000)11, 23264854G, 23272034A, and 23273471T.
13. A method for identifying an animal from the Bovidae family having a higher growth rate potential and/or higher liveweight potential, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide sequence at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989,23215765,23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the presence of one or more of 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T infers a higher growth rate and/or higher liveweight potential. Received at IPONZ on 24.02.2012 50764 7NZPR or 302016061 O-5
14. A method for identifying an animal from the Bovidae family having a lower liveweight potential and/or a lower growth rate potential, the method comprising at least the step of analysing a nucleic acid from said animal to determine the nucleotide sequence at one or more positions corresponding to the following positions on chromosome 14 of Bos taurus; 22986260, 23186380, 23186648, 23187238, 23205280, 23211989, 23215765, 23218552, 23228064, 23232324, 23235239, 23243588, 23243734, 23264810, 23264854, 23272034, and 23273471, or to determine the nucleotide sequence of one or more genetic marker in linkage disequilibrium with one or more thereof, wherein the presence of one or more of 22986260C, 23186380C, 23186648G, 2318723 8G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and 23273471C infers a lower liveweight and/or lower growth rate potential.
15. A method for breeding animals from the Bovidae family to produce offspring, which comprises selecting at least a first animal that has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T wherein the position of the genetic marker is read in relation to chromosome 14 of Bos taurus, and mating said first animal with a second animal to produce offspring.
16. A method as claimed in claim 15 wherein the second animal is selected on the basis it has one or more of the following genotypes 22986260T, 23186380T, 23186648A, 23187238A, 23205280A, 23211989(deletion), 23215765C, 23218552T, 23228064G, 23232324A, 23235239C, 23243588 (deletion), 23243734A, 23264810(CCG)11, 23264854G, 23272034A, and 23273471T.
17. A method for breeding animals from the Bovidae family to produce offspring, which comprises selecting at least a first animal that has one or more of the following genotypes 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and Received at IPONZ on 24.02.2012 507647NZPR Q/; 502016061 oO 23273471C wherein the position of the genetic marker is read in relation to chromosome 14 of Bos taurus, and mating said first animal with a second animal to produce offspring.
18. A method as claimed in claim 17 wherein the second animal is selected on the basis it has one or more of the following genotypes 22986260C, 23186380C, 23186648G, 23187238G, 23205280G, 23211989G, 23215765G, 23218552C, 23228064T, 23232324G, 23235239G, 23243588TT, 23243734G, 23264810(CCG)9, 23264854A, 23272034G, and23273471C.
19. A method as claimed in any one of claims 1 to 18 wherein analysis of the nucleotide sequence of the one or more genetic markers occurs using one or more of: polymerase chain reaction (PCR); gel electrophoresis; Southern blotting; nucleic acid sequencing; restriction fragment length polymorphism (RFLP); single-strand confirmation polymphism (SSCP); LCR (ligase chain reaction); denaturing gradient gel electrophoresis (DGGE); allele-specific oligonucleotides (ASOs); proteins which recognize nucleic acid mismatches; RNAse protection; oligonucleotide array hybridisation; denaturing HPLC (dHPLC); high resolution melting (HRM); and, matrix-assisted laser desorption/ionisation time-of-flight mass spectroscopy (MALDI-TOF MS).
20. A method for inferring the size potential of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal, wherein the level infers the size potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
21. A method for inferring the growth rate and/or the liveweight potential of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or Received at IPONZ on 24.02.2012 50764 7NZPR 07 302016061 0/ more nucleic acids encoding the one or more proteins, precursors, isoforms, or fragments in said animal, wherein the level infers the growth rate potential and/or the liveweight potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
22. A method for selecting or rejecting an animal from the Bovidae family, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, and selecting or rejecting an animal based on the level of one or more of said proteins, precursors, isoforms or fragments thereof, and/or nucleic acids encoding same, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
23. The method of claim 22 wherein the method comprises comparing the level of expression of one or more of said proteins, precursors, isoforms, fragments, and/or nucleic acids against a standard, wherein a difference in the level of expression in the sample compared to the standard infers the size potential, growth rate potential and/or liveweight potential of the animal and/or its offspring.
24. A method for estimating the worth of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, wherein the level infers the size potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CIICIID7. RDHE2, SDR16C6, and PENK.
25. A method for estimating the worth of an animal from the Bovidae family and/or its offspring, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or 507647NZPR 302086289 88 received 16 april 2012. more proteins, precursors, isoforms or fragments in said animal, wherein the level infers the growth rate potential and/or the liveweight potential of the animal and/or its offspring, and wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK,
26. A method of identifying an animal from the Bovidae family which is more or less likely to produce one or more desirable trait selected from size, growth rate and liveweight, the method comprising at least the step of observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK, wherein the level infers the animal is more likely or less likely to produce the one or more desirable trait.
27. A method as claimed in claim 26 wherein the method comprises at least the steps of: a) observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK; and, b) identifying whether the animal is more likely or less likely to produce one or more desirable trait on the basis of the level of the one or more protein, isoform, precursor, fragment or nucleic acid encoding any one thereof.
28. A method for identifying an animal from the Bovidae family having a higher growth rate potential and/or a higher liveweight potential, the method comprising at least the step of observing the level of expression of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, wherein an increase in the level of expression compared to a standard infers a higher growth rate potential and/or higher Received at IPONZ on 24.02.2012 507647NZPR OQ 302016061 oy liveweight potential, and wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
29. A method for identifying an animal from the Bovidae family having a lower growth rate potential and/or a lower liveweight potential, the method comprising at least the step of observing the level of expression of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in said animal, wherein a decrease or substantially no increase in the level of expression compared to a standard infers a lower growth rate potential and/or lower liveweight potential, and wherein the one or more protein is selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
30. A method for breeding animals from the Bovidae family to produce offspring, the method comprising observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in at least a first animal, selecting the first animal where there is an increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, and mating said first animal with a second animal to produce offspring, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
31. A method as claimed in claim 30 wherein the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in the second animal, and selecting the second animal where there is an increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more Received at IPONZ on 24.02.2012 507647NZPR q/\ 302016061 "U proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
32. A method for breeding animals from the Bovidae family to produce offspring, the method comprising observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in at least a first animal, selecting the first animal where there is a decrease or substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, and mating said first animal with a second animal to produce offspring, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
33. A method as claimed in claim 32 wherein the method further comprises observing the level of one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and/or one or more nucleic acids encoding the one or more proteins, precursors, isoforms or fragments in the second animal, and selecting the second animal where there is a decrease or substantially no increase in the expression of one or more of said proteins, precursors, isoforms, fragments and/or nucleic acids encoding same compared to a standard, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK.
34. A method as claimed in any one of claims 20 to 33 comprising at least the steps of: a) taking a sample from an animal; b) detecting one or more proteins, one or more precursors thereof, one or more isoforms thereof, one or more fragments thereof, and one or more nucleic acids encoding any one or more of the proteins, precursors, isoforms or fragments in the Received at IPONZ on 24.02.2012 507647NZPR Q1 302016061 sample, wherein the one or more proteins are selected from RPS20, MOS, PLAGl, CHCHD7, RDHE2, SDR16C6, and PENK; and, c) comparing the level of the one or more proteins, precursors, fragments, isoforms and/or nucleic acids against a standard.
35. A method as claimed in any one of claims 20 to 34 wherein at least an approximately 1.2 fold increase in the level of one or more of the proteins, a precursor thereof, an isoform thereof, a fragment thereof, and/or a nucleic acid encoding same, is indicative of a higher growth rate potential and/or a higher liveweight potential of the animal and/or its offspring.
36. A method as claimed in any one of claims 20 to 35 wherein the level of the one or more proteins, precursors, fragments, isoforms and nucleic acids encoding same is determined using an immunoassay, separation based on characteristics such as molecular weight and isoelectric point, gel electrophoresis, Western Blotting or mass spectroscopy.
37. A method as claimed in any one of claims 1 to 36 wherein the method further comprises analysis of one or more additional biological markers.
38. A method of any one of claims 1 to 37 wherein the animal is bovine.
39. A method as claimed in claim 38 wherein the bovine animal is Bos taurus or Bos indicus.
40. A method as claimed in claim 39 wherein the animal is chosen from the group consisting Jersey, Holstein, Friesian, Simmental or crossbred dairy cattle.
41. A method as claimed in any one of claims 1, 2, 3, 5, 6, 9,13,14,15,17, 20, 21, 22, 24,25, 26, 28, 29, 30 or 32 substantially as hereinbefore described with particular reference to any one of the examples, figures and/or sequence listings.
NZ57945610A 2010-08-31 2010-08-31 Biological markers on Bovine Chromosome 14 for size and uses therefor NZ579456A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
NZ57945610A NZ579456A (en) 2010-08-31 2010-08-31 Biological markers on Bovine Chromosome 14 for size and uses therefor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
NZ57945610A NZ579456A (en) 2010-08-31 2010-08-31 Biological markers on Bovine Chromosome 14 for size and uses therefor

Publications (1)

Publication Number Publication Date
NZ579456A true NZ579456A (en) 2012-05-25

Family

ID=46125079

Family Applications (1)

Application Number Title Priority Date Filing Date
NZ57945610A NZ579456A (en) 2010-08-31 2010-08-31 Biological markers on Bovine Chromosome 14 for size and uses therefor

Country Status (1)

Country Link
NZ (1) NZ579456A (en)

Similar Documents

Publication Publication Date Title
US20090269741A1 (en) Method for assessing traits selected from longissimus dorsi peak force, intramuscular fat, retail beef yield and net feed intake in bovine animals
EP3060682B1 (en) Genetic markers and uses therefor
US7947444B2 (en) Leptin promoter polymorphisms and uses thereof
US20120012065A1 (en) Systems and Methods for Improving Protein and Milk Production of Dairy Herds
EP4273268A9 (en) Genetic markers and uses therefor
WO2011028134A1 (en) Biological markers and uses therefor
AU2009275988B2 (en) A genetic marker test for Brachyspina and fertility in cattle
US20060275793A1 (en) Association between markers in the leptin gene and carcass traits in commercial feedlot steer and heifers
KR102235340B1 (en) SNP marker set for predicting growth traits of Korean native chicken and uses thereof
US20080096207A1 (en) Leptin and Growth Hormone Receptor Gene Markers Associated with Rearing, Carcass Traits and Productive Life in Cattle
KR101663404B1 (en) Single nucleotide polymorphism marker in TDRKH gene for diagnosis of meat quality in Hanwoo and method for diagnosis of meat quality in Hanwoo using same marker
AU2007286129A1 (en) Association of single nucleotide polymorphisms, dairy form and productive life
NZ579456A (en) Biological markers on Bovine Chromosome 14 for size and uses therefor
US20210062260A1 (en) Genetic markers and uses therefor
WO2004070055A1 (en) Dna markers for marbling
WO2019125188A1 (en) Genetic markers and uses therefor
EP2149611A1 (en) A genetic marker test for brachyspina and fertility in cattle
WO2014193247A1 (en) Dwarfism genetic markers and uses therefor

Legal Events

Date Code Title Description
S883 Correction of error according to section 88(3) (mistake in register caused on part of patentee or applicant)

Free format text: THE INVENTORS (72) HAVE BEEN CORRECTED

PSEA Patent sealed
RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 3 YEARS UNTIL 31 AUG 2017 BY BALDWINS INTELLECTUAL PROPERTY

Effective date: 20140818

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2018 BY BALDWINS INTELLECTUAL PROPERTY

Effective date: 20170831

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2019 BY BALDWINS INTELLECTUAL PROPERTY

Effective date: 20180828

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2020 BY CATALYST INTELLECTUAL PROPERTY

Effective date: 20190828

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2021 BY CATALYST INTELLECTUAL PROPERTY

Effective date: 20200828

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2022 BY CATALYST INTELLECTUAL PROPERTY

Effective date: 20210816

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2023 BY CATALYST INTELLECTUAL PROPERTY LIMITED

Effective date: 20220824

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2024 BY CATALYST INTELLECTUAL PROPERTY LIMITED

Effective date: 20230725

RENW Renewal (renewal fees accepted)

Free format text: PATENT RENEWED FOR 1 YEAR UNTIL 31 AUG 2025 BY CATALYST INTELLECTUAL PROPERTY LIMITED

Effective date: 20240717