Agricultural and Food Sience, Vol. 16 (2007): 170-176 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 170 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 171 Quantitative trait loci for udder conformation and other udder traits in Finnish Ayrshire cattle Nina F. Schulman, Sirja M. Viitala, Johanna H. Vilkki MTT Agrifood Research Finland, Biotechnology and Food Research, FI-31600 Jokioinen, Finland, e-mail: nina.schulman@mtt.fi Udder traits are important due to their correlation with clinical mastitis which causes major economic losses to the dairy farms. Chromosomal areas associated with udder conformation traits, milking speed and leakage could be used in breeding programs to improve both udder traits and mastitis resistance. Quantitative trait loci (QTL) mapping for udder traits was carried out on bovine chromosomes (BTA) 9, 11, 14, 18, 20, 23, and 29, where earlier studies have indicated QTL for mastitis. A granddaughter design with 12 Ayrshire sire families and 360 sons was used. The sires and sons were typed for 35 markers. The traits analysed were udder depth, fore udder attachment, central ligament, distance from udder to floor, body stature, fore teat length, udder balance, rear udder height, milking speed, and leakage. Associations between markers and traits were analysed with multiple marker regression. Five genome-wise significant QTL were detected: stature on BTA14 and 23, udder balance on BTA23, rear udder height on BTA11, and central ligament on BTA23. On BTA11 and 14 the suggested QTL positions for udder traits are at the same position as previ- ously detected QTL for mastitis and somatic cell count. Key-words: Quantitative trait loci, udder conformation, dairy cattle © Agricultural and Food Science Manuscript received January 2007 Introduction Dairy cattle breeding has traditionally focused on the improvement of production traits. Recently other traits such as udder health, fertility, and conforma- tion have become more important since mastitis and fertility disorders are increasing among high yielding cows causing major economic losses to the dairy industry. Among all culling cases in Finnish Ayrshire 22.1% was due to mastitis, 21.4% due to fertility disorders, and 7.8% due to poor udder conformation in 2006 (Jukka Pösö, personal com- munication, 2007). When the replacement costs are estimated to be close to 30% of the total variable A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 170 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 171 costs for milk production (Enroth 2006), clear sav- ings could be obtained with improvement of these traits. Udder conformation traits have been shown to be related to udder health (Rogers and Hargrove 1990, Rupp and Boichard 1999). For example cows with higher udders, closer teat placement, and shorter teats are less likely to have high somatic cell count (SCC) (Rogers and Hargrove 1990). Stronger cen- tral ligament, higher udder depth, higher fore udder attachment and better udder balance are associated with less clinical mastitis (Rupp and Boichard 1999, Lund et al. 1994). Good udder conformation is also important for an easy and uniform milking routine, which is especially needed with automatic milking. Milking speed and leakage are other traits of the ud- der that also are correlated with udder health. Fast milking speed is associated with increased SCC but the relationship with mastitis is unclear (Rupp and Boichard 1999). During the last 15 years several quantitative trait loci (QTL) mapping experiments have been under- taken (Khatkar et al. 2004). The final objective of these studies is to locate genes or linked markers that can be used in the breeding scheme to improve the economically important traits. Most of the de- tected QTL are affecting milk production traits and a few papers are dealing with other traits of dairy cattle. Especially the functional traits would benefit from marker-assisted selection because they are of- ten traits that have low heritability and /or are diffi- cult to measure and the genetic gain with traditional breeding schemes has been moderate. Udder con- formation QTL have earlier been detected mainly in Holstein cattle (Schrooten et al. 2000, Boichard et al. 2003, Hiendleder et al. 2003). Quantitative trait loci for several conformation traits seem to be clustering on the same chromosomes, especially BTA6, 19, and 28. However, the QTL are not very consistent between studies. Putative QTL for milk- ing speed have also been found on many chromo- somes but none of them are observed in more than one investigation (Schrooten et al. 2000, Boichard et al. 2003, Hiendleder et al. 2003). In a previous analysis with an Ayrshire cattle granddaughter design QTL for SCC and mastitis were detected (Schulman et al. 2004). The objective of the present study was to investigate if QTL for udder conformation, milking speed and leakage of milk are located on the same chromosomes where udder health QTL were earlier detected. Material and methods Phenotypes and animals The analysed traits were stature (ST), fore teat length (FTL), udder balance (UB), rear udder height (RUH), udder depth, distance between udder and floor (DUF), central ligament (CL), fore udder attachment, milk- ing speed (MS), and leakage. Udder balance means the balance of the udder from the side view, RUH is the distance of the udder to the vulva, udder depth is the relative distance of the udder from the tarsal joint, CL is the depth of the central ligament of the udder, and fore udder attachment is measuring the angle between the udder and the abdomen. Among the 10 traits analysed FTL, UB, RUH, udder depth, DUF, and CL are measuring different aspects of udder conformation. Stature is measuring the size of the cow, and it is related to the udder conforma- tion traits. For example a bigger cow tends to have longer distance from udder to floor. Milking speed and leakage are functional traits of the udder. The conformation traits were scored subjectively by field advisors of the breeding organization. The scoring is done on a linear scale except for distance between udder and floor which is measured directly as the distance in cm from the rear teats to the floor. Milking speed indicates how fast the cow can be milked and leakage the tendency to leak milk before milking. The scoring of these two traits is based on the judge- ment of the farmer. Heritabilities of the traits range form 0.07 to 0.61 (J. Pösö, Finnish Animal Breeding Association, personal communication, 2007) (Table 1). Estimated breeding values (EBV) from the spring evaluation 2004 were obtained from the Finnish Animal Breeding association. Estimated breeding value for conformation traits were computed using animal model, repeatability model BLUP. Estimated breeding value for milking speed and leakage were estimated using a BLUP sire model. A G R I C U L T U R A L A N D F O O D S C I E N C E Schulman, N. et al. Udder trait QTL in Finnish Ayrshire 172 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 173 A granddaughter design with 12 sire families was used in this study except for BTA20 where only 11 families were genotyped. These families are closely related to the families analysed in a previous whole genome scan (Viitala et al. 2003, Schulman et al. 2004). The total number of sons with breeding values for ST, FTL, UB, and RUH was 291 ranging form 18 to 40 sons per sire with an average of 24 sons. For the other traits the number of sons avail- able was 360 ranging from 21 to 48 sons with an av- erage of 30 sons. The number of daughters per son differed between the traits. For ST, FTL, UB, and RUH the median of the number of daughters was 40 and the mean 97, for udder depth, DUF, CL, and fore udder attachment the median was 34 and mean 80, for MS and leakage the median was 81 and the mean 98. The minimum number of daughters per son varied between different grandsire families from 6 to 37. Semen samples for DNA extraction were obtained from Finnish AI stations. Markers and genotypes Seven bovine chromosomes, BTA9, 11, 14, 18, 20, 23 and 29 were chosen for this study because QTL affecting mastitis and/or SCC have been detected on these chromosomes in previous studies (Boichard et al. 2003, Schulman et al. 2004, Holmberg and Andersson-Eklund 2004) and udder health traits are shown to be correlated with udder conforma- tion traits (Rupp and Boichard 1999). Two SNPs (growth hormone receptor and prolactin receptor) and 33 microsatellite markers were typed and link- age maps were computed using CRIMAP (Green et al. 1990). Number of markers per chromosome varied from 3 to 7 with an average spacing of 11 cM. The average number of informative markers per chromosome ranged from 2 to 6. The marker tables are available at http://www.mtt.fi/julkaisut/ cattleqtl. Methods of DNA extraction, PCR reaction protocols and electrophoresis are described by Vilkki et al. (1997) and Viitala et al. (2003, 2006). Statistical analysis The data was analysed using multiple linear regres- sion (Knott et al. 1996). The analysis was done both between and within families. Chromosomes were analysed separately. The model was: yij = ai + bixij + eij, where yij is the breeding value of bull j, who belongs to family i, ai is the polygenic effect for half-sib family i, bi is the allele substitution effect for a QTL within family i, xij is the conditional probability for bull j of inheriting the first haplotype from sire i, and eij is the residual. The breeding value was regressed on the probability at 1 cM intervals. The analysis was nested within families and the test statistic was an F-ratio. The breeding values were weighted for the number of daughters and the heritability of the trait. Significance thresholds and P-values for the F statistic were obtained by permutation, which was repeated 10,000 times for each trait and chromo- some separately (Churchill and Doerge 1994). Genome-wise P-values were obtained by Bonfer- roni correction Pgenome=1-(1-Pchromosome) 29, where 29 is the number of bovine autosomes. Results and discussion Ten udder traits and 7 chromosomes were analysed in this mapping study. Five genome-wise significant Trait Heritability Stature 0.61 Front teat length 0.38 Udder balance 0.22 Rear udder height 0.27 Udder depth 0.32 Distance udder to floor 0.42 Central ligament 0.22 Fore udder attachment 0.28 Milking speed 0.20 Leakage 0.07 Data obtained from FABA Breeding (Jukka Pösö, personal communication). Table 1. Heritabilities of the analysed traits. A G R I C U L T U R A L A N D F O O D S C I E N C E Schulman, N. et al. Udder trait QTL in Finnish Ayrshire 172 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 173 QTL were detected: ST on BTA14 and 23, UB on BTA23, RUH on BTA11, and CL on BTA23. In addition, eight 5% and five 1% chromosome-wise significant QTL are of a suggestive nature. Genetic correlations between udder conformation traits are reported to be only moderate (Perttula 2007). Regarding this, it was somewhat surprising that several of the QTL were clustering on BTA11, 14 and 23 (Table 2). No QTL were found for udder depth and fore udder attachment. Also, no QTL were detected on BTA18. The false discovery rate of the analysis was calculated to be 19.4% ((0.05 × 10 × 7)/18 ) which means that 1/5 or 3.5 of the QTL are found by chance alone. The detected QTL should be regarded as preliminary results which should be confirmed in another independent experiment (Weller et al. 1998). Quantitative trait loci for five of the analysed traits (FTL, UB, RUH, DUF, and CL) were locat- ed on BTA11 close to the markers BMS2569 and INRA177. The distance between these two markers was 25 cM. The QTL for RUH was genome-wise significant. There were 2 to 4 segregating families per trait on this chromosome. One family had sig- nificant F-values for FTL, UB, and RUH. The posi- tion of the highest value of the test statistic for this family was at the distal end of the chromosome for UB and RUH and in the middle of the chromosome for FTL. The size of the QTL effects of the different traits among the significant families ranged from 0.52 to 1.58 standard deviations of EBV. The larg- est effects were detected for RUH and smallest for UB. No udder conformation trait QTL have previ- ously been detected on BTA11. On BTA14 we found QTL for ST, FTL, RUH, and DUF. The QTL for ST was genome-wise sig- nificant. QTL for ST and DUF seem to be located close to BMS740 and QTL for FTL and RUH be- tween the markers BMS1747 and RM11. There were 2 to 3 families with significant QTL effects per trait on BTA14. The size of the QTL effects ranged from 0.68 to 1.55 standard deviations of EBV. The largest effects were detected for FTL and the smallest for ST. A QTL for ST has pre- viously been detected on BTA14 close to marker BMS1941 (Spelman et al. 1999) which is 19 cM from BMS740 on the MARC map (http://www. marc.usda.gov/genome/genome.html ). However, the confidence intervals are likely to be large in both studies and it may be the same QTL. Quan- titative trait loci for ST, FTL, UB, DUF, CL, and MS were located on BTA23. The QTL for ST, UB, and CL were genome-wise significant. There were 1 to 5 segregating families per trait on this BTA11 BTA14 BTA23 Trait Pchr 1 Pgen 2 Markers Pchr 1 Pgen 2 Markers Pchr 1 Pgen 2 Markers Stature <0.01 0.01 BMS740 <0.01 <0.01 BM1818- BM1443 Front teat length 0.03 ns BMS2569- INRA177 0.01 ns BMS1747- RM11 <0.01 ns BM1818- BM1443 Udder balance <0.01 BMS2569- INRA177 <0.01 0.01 RM185 Rear udder height <0.01 <0.01 BMS2569 <0.01 ns BMS1747- RM11 Distance udder to floor 0.05 ns BMS2569- INRA177 0.02 ns RM11- BMS740 0.01 ns RM185- BM1818 Central ligament 0.01 ns INRA177 <0.01 0.04 BM1818- BM1443 Milking speed 0.02 ns RM185 1 chromosome-wise significance level 2 genome-wise significance level Table 2. Significance levels and locations of udder trait quantitative trait loci. A G R I C U L T U R A L A N D F O O D S C I E N C E Schulman, N. et al. Udder trait QTL in Finnish Ayrshire 174 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 175 chromosome. The size of the QTL effects ranged from 0.02 to 0.38 standard deviations of EBV. In a previous study with related animals, a QTL for live weight was detected between markers BM1258 and BoLA DRBP1 and a QTL for MS was suggested at the marker BoLA-DRBP1 (Elo et al. 1999). Live weight is related to ST because larger cows tend to have a higher live weight. These may be the same QTL as detected in the present study for ST and MS because the marker map in this study did not extend to the region of the chromosome with the earlier findings and because confidence intervals are likely to be large. Also, Ashwell et al. (2005) found a chromosome-wise suggestive QTL for ST in the proximal end of the chromosome and Schrooten et al. (2000) found a chromosome-wise significant QTL for MS between the markers BM1258 and GBCYP21. Further, Hiendleder et al. (2003) found a suggestive QTL for teat length at the distal end of BTA23. Quantitative trait loci were also suggested on BTA9 for CL, on BTA20 for UB, and on BTA29 for leakage. All these putative QTL were located at the end of the linkage groups and they were significant at the chromosome-wise level. A QTL for UB has previously been suggested on BTA20 but closer to the proximal end of the chromosome (Boichard et al. 2003). As udder conformation traits are reported to be correlated with udder health traits (Rupp and Boi- chard 1999), it is of interest to see if QTL for both trait groups are detected on the same chromosome locations. When comparing the locations of the udder conformation QTL detected in the present study with QTL for udder health traits from a previ- ous granddaughter design in the same breed, some overlapping areas were found (Fig. 1). A QTL for mastitis was detected on BTA11 close to the mark- er BM716 (Schulman et al. 2004). This marker is close to (2 cM apart on the MARC map) BMS2569 where a QTL for RUH was detected in the present study. Also, QTL for FTL, UB, and DUF were sug- gested in the marker interval next to BMS2569. In the previous study a QTL for SCC was detected close to the marker HELMTT41 which is 43 cM from BMS2569 on the MARC map. On BTA14 a QTL for mastitis was previously detected between markers BMS1747 and RM011. This is the same area where a QTL for FTL was suggested in the present study and might indicate that the QTL that affects teat length indirectly affects also mastitis susceptability. High correlation between teat length and mastitis has been reported in the literature (Lund et al. 1994). A QTL for SCC was previ- ously detected on BTA14 at the marker BMS740 which is close to the same position where QTL for ST and DUF were found in the present study. On BTA23 a QTL for SCC was suggested close to the marker RM033 in the previous granddaugh- ter design, and in the present study, QTL for UB and MS were found close to marker RM185. The marker RM033 is closer to the proximal end of BTA23 on the MARC map than RM185 which is the first marker on the map of the present study. Therefore this same QTL may have an effect on UB and MS. In conclusion, genome-wise significant QTL were detected on 3 of the 7 analysed chromosomes: stature on BTA14 and 23, udder balance on BTA23, rear udder height on BTA11, and central ligament on BTA23. In addition 13 chromosome-wise sig- Fig. 1. Locations of quantitative trait loci (QTL) for ud- der traits and udder health traits on BTA11 and BTA14. The udder health QTL were detected in a previous whole genome scan of Finnish Ayrshire (Schulman et al. 2004). The abbreviations of the traits are CM = clinical mastitis, RUH = rear udder height, UB = udder balance, DUF = dis- tance udder to floor, FTL = front teat length, CL = central ligament, SCC = somatic cell count, ST = stature. A G R I C U L T U R A L A N D F O O D S C I E N C E Schulman, N. et al. Udder trait QTL in Finnish Ayrshire 174 A G R I C U L T U R A L A N D F O O D S C I E N C E Vol. 16 (2007): 170-176 175 nificant QTL were found. Several of the QTL were clustering on BTA11, 14 and 23. On BTA11 and 14 the suggested QTL positions for udder traits are overlapping with previously detected QTL for mas- titis and SCC. Fine mapping of these areas could lead to the finding of new tools to improve both udder conformation and health. Acknowledgements. We would like to thank FABA Breed- ing for providing the material for the study. Jukka Pösö from FABA Breeding is especially acknowledged for valu- able information regarding the studied traits. This project has been partly funded by the Ministry of Agriculture and Forestry in Finland (project 310145, Mapping and selec- tion of genes affecting productive age of dairy cows) and FABA Breeding. References Ashwell, M.S., Heyen, D.W., Weller, J.I., Ron, M., Sonste- gard, T.S.,Van Tassell, C.P. & Lewin, H.A. 2005. De- tection of quantitative trait loci influencing conforma- tion traits and calving ease in Holstein-Friesian cattle. Journal of Dairy Science 88: 4111–4119. Boichard, D., Grohs, C., Bourgeois, F., Cerqueira, F., Faugeras, R., Neau, A., Rupp, R., Amigues, Y., Boscher, M.Y. & Leveziel, H. 2003. Detection of genes influencing economic traits in three French dairy cattle breeds. Genetics Selection Evolution 35, 1: 77–101. Churchill, G.A. & Doerge, R.W. 1994. Empirical thresh- old values for quantitative trait mapping. Genetics 138: 963–971. Elo, K.T., Vilkki, J., de Koning, D.J., Velmala, R.J. & Mäki- Tanila, A.V. 1999. A quantitative trait locus for live weight maps to bovine chromosome 23. Mammalian Genome 10: 831–835. Enroth, A. 2006. Mallilaskelmia maataloudesta 2006. Pro Agria Maaseutukeskusten Liiton julkaisuja 1030. 46 p. Green, P., Falls, K. & Crooks, S. 1990. Documentation of CRIMAP, version 2.4. Washington University School of Medicine, St. Louis, MO. Hiendleder, S., Thomsen, H., Reinsch, N., Bennewitz, J., Leyhe-Horn, B., Looft, C., Xu, N., Medjugorac, I., Russ, I., Kühn, C., Brockmann, G.A., Blümel, J., Brenig, B., Reinhardt, F., Reents, R., Averdunk, G., Schwerin, M., Förster, M., Kalm, E. & Erhardt, G. 2003. Mapping of QTL for body conformation and behavior in cattle. Jour- nal of Heredity 94, 6: 496–506. Holmberg, M. & Andersson-Eklund, L. 2004. Quantitative trait loci affecting health traits in Swedish dairy cattle. Journal of Dairy Science 87: 2653–2659. Khatkar, M.S., Thomson, P.C., Tammen, I. & Raadsma, H.W. 2004. Quantitative trait loci mapping in dairy cat- tle: review and meta-analysis. Genetics Selection Ev- olution 36: 163–190. Knott, S.A., Elsen, J.M. & Haley, C.S. 1996. Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theoretical and Applied Genet- ics 93: 71–80. Lund, T., Miglior, F., Dekkers, J.C.M. & Burnside, E.B. 1994. Genetic relationships between clinical mas- titis, somatic cell count, and udder conformation in Danish Holsteins. Livestock Production Science 39: 243–251. Perttula, H. 2007. Ayrshire-ensikoiden maidontuotanto- ja utarerakenneominaisuuksien perinnölliset yhteydet. Pro gradu. Helsingin Yliopisto. 33. Rogers, G.W. & Hargrove, G.L. 1990. Correlations among linear type traits and somatic cell counts. Journal of Dairy Science 74: 1087–1091. Rupp, R. & Boichard, D. 1999. Genetic parameters for clinical mastitis, somatic cell score, production, udder type traits, and milking ease in first lactation Holsteins. Journal of Dairy Science 82: 2198–2204. Schrooten, C., Bovenhuis, H., Coppieters, W. & van Aren- donk, J. A. 2000. Whole genome scan to detect quanti- tative trait loci for conformation and functional traits in dairy cattle. Journal of Dairy Science 8: 795–806. Schulman, N.F., Viitala, S.M., de Koning, D.J., Virta, J., Mäki-Tanila, A. & Vilkki J.H. 2004. Quantitative trait loci for health traits in Finnish Ayrshire cattle. Journal of Dairy Science 87: 443–449. Spelman, R.J., Huisman, A.E., Singireddy, S.R., Coppi- eters, W., Arranz, J., Georges, M. & Garrick D.J. 1999. Short Communication: Quantitative trait loci analysis on 17 nonproduction traits in the New Zeeland dairy popu- lation. Journal of Dairy Science 82: 2514–2516. Viitala, S.M., Schulman, N.F., de Koning, D.J., Elo, K., Ki- nos, R., Virta, A., Virta, J.,Mäki-Tanila, A. & Vilkki, J.H. 2003. Quantitative trait loci affecting milk production traits in Finnish Ayrshire dairy cattle. Journal of Dairy Science 86: 1828–1836. Viitala, S., Szyda, J., Blott, S., Schulman, N., Lidauer, M., Mäki-Tanila, A., Georges, M. & Vilkki, J. 2006. The role of the bovine growth hormone receptor and prolactin re- ceptor genes in milk, fat and protein production in Finn- ish Ayrshire dairy cattle. Genetics 173: 2151–2164. Vilkki, H.J., de Koning, D.J., Elo, K., Velmala, R. & Mäki- Tanila, A. 1997. Multiple marker mapping of quantita- tive trait loci of Finnish dairy cattle by regression. Jour- nal of Dairy Science 80: 198–204. Weller, J.I., Song, J.Z., Heyen, D.W., Lewin, H.A. & Ron, M. 1998. A new approach to the problem of multiple comparisons in the genetic dissection of complex traits. Genetics 150: 1699–1706. A G R I C U L T U R A L A N D F O O D S C I E N C E Schulman, N. et al. Udder trait QTL in Finnish Ayrshire 176 SELOSTUS Utareen rakenteeseen ja käyttöominaisuuksiin vaikuttavia geenialueita suomalaisessa ayrshirepopulaatiossa Nina F. Schulman, Sirja M. Viitala ja Johanna H. Vilkki MTT Biotekniikka- ja elintarviketutkimus Lehmän utareen rakenne- ja käyttöominaisuudet ovat tärkeitä, koska ne vaikuttavat utareterveyteen ja hyvän muotoinen utare myös helpottaa lypsyä. Näiden omi- naisuuksien geenialueita voitaisiin käyttää jalostus- ohjelmissa sekä utareen rakenteen että utaretulehdus- alttiuden parantamiseksi. Niinpä tässä tutkimuksessa kartoitettiin naudan kromosomit 9, 11, 14, 18, 20, 23 ja 29, joista aikaisemmassa tutkimuksessa oli löydetty utareterveyteen vaikuttavia geenialueita. Aineistona oli 12 suomalaista ayrshireperhettä, yhteensä 372 yksilöä. Kultakin yksilöltä tyypitettiin 35 geenimerkkiä. Tutkitut ominaisuudet olivat: utareen muoto, etukiinnitys, keskiside, maavara, takakorke- us, etuvedinten pituus, tasapaino, takakiinnityksen korkeus, lypsettävyys ja vuoto. Ominaisuuksien ja geenimerkkien väliset yhteydet analysoitiin monen markkerin regressiomenetelmällä. Tutkimuksessa löydettiin viisi genomitasolla merkitsevää geenialuetta: takakiinnityksen korkeus kromosomissa 11, takakorkeus kromosomeissa 14 ja 23, sekä tasapaino ja keskiside kromosomissa 23. Lisäksi löydettiin 13 kromosomitasolla merkitsevää geenialu- etta. Kromosomeissa 11 ja 14 löydetyt geenialueet ovat samoilla kohdilla kuin aiemmassa kartoitustutkimuk- sessa löydetyt utareterveyteen vaikuttavat geenialueet. Hienokartoittamalla nämä alueet saataisiin mahdol- lisesti uusia keinoja, joilla edistää utarerakennetta ja -terveyttä. Quantitative trait loci for udder conformation andother udder traits in Finnish Ayrshire cattle Introduction Material and methods Results and discussion References SELOSTUS