This is to certify that the thesis entitled BULK SEGREGANT ANALYSIS FOR BLOOM TIME QTL IN SOUR CHERRY (Prunus cerasus L.) presented by Ann Marie Bond has been accepted towards fulfillment of the requirements for the MS. degree in Plant Breedirg and Genetics/Horticulture @9530?“ MaerofeW ignature 777 // 4009‘ 3‘! Date MSU is an All'innafive Action/Equal Opportunity Institution LIBRARY Michigan State University PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 6/01 c:/CIRC/DateOuetp65-p.15 BULK SEGREGANT ANALYSIS FOR BLOOM TINIE QTL 1N SOUR CHERRY (Prunus caasus L.) By Ann Marie Bond A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Horticulture and Genetics 2004 5 Li; , {{[Ellrllr I\ll..| fill IL. i’li ABSTRACT BULK SEGREGANT ANALYSIS FOR BLOOM TIME QTL IN SOUR CHERRY (Prunus cerasus L.) By Ann Marie Bond Spring freeze damage to sour cherry (Prunus cerasus L.) flower buds is the major limiting factor to sour cherry production in the US. Significant crop reductions from spring freezes occur approximately every three years. In 2002, the losses to spring freeze damage were particularly devastating with 95% of the sour cherry crop in Michigan destroyed. One breeding approach to minimize the potential for freeze injury to sour cherry flowers is to develop late-blooming cultivars that would have an increased chance of avoiding spring freeze damage. Fortunately there is extensive variation for bloom time in sour cherry with extremely late blooming gerrnplasm available. The goal of this project was to investigate the inheritance of bloom time in tetraploid sour cherry (2n = 4x = 32) using a QTL approach. The progeny population used was from a cross between ‘Balaton’ and the late blooming sour cherry cultivar ‘Surefire’. QTL discovery was done using a bulk segregant approach comparing late and early blooming progeny individuals using SSR and AFLP markers. Using this approach, an AFLP marker was identified that was significantly linked to a putative QTL, termed blm3, controlling bloom time. The late allele, contributed by the late blooming parent ‘Surefire’, explained 16.5% of the phenotype variance and delayed bloom time by approximately 15 degree days. Future characterization of this QTL and other QTL controlling bloom time may suggest ways to delay bloom in sour cherry and other Rosaceous species, thereby reducing the probability of spring freeze damage and crop loss. DEDICATION I would like to dedicate my thesis to the men and women of the United States Armed Forces who unselfishly risk their lives, home and abroad, to preserve the way of life that I have the liberty and opportunity to enjoy. iii ACKNOWLEDGEMENTS I am very grateful to my advrsor, Dr. Amy Iezzoni, for her guidance, financial support, and editorial assistance. I would like to thank my committee, Dr. Dechun Wang, Dr. Jim Hancock, and Dr Dave Douches, for their patience, understanding, and last minute advice. Also, I would like to extend countless thanks to members of both the Iezzoni lab and the Hancock lab -— your wonderful sense of humor, words of advice, and constant guidance were of immeasurable value. To my roommates Lee Ann Pramuk and Janelle Glady - without whom I’d never have gotten his far — thanks so much for strawberries and coconuts! I would also like to thank the 7th Michigan Volunteer Infantry Company B for their constant support and love — you have become my Michigan family and will forever have a special place in my heart. And finally, Mom, Dad, Sarah, and Seth — Thanks for all the constant support, patience, and love. You have helped me attain the seemingly impossible by believing in me when I needed it the most! Bluebirds and Dandelions to you always! iv TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES I. INTRODUCTION II. LITERATURE REVIEW. A. Origin and history of sour cherry B. Floral development and freeze susceptibility C Variation 1n bloom time ... D. Genetic analysis of bloom time in Prunus ... III. MATERIALS AND METHODS A. Plant material B. Evaluation of bloom time . . C. Determination of early and late blooming bulk populations D. DNA isolation E. SSR analysis F. AFLP analysis G. QTL analysis IV. RESULTS A. Variation 1n bloom time B. Selection of early and late blooming bulk populations C. SSRanalyses D. AFLP analyses E. QTL analysis of marker E-AAA/CGT265 with bloom time V. DISCUSSION VII. APPENDICES .. A. Abbreviated SAS input files B. Population data for 1999, 2000, 2001 C. Bulk selection D. Degree Day data E. SSR raw data F. AFLP raw data G. Statistical comparisons VI. REFERENCES macaw; 10 10 12 12 14 16 17 20 25 33 34 36 37 47 48 vii 14 29 33 49 LIST OF TABLES Table l. Prunus linkage group(s) and reference information for the twenty-eight SSR primer pairs used in this study - Linkage groups are numbered as in Joobeur et al. (2000) ... ... 11 Table 2. Selective nucleotide sequences used for AFLP analyses with the restriction enzymes Eco RI and Mse I. 13 Table 3. Bloom data code and degree days (°C) to reach 50% bloom for the nine progeny selected for the early bulk ( = no data, letters refer to bloom code as described in Materials and Methods). . .. . .. . .. . . 16 Table 4. Bloom data code and degree days (°C) to reach 50% bloom for the nine progeny selected for the late bulk (- = no data, letters refer to bloom code as described in Materials and Methods). ... 17 Table 5. SSR primer pairs, fragments observed, and percent polymorphic fragments for Balaton, Surefire, and the 18 individuals in the early and late bloom bulk populations" (see Appendix F for SSR genotyping data). Table 6. AFLP analysis of Balaton, Surefire, and the bulk progeny individuals utilizing primers Eco RI +ANN and Mse I + CNN. .. . 21 Table 7. AFLP analysis of Balaton, Surefire, and the bulk progeny individuals utilizing primers Eco RI +AN and Mse 1+ CNN. Table 8. Presence or absence of E-AAA/M-CGT265 in Balaton, Surefire and the bulk progeny individuals. (na = no amplification in this particular PCR reaction). Table 9. Chi Square Analysis of the segregation of marker E-AAA/M-CGT265 in the Balaton and Surefire population to a 1:1 and 1:2. (*Significance at P < 0.05 for deviation from expected ratio, as = not significant). 25 vi LIST OF TABLES (continued) Table 10. Analysis of variance for the association of marker state for E-AAA/M- CGT265 with heat unit accumulation to 50% bloom time among 151 progeny evaluated in 2002 and 2003. 27 Table l 1. Mean heat unit accumulation for 50% bloom for progeny separated by the presence or or absence of blm3. (a and b denote significant difference within columns at P < 0.05). ' 27 APPENDIX C. Table 1. Day and degree day comparisons for 50% bloom to identify early and late blooming individuals for the early and late bulks (- = no assignment). D. Table 1. Growing degree days data (°C) for all Balaton x Surefire progeny individuals. E. Table 1. Presence or absence data for tested SSR pairs including parents and progeny evaluated. 44 F. Table 1. Presence or absence for E-AAA/M-CGT265 fi'om all the parents and progeny evaluated. 47 G. Table 1. SAS output for “group” comparison where group a is absent for E-AAA/M- CGT265 and b is present (therefore representing blm3 and Blm3 respectively, includes 2002 and 2003). 48 vii LIST OF FIGURES Figure 1. Frequency distributions of 50% bloom for the progeny population in 2002 (A) and 2003 (B). Bloom dates for the parents Balaton and Surefire are identified by B and S, respectively. Degree days were calculated in Celsius with a base temp of 4.4°C. 15 Figure 2. Frequency distributions of the degree days needed for the progeny individuals to reach 50% bloom in 2002 (A) and 2003 (B). The degree days for the parents Balaton (B) and Surefire (S) are shown as inserted letters respectively. White bars indicate degree days for progeny within the early bulk and dark grey bars for the respective late bulk. ... 18 Figure 3. Segregation for E-AAA/M-CGT265 among the parents and bulk progeny in this particular PCR reaction. B= Balation, S=Surefire, Early bulk = 2(6), 2(61), 3(2), 3(24), 3(25), 3(46), 4(31), 4(34), Late bulk = 2(7), 3(26), 3(51), 4(22), 4(24), 4(45). 24 APPENDIX B. Figure 1. Frequency distributions of 50% Bloom for the Balaton x Surefire progeny population over the first 3 years of bloom: 1999 (A), 2000 (B), and 2001 (C). Bloom date for the parents Balaton (B) and Surefire (S) are shown as inserted letters correspondingly. 34 viii I. INTRODUCTION Spring freeze damage to sour cherry (Prunus cerasus L., 2n = 4x) flower buds is the major limiting factor to sour cherry production in the United States. (Ricks, 1992). Significant crop reductions from spring freeze events occur approximately every three years. In 2002, the losses to spring freeze damage were particularly devastating with ~95% of the Michigan sour cherry crop destroyed (Kleweno and Matthews, 2003). The genetic uniformity of sour cherry production in the U. S., a monoculture of the 400 year old cultivar Montmorency, also set the stage for such a disaster. Within the flower, the pistil is the most sensitive tissue, and in the vast majority of freeze events, the pistil tissue is killed while the other flower tissues remain un-injured (Dennis and Howell, 1974). For example, sour cheny leaves can cold acclimate to survive temperatures of -13 2°C, where as cherry pistils are killed at temperatures between -2 to -3°C (Owens et al., 2001). Initial investigations suggest that cherry pistils may be deficient in their ability to cold acclimate due to absence in CBFI activity (Owens et al., 2002). CBFI is a transcription factor that has been demonstrated to be required for cold acclimation in other plant species (Jaglo et al., 2001). One breeding approach to minimize the potential for freeze injury to sour cherry flowers is to develop late-blooming cultivars. Since flowers become increasingly sensitive to freeze injury as they develop in the spring (Ballard et al., 1971), late blooming cultivars have an increased chance of avoiding spring freeze damage. Fortunately there is extensive variation in sour cherry for bloom time, as well as available gerrnplasm that blooms significantly later than Montmorency (Iezzoni and Hamilton, 1985; Iezzoni and Mulinix, 1992). Previous studies in Prunus species identified Quantitative Trait Loci (QTL) for bloom time. In an inter-specific peach x almond cross, a major QTL for bloom time (called Lb) was identified on linkage group 4 (Arus et al., 1999; Ballester et al., 1998). In a peach cross, 6 QTL for flowering time were identified with a large QTL on Group S (Yamamoto et al., 2001). In the sour cherry cross Rheinische Schattenmorelle x Erdi Boterrno, two QTL for bloom time were identified (Wang et al., 2000). However, the alleles identified at these loci, blml and blm2, conferred early bloom. A second cherry cross, Balaton x Surefire, was made to develop a pseudo-testcross population for the identification of late blooming QTL allele. Surefire was chosen as a parent because it is one of the latest blooming sour cherry selections. Since bloom time in sour cherry has a high heritability (Wang et al., 2000) and QTL analyses have been successfirl in Prunus, we hypothesize that it should be possible to identify additional or orthologous QTL for bloom time in the Balaton x Surefire cross. Our objective was to identify QTL controlling bloom time and the alleles at these QTL that contribute to late bloom. QTL discovery was done using a bulk segregant approach (Mikas et al., 1996; Michelrnore et al., 1991) comparing late and early blooming progeny individuals with_Simple Sequence Repeat (SSR) and Amplified Fragment Length Polymorphism (AFLP) marker analysis. II. LITERATURE REVIEW A. Origin and histog of sour chem Cherries are members of the Rosaceae family, which as a whole, is economically important in the globe’s temperate regions for both their fi'uits and ornamentals. Cherries have been valued since ancient times as one of the first tree fi'uits to ripen in temperate growing regions. The two main groups of cherries, the diploid sweet (Prunus avium, 2n = 2x = 16) and tetraploid sour cherries (syn. tart cherries, 2n = 4x = 32), are reported to have originated in an area that includes Asia minor, Iran, Iraq, and Syria (Watkins, 1976). The first diploid Prunus species arose in central Asia: sweet, sour, and ground cherry Prunusfi-uricosa, (a low growing bush cherry native to Russia) were early derivatives of this ancestral Prunus. Sour cherry is believed to have arisen through natural hybridizations between ground cherry and sweet cherry (Olden and Nybom, 1968). B. FIoLaLdevelopment and freeze susceptibilig Cheny production is limited to temperate regions that experience moderately cold winter temperatures. The cherry tree requires a dormancy period each year that begins with defoliation in the fall. Spring growth starts when the dormancy or chilling requirement has been satisfied and temperatures increase sufiiciently in the spring. In colder regions, cherry production is limited by cold mid-winter temperatures. For example, temperatures below -3 0°C can result in wood or flower bud injury, particularly in the more cold-sensitive sweet cherries. However, in some chen'y growing regions, low-temperature damage to flower buds is the most important factor limiting yields. Vegetative bud burst and flowering occur in the spring in response to warm temperatures favorable for rapid plant growth (above 10°C). Individual cherry flowers are perfect, consisting of sepals, petals, stamens, and a pistil. The stamen is made up of a filament and a pollen-bearing anther where as the pistil is comprised of the stigmatic surface, the style and the ovary containing a pair of ovules. In Prunus, one ovule typically develops into the seed while the other ovule aborts very early. The ovary develops into the fleshy pericarp of the drupe fruit. Since flowers become increasingly sensitive to freeze injury as they develop in the spring (Ballard et al., 1971; Dennis and Howell, 1974), late blooming cultivars have an increased chance of avoiding spring fi'eeze damage. Therefore, late bloom is a common goal in Prunus breeding programs (Bailey and Hough, 1975). Difi‘erences in bloom time among selections are fi’equently compared using heat accumulation above a given threshold to reach a designated pheonological state (Baskerville, 1968; Baskerville and Emin, 1969). The resulting calculated “degree-days” is therefore a more usefirl comparative measure for bloom time than calendar days. Since all the flower buds on a tree do not open at the same time, bloom time is sometimes evaluated when 50 percent of the blossoms are estimated to be open, is. 50% bloom (Wang et al., 2000). C. Variation in bloom time Sour cherries offers a unique opportunity to study late bloom in Prunus, as sour cherry and its parental species, P. fruticosa, are the latest blooming commercially important Prunus species (Iezzoni et al., 1990). Presumably P. fi-uticosa evolved extremely late bloom time to avoid spring freeze damage in its native habitat in central Russia. Due to continued introgression between sour cherry and P. fruticosa, there is extensive variation in sour cherry for bloom time, including germplasm that blooms significantly later than Montmorency (Iezzoni and Hamilton, 1985; Iezzoni and Mulinix, 1992) D. Genetic analysis of bloom time in Prunus Various linkage maps have been created in Prunus species; however, most of the linkage maps have been constructed from diploid crosses among peach selections or between peach and almond (Arus et al., 2003; Ballester et al., 1998; Yamamoto et al., 2001). Recently, SSRs (syn. micro satellite markers) have been suggested to be a valuable marker resource for comparative mapping (Aranzana et al., 2001). For example, 109 SSRs that were originally developed for peach and cherry, are being mapped on an almond x peach F2 population. They tested the marker variability between 25 cultivars (14 peach and 11 nectarine) and found 24 of the 35 SSRs to be polymorphic and consequently could distinguish between all of the individuals in the test (Aranzana et al., 2001). Testolin (2000) did similar research studying genetic origin and fingerprinting in peach, and all were polymorphic. Of the 26 SSR primer pairs they tested, 17 exhibited Mendelian inheritance. SSRs fiom peach and cherry were also found to be useful to distinguish among sour cherry germplasm accessions (Cantini et al., 2001). Additionally, AFLPs have been used to increase marker coverage and gene tagging. Unfortunately publication of linkage maps from different researchers originally resulted in different numbers associated with corresponding linkage groups. For example, numbered consecutively, the nomenclature used by European researchers for chromosomes 1, 2, 3, 4, 5, 6, 7, and 8 (Aranzana et al., 2003), are numbered 5, 2, 4, 6, 8, 3, 1, and 3, in a Japanese peach map [Groups 8 and 6 are shown as one linkage group (3)] (Yamamoto et al., 2001). This complicates the comparative QTL analysis. Fortunately, since 2002, the nomenclature of the European groups, as first published by J oobeur et al. (2000), has been adopted as the oflicial nomenclature for Prunus. - Despite the preliminary stages of linkage mapping in diploid Prunus, QTL for bloom time have been reported independently by three groups (Ballester et al., 1998; Yamamoto et al., 2001; Wang et al., 2000). In an inter-specific peach x almond cross, a major QTL for bloom time (Lb) was identified on linkage group 4 (Arus et al., 1999; Ballester et al., 1998). In a peach cross, 6 QTL for flowering time were identified with a large QTL on the Japanese peach map, Group S (Yamamoto et al., 2001). In the sour cherry cross Rheinische Schattenmorelle x Erdi Botermo, two bloom time QTL were identified (Wang et al., 2000); however, the alleles blml and blm2 conferred early bloom. The objective of this study is to identify QTL controlling bloom time, but also to identify the alleles at these loci that may contribute the desired late bloom phenotype. Because bloom time in sour cherry has a high heritability (Wang et al., 2000) and QTL analyses have been successful in Prunus, we hypothesize it should be possible to identify additional or orthologous QTL for bloom time in a second sour cherry cross. This second cherry cross, Balaton x Surefire, was designed to assist in the identification of late blooming alleles as Surefire is one of the latest blooming sour cherry selections. Marker analyses of the progeny of this cross would initially consist of markers previously been demonstrated to be linked to published bloom time QTL (Arus, 2003; Bliss et al., 2002; Wang et al., 1998; Yamamoto et al., 2001). Most of these markers are SSRs. Since these markers have been previously mapped in other Prunus species, it should be possible to determine if marker order and distance have been conserved between sour cherry and the other Prunus species. Ifnecessary, genome coverage to look for additional QTL could be accomplished with AFLP analyses as was performed in peach to map a gene associated with nematode resistance (Lu et al., 1998). - QTL identification in a heterozygous polyploid crop such as sour cherry is more challenging than in a diploid. Only markers that meet the criteria of a single dose restriction fragment can be placed on the linkage map, and most loci are duplicated in a polyploid. This requires that only markers in simplex in one or both parents can be mapped. QTL alleles must also meet this criteria for linkage mapping. Because identifying QTL in a polyploid species can be quite a challenge, the strategy is to use a bulk segregant analysis (Michelmore et al., 1991), whereby the initial goal is to identify marker polymorphisms that differ between a group of late blooming and early blooming selections. Bulk segregant analysis has previously been used to identify markers associated with QTL for quantitative traits (Mansur et al., 1993; Mikas et al., 1996). The association of marker(s) with a QTL for bloom time would subsequently be tested by scoring all progeny for the marker of interest. 1]]. MATERIALS AND METHODS A. Plant material. The two sour cherry cultivars used as parents to produce the population for QTL analysis were ‘Balaton’ and ‘Surefire’. Balaton is a landrace Hungarian variety called Ujfehertoi Furtos. Surefire is a cultivar released from the New York State Agricultural Experiment Station, Cornell University, from a cross between Borchert Black Sour x NY 6935 (Richmorency x Schattenmorelle) (R. Andersen, pers. comm.) Surefire was chosen as the late blooming parent as it is one of the extremes for this trait. The Balaton x Surefire cross was made in 1996. One hundred ninety-seven (197) seedlings were planted at the Michigan State University Clarksville Horticultural Experimental Station (CHES), Clarksville, Michigan in the spring of 1998. The seedlings were planted in field number 27 e and spaced within a row 1.5 meters apart, with the rows 6.1 meters apart. The following progeny individuals were not included in this study as they are suspected to have resulted fi'om out-crossing or self-pollination [designated by row(tree)]: 1(65), 2(20), 2(34), 2(54), 3(06), 3(11), 3(35), 4(08), 4(23), 4(36), 4(39), 4(47) (A. Sebolt and N. Hauck, pers. comm.) With these individuals eliminated, the final population size was 185. The following progeny have not yet been genotyped for their S-alleles: 2(5), 2(10), 2(14), 2(21), 2(22), 2(28), 2(29), 2(30), 2(31), 2(46), 2(49), 2(60), 3(9), 3(10), 3(13), 3(23), 3(32), 3(43), 3(50), 3(51), 3(54), 3(56), 3(60), 4(1), 4(4), 4(5), 4(17), 4(18), 4(19), 4(20), 4(30), 4(33), 4(3 7), 4(50), 4(51), 4(52), 4(53), 4(54), 4(59), 4(60), 4(63), 4(64), 4(65), and 4(66). B. Evaluation of bloom time. Bloom time observations were made daily (at 10AM beginning mid April) to record the date when approximately 50% of the flowers had opened for each seedling tree. Time to bloom was expressed as degree days (DD) from January 1 with a base temperature of 4.4°C using hourly temperature readings collected at the Automated Weather Station located at CHES (Baskerville and Emin, 1969). Daily heat unit accumulation was calculated by summing the positive differences of hourly temperature readings minus 4.4°C, then dividing by 24. On the date of 50% bloom, heat unit accumulation was calculated to hour 10, the approximate time data was recorded. Bloom time was evaluated over five years (1999 - 2003). C. Determination of early and late blooming bulk populations. In order to divide the population into two extreme phenotypic groups to form the early and late bulks, the progeny were organized chronologically for each year using the date of 50% bloom for each individual. For example, in 2003, the earliest blooming progeny [individuals 3(24) and 3(25)] reached 50% bloom on April 28tll or 177.2 DD. This day was recorded as day “a”. April 2911. was recorded as day “b” and so on until day “c”. This was performed for each year of available bloom data. A table was made including alist of progeny with their bloom date for each year and the corresponding bloom day letter code. Progeny with consistent early blooming dates and late blooming dates were chosen to form the “early bulk” and “late bulk”, respectively. The individuals in the early and late bulks were screened using the SSR and AFLP markers described below. Any marker(s) that was found to be present in one bulk but absent in the other, was subsequently genotyped for the entire progeny population of 151. D. DNA isolation: Young leaves were collected from the parents and each progeny individual, placed immediately on dry ice, transported back to the laboratory, and placed in the -80°C freezer for 24 hours. The samples were then fi'eeze-dried for at least 48 hours and then frozen at -20°C for storage. DNA was isolated using the CTAB method described by Stockinger et al. (1996). E. SSR analysis: Twenty-eight SSR primer pairs from both cherry and peach with known Prunus linkage map positions were used to screen the parents and bulk progeny individuals (Aranzana et al., 2001; Cantini et al., 2001; Cipriani et al., 1999; Dirlewanger et al., 2002; Sosinksi et al., 2000; Testolin et al., 2000; Wunsch et al., 2002) (Table 1). In an attempt to target those regions where bloom time QTL had previously been reported in Prunus, 10, 5, and 3 SSR primer pairs that map to Prunus linkage groups 1, 4, and 2, were tested (respectively). Primer pairs were selected by proximity (within 20 cM) to the published QTL. PCR amplifications were performed in a Perkin Elmer Cetus DNA Thermocycler 480 following procedures described in Cantini et al. (2001). To test the success of the PCR reaction, PCR products were separated by electrophoresis in 1% agarose gels with 0.1% ethidium bromide. For genotyping, the PCR products were mixed with 7.6 ul of fonnamide dye (10 ml 98% forrnamide, 200 ul 0.5mM EDTA, 10 mg 10 Table 1. Prunus linkage group(s) and reference information for the twenty-eight SSR primer pairs used in this study — Linkage groups are numbered as in Joobeur et al. (2000). Primer Pair Linkage Group (5) Reference BPPCT008 G6 Dirlewanger et al. (2002) BPPCTOI 5 G4 Dirlewanger et al. (2002) BPPCT021 Gl Dirlewanger et al. (2002) BPPCT023 G4 Dirlewanger et al. (2002) BPPCT038 G5 Dirlewanger et al. (2002) BPPCT04O G4 Dirlewanger et al. (2002) CPPCT002 G3 Aranzana et al. (2001) CPPCT003 G1 and G4 Aranzana et al. (2001) CPPCT026 G1 Aranzana et al. (2001) CPPCT027 G1 Aranzana et al. (2001) CPPCTO34 G1 Aranzana et al. (2001) PceGA25 G5 Cantini et a1. (2001) PceGA59 G1 Cantini et al. (2001) PceGA34 G2 A. Iezzoni (pers. comm.) pchgms3 G1 Sosinski et al. (2000)y pchgmsS G6 Sosinski et al. (2000)y PMSZ G7 Cantini et al. (2001) PMS3 G4 Cantini et al. (2001) PMS67 G1 Cantini et al. (2001) UDP96—008 G3 Cipriani et al. £999) 2 UDP96-018 G1 Cipriani (1999z UDP97-403 G3 Cipriani (1999) z UDP98-022 G1 Testolin et al. (2000) y UDP98-025 G2 Testolin et al. (2000)y UDP98-405 G7 Cipriani (1999)z UDP98-409 G8 @riani (1999)z UDP98-411 G2 Testolin et al. (2000) y UDP98-412 G6 Testolin et al. (2000) ’ 2 Further used in Prunus persica and P. avium by Testolin et a1. (2000), Wunsch et al. (2002). y Further used in P. persica L. and P. avium L. by Wunsch et al. (2002). 11 bromophenol blue, 10 mg Xylene Cyanol), heated at 95°C for 5 minutes and immediately placed on ice. Five microliters of each sample were loaded on a denaturing 6% polyacylamide gel. The samples were then electrophoresed at 80 watts for 2.5 hours and stained using the Silver Staining Protocol from the Promega Technical Manual. (Promega Corporation; Madison, Wisconsin) F. AFLP analysis: DNA digestion, adaptor ligation, pre-selective amplification and selective amplifications were carried out according to standard procedures described by Vos et al. (1995). Primer combinations included three selective bases for one primer (Mse I) and both two and three selective bases for the other (Eco RI) (Table 2). PCR products were first checked for successful amplification by electrophoresis in 1% agarose gels stained with 0.1% ethidium bromide. The PCR products were then prepared for polyacrylamide gel electrophoresis by mixing the reaction products with 7.6 ul of formamide dye (10 ml 98% formamide, 200 111 0.5mM EDTA, 10 mg bromophenol blue, 10 mg Xylene Cyanol), heated at 95°C for 5 minutes and immediately placed on ice. Five microliters of each sample were loaded on a denaturing 6% polyacrylamide gel. The samples were then electrophoresed at 80 watts for 2.5 hours and stained using the silver staining protocol according to Promega (Promega Corporation; Madison, Wisconsin). G. QTL analysis: A Nested Factor Design was used for a Single Marker QTL Analysis to determine if any marker is associated with bloom time. The linear additive model for the Analysis 12 of Variance is presented below. The two years represented the replications for this model. Computations were done using SAS (SAS Institute, Inc., 1999). Input files can be found in Appendix A. Yijkl = 11 + 8T00P1+ yearj + tredgroum + year’SIOUPij + $in where Yijk] = date of bloom in degree days 11 = grand mean _ group; = marker present/absent by genotype yearj = replications, 2002 and 2003 tree(group)k = progeny nested within the group year*groupij = interaction $in = years Table 2. Selective nucleotide sequences used for AFLP analyses with the restriction enzymes Eco RI and Mse I. Two base pairs Three base pairs EcoRI+AA/MseI+CAA EcoRI+AAA/MseI+CGT EcoRI+AA/MseI+CAC EcoRI+AGG/MseI+CAC EcoRI+AC/MseI+CAA EcoRI+ATA/MseI+CAC EcoRI+AC/MseI+CAC EcoRI+ATA/MseI+CCG EcoRI+AC/MseI+CAG EcoRI+ATA/MseI+CCT EcoRI+AG/MseI+CAA EcoRI+AT/MseI+CAC EcoRI+AG/MseI+CAC EcoRI+AT/MseI+CAA EcoRI+AT/MseI+CAC 13 IV. RESULTS Maligtion in bloom time Bloom time was recorded for all blooming progeny in the Balaton x Surefire population over five years (A. Iezzoni, unpublished data) (Appendix E). The number of blooming progeny reached a maximum of 175 individuals in 2002. In 1999, only 33 of the 3 year old seedlings bloomed followed by 85 individuals in 2000, 143 individuals in 2001, 175 individuals in 2002 and finally, 170 individuals in 2003. This trend likely represents the differences in the length of juvenility period among the seedlings with all the seedlings flowering by age seven. To maximize the number of progeny evaluated - only year 2002 and 2003 bloom data was used in the Single Marker Analysis (Appendix A for SAS input code). Bloom data was converted fi'om calendar days to degree days in order to compare bloom time over years on a heat unit basis. The progeny distribution for bloom time spanned 90.9 heat units in 2002 and 109.1 heat units in 2003 (Figure 1). The heat units accumulated for 50% bloom differed significantly between the two years; however, the year x progeny interaction was not significantly different over the two years. Balaton reached 50% bloom earlier than Surefire, exhibiting differences of 44.5 degree days in 2002 and 40.9 degree days in 2003. The progeny exhibited transgressive segregation for both early and late bloom time (Figure 1). In 2002, 28% of the progeny bloomed earlier than the early blooming Balaton parent and 5% of the progeny bloomed later than the late blooming Surefire parent. In 2003, 34% of the progeny bloomed earlier than Balaton and 6% of the progeny l4 2002 Bloom Distribution Number of Trees In Bloom Number of Trees In Bloom Figure 1. Frequency distributions of 50% bloom for the progeny population in 2002 (A) and 2003 (B). Bloom dates for the parents Balaton and Surefire are identified by B and S, respectively. Degree days were calculated in Celsius with a base temp of 4.4°C. 15 bloomed later than Surefire. These early and late blooming progeny individuals provided the opportunity to select individuals for bulk segregant analysis. B. Selection of early and late blooming bulk populations In order to divide the population into two extreme phenotypic groups to form the early and late bulks, the progeny were organized chronologically for each year using the date of 50% bloom for each individual. The population was divided into two extreme groups (denoted early and late) and the bloom data was organized chronologically for each year, the date of 50% bloom for each individual. Nine progeny with consistent early blooming dates were chosen to be a part of the “early bulk” for marker analysis: 2(6), 2(61), 3(2), 3(24), 3(25), 3(46), 3(54), 4(31), and 4(34) (Table 3). The range in heat units to reach 50% bloom of the individuals chosen for the early bulk was 177.2 to 216.6 degree days in 2003. The progeny population mean was 232.9 degree days. These selected progeny had noticeably earlier bloom time in comparison to the population mean over all years. Table 3. Bloom data code and degree days (°C) to reach 50% bloom for the nine progeny selected for the early bulk (— = no data, letters refer to bloom code as described in Materials and Methods). Early bulk 1999 2000 2001 2002 2003 rank hours Rank hours rank hours rank hours rank hours 2(6) f 226.7 c 221.8 b 188.7 - 199.7 - 211.2 2(61) a 187.2 c 221.8 b 188.7 a 179.8 c 195.2 3(2) - - c 221.8 b 188.7 f 193.3 - 216.6 3(24) b 192.4 a 211.8 b 188.7 c 189.1 a 177.2 3(25) - - - - - - f 193.3 a 177.2 3(46) - - - - b 188.7 c 189.1 b 186.9 3(54) - - - - c 205.3 c 189.1 e 208.1 4(31) - - e 234.5 a 178.7 c 189.1 c 195.2 4(34) - - (1 228.1 b 188.7 f 193.3 a 177.2 l6 Similarly for late bloom, the last day of 50% bloom for the population in 2003 was May 13th or 286.3 DD. This day was recorded as day “2”, May 12th was recorded as “y” and so on until day “v”. Nine progeny with consistent late blooming dates were chosen to be part of the “late bulk” for marker analysis: 2(7), 2(39), 3(26), 3(51), 4(22), 4(24), 4(32), 4(45), and 4(66) (Table 4). The nine progeny individuals selected for each bulks represented 4.9% of the total population. The range in heat units to reach 50% bloom of the individuals chosen for the late bulk was 246.6 to 286.3 degree days in 2003. The progeny population mean was 232.9 degree days. Table 4. Bloom data code and degree days (°C) to reach 50% bloom for the nine progeny selected for the late bulk (- = no data, letters refer to bloom code as described in Materials and Methods). Late bulk 1999 2000 2001 2002 2003 rank hours Rank hours rank hours rank hours rank hours 2(7) v 284.5 w 289.7 - 239.3 x 263.3 w 266.6 2(39) - - 2 343.4 2 274.8 y 268.0 y 283.6 3(26) - - - - - - - - 2 286.3 3(51) - - - - x 255.9 y 268.0 - 246.6 4(22) 2 318.8 2 343.4 x 255.9 w 257.8 2 286.3 4(24) - - y 324 1 2 274.8 2 270.7 2 286.3 4(32) - - - - - 221.8 w 257.8 w 266.6 4(45) - - x 306.6 x 255.9 v 251.1 x 280.6 4(66) - - x 306.6 - - y 268.0 2 286.3 C. SSR analyses: The SSR genotypes for 28 primer pairs were determined for the two parents and, individually, for the 17 progeny in the early and late bulks. No fragments were amplified with five of the primer pairs (Table 5). Considering the remaining primer pairs tested, the number of fragments amplified from these sour cherry individuals ranged from 0 to 9 17 2002 Bloom Distribution s68 odow «snow «.53 engage... v. EN m.m..mm§%mma%a§ N NnN m ......,. .....e..as.z . - III. m m gawmmmmmemmfifiamfim Nd «N B .fl III . m Baggage...“ a mom .n III . a t B aaafieaae o com M. W B v.mo~ D m 0 _ 32 m m .l e B o m3 m 3 - 0 n . 0 - m we 2 «N? T ado. 33 - . - D a 8. T 5.5? fl- . «NC a n at 5 0 505m5050 mamfiw 3 3 2 1 1 :30. c 80. E: EBOEQOO-PDEE-dz O - 5.3.8 2 18 Degree Days at 50% Bloom Figure 2. Frequency distributions of the degree days needed for the progeny individuals to reach 50% bloom in 2002 (A) and 2003 (B). The degree days for the parents Balaton (B) and Surefire (S) are shown as inserted letters respectively. White bars indicate degree days for progeny within the early bulk and dark grey bars for the respective late bulk. No-35: as. «225.: .8888 .8882 2852:: can 85E .2 8.53.. a; 82823.5 8... - 4N: .NN: .8: .2: .Na 4 233.25 - N2 .2: .2: .N2 4 :33me 8d 1 MN ON: .N: N: m 2:13.25 23 u 3. 2: .2: .2: .2: 2: m 2:18.25 28 u 5 Na .8 .3 2: .2: .2: .3 N N82229: mg n a: 2: N: .2: .N: 4 2:18.25 - EN .NoN .02 .EN .EN .NON 0 223.53 - m2 .3 .02 .mN: a. ”8-8.25 2: u 2m 8: .2: .2: .8: .N: m 82.: coon 0:. m: .E .N: 62 N2 .3 e Nmsa - o: .8: N 2.23.: 8d u we 2: .3: .2: a: 2: .w: o 2.522 3.0 u 3 SN .NNN .N: 3: .8: .3: .2: N 0208.: - o: .2: N 358.: 8° 1 2... N2 .8: a: ...: 2: .N: 0 48.526 - 3 .3. .8 .8 4 $385 one u N: m: 5 N 085.50 owe n we om .2.N 58 SN .OmN .NVN c 88026 - 3.22.8:2:23:18: N «85.2mm - 3N 2... SN N 388% mg n E... S. 2... .Nfl .3 o: 4 288.5 Moo 1 2m m: 2... .ma ..a N: 2... .2: .2: .2: .2: .3 0 8220.75 ENEEOEbom AENQ 33 E 35 $5592.“ GEmmlomNo E 3:3 3:08th 3:08me coca—25:80 Eoocom oEEoEboa 33030 oEEOEocoE 32030 mo .02 SEE ASE mEnbocom Mmm com m 265%? 83 ... gown—Naom— e==n Eco—n 8e— uee been 2: E £52222: 2 05 can .oEoBm deem—am 8m manage—wet 0399523 .583 wee 602030 Eco—HES.“ .Eeq SEE mwm .m «Sun. 19 (Table 5). Nine SSR primer pairs produced only monomorphic fragments, while the other 14 primer pairs amplified fragments were polymorphic among the parents and progeny tested (Appendix F). None of the polymorphic SSR fragments were consistently present in individuals from one bulk and absent from individuals in the other bulk. Therefore, this SSR analysis did not reveal a marker that could be tested as a candidate for linkage to a bloom time QTL. D. AFLP analyses: Fifteen AFLP primer combinations were used to screen the parents and individuals from both the early and late bulks. All fragment sizes were scored for all the individuals tested. The monomorphic and polymorphic fragments were summarized to evaluate the feasibility of firture AFLP linkage mapping in this population (Table 6 and 7). Initially, primer combinations with three selective nucleotides for both Eco RI and Mse I were used, which exhibited an average of 21 bands per gel (averaging 30% polymorphism). Based upon Vilanova et al. (2003), two selective nucleotides were used with the Eco RI primer producing on average of 39 bands per gel (averaging 14% polymorphism). Only one fragment was found that was present in one parent, absent in the other parent and differentially present in the bulk progeny. The candidate marker of ~265 bp was identified with AFLP primers Eco R1 + AAA and Mse I + CGT, and therefore the marker name assigned to this fragment was E-AAA/M-CGT255 (Figure 3). This primer pair was initially screened with both parents, eight early bulk individuals, and six late 20 22“ A .23.. 2 es. .2 m can A vs... _ as. .oaN .RN .NNN 53 .2.N .NNN .EN you + _ a: :.o .1. NS .22“ .SN .ooN .8: .2: .2: .8: d: 62 N 2.5.. + 2 82 80 + H 8: 85 u we o8 .2 N .03 .NN ON... .8... .03 .EN N 2.5. + E 82 can A 22.3 _ Ea oNN.oNN.NNN .oNN.NeN.eeN.omN.oNN.N_N.ooN on... A EB _ 2:. .Ne .Ne .8: .2: .8: .N: .8: 05 + _ a: «N6 n 3» ...NN .NNN .BN .2; .2: ...: .2: .5 .2: .2: .NE .2: .2: .2: mm 2.5. + E 82 23 A 853 _ Ea o? A as... _ e5 .oaN .oWN .QN .2N 05 + _ a: Sac n 2R .on .mNN .2: .Nm: .2; .NN: 53 .SN .8: .8: .02 .42 £2 2 $9.. + 2 82 98 A 355 m 23 2H A 2.3 _ as . 22 .23 .aNN .mNN 50 + _ 2.: eNd n 2% .SN 68 .2N .2: .NNN .omN .2.N .NNN .BN .2: .o: 2 2.5. + 2 82 EmEBoEbom Emma 83 E one.” €qu 023 E SE 3538.: 532588 «:38.— mEoEmEm 0282.:an 3202.0 fleas—wet 02808288 83890 mo .02 SEE .220 + _ a: 23 22$ 2 8m mafia 2.32: .3333 282...: :2. 2: 23 .0525 2.923% 3%.... .35.. .e ......c. 21 on: A :33 w 3:: .32” v 3:: :3 A 3:3 n 63 v 3:: com A 3:3 2 .::_.3:_.o:_ 0? v 3:: own A 3:3 _ .3: .32 .NS .32 .Nm_ 62 63 <2 + E Sm 3?. A 3:3 3 3:: .02 v 3:: own A 3:3 2 32 v 3:: 63 A 3:3 _ .33 v 3:: com A 3:3 Z .33 v 3:: com A 3:3 _ .32 .ofi 3:— .m2 .2; £2 .3. <00 .0 _ 0.00N 00.0.30 5.05N N03 .30 0.05N .0.0030 0.000 00.5030 0.0.0 00.0.30 .030. 0.503 00.0030 N. .0N N0.5030 0.00N .0.0030 5. .0N 00.0030 0.05N 003030 >00 0?. 0.NON 00.0030 ..0NN N03030 0.0 .N .0. .030 0.NON 00.330 0.00N 00.0030 00E0>0 0.00N 003030 0.00 . N0. .030 .. .0N 3.00.00 0.NON 00.5N.00 0.0 .N 00.00.00 >00 .03 N .55 . 00.0N.00 0.05. N0.0N.00 5.05 . .0.0N.00 0. . .N 00.0N.00 N .50 . 00.5N.00 .00....00 0.00N 00.0.30 0.00N N0.. .30 0.000 0030.00 00 0 0.00N 00.0030 5.0 .N N03030 0.00N .0.0030 00 0 0.00N 00.0030 5.0 .N N03030 0. .NN .03030 0.00N 0030.00 00 0 5.00N 00.0030 0.00N N0.5030 0.030.. 30. 00 . 0 5.00N 00.0030 0.00N N0.5030 0.00N .0.0030 N0 0 0.00N 00.0.30 5.0 .N N03030 0. .NN .03030 0.00N 00. .030 .0 0 0.0.0 0.00N N0.0030 0. .NN .03030 0.00N 00.330 00 0 N00. 00.00.00 5.00. N0. .030 0.030.. 0: 00 0 0.00N 00.. .30 .. .0N N0.0030 0.00N .0.0030 0.000 00.0030 00 0 0.00N 00.. .30 N.NON N0.0030 0.00N 3.0030 50 0 0.0 .N 00.0030 5.00. N0. .030 0.00N 3.330 0. .NN 00.0N.00 0.00N 00.0030 00 0 5.0NN 003030 5.0 .N N03030 0. .NN .03030 0.05N 00.0030 00 0 0.0.0 N.NON N0.0030 0.0.0 5.00N 00.0030 0 .00N 00.0030 00 0 0.032. 0: N.NON N0.0030 0.030.. 0.. 00 0 0.00N 00.0030 N .NON N0.0030 0.030.. 0: N0 0 N00. 00.00.00 5.00. N0. .030 0.030.. 0.. .0 0 N00. 00.00.00 0.00N N0.00.00 0.030.. 0: 00 0 ..NON 00.330 0.00 . N0.00.00 0.00N 3.330 00 0 N. . .N 00.0030 5.0 .N N03030 0L032. 0: 00 0 5.00N 00.0030 N.NON N0.0030 0. .NN .03030 5.00N 00.0030 0.00N 003030 00 0 0.00N 00.. .30 .. .0N N0.0030 0.00N .0.0030 0.000 0030.00 00 0 000 .0 $00 000 .0 $00 000 .0 $00 000 500.0 000 500.0 00.... 30¢ 0 00N NO0N .00N 000N 000 . 43 APPENDIX E - - - - .8. .8. 2 a: ..8-8-8 8.58.. . .8. .8. ... .8. 8. .8.. 8. .8.. .8. .8.. .8. .8.. .88. .8. .88. .8. .88. .8. - - 8...-8.8 v8.0.0.0 8. 8. 88 8. 88.8.... +8.8... 88.8.... 88.8.... - - 2 8-8-8 88.0.8 .28. .8188. .8188. .8. - - 8...... 8..-8-8 8.8.0....0 8. 8. 8. .8. 8. .8. .8. .8. .8. .8. 88.8. 8...... ......8. 88.8. - - .8182 “...-8.8 880.8 - - - - 2 8 2. ...-8-8 .85.... 2 2 2 2 - - 8 83.8.8 895...... 8.. .8. .8. .8.. .... - - - - 2 .8.. ..o. a. 2-8-8 88.0.... 2 .... 2 2 - - 2 83.8.8 888...... v8.8. 8...... 2 v8.8. - - .88... 83.8.8 .85.... 8.. .8. .88. ..8. .8. .8. .8. .8. .88. .8. 8..... .8.... 8.8. 8. - - 8..... 8.8.8 8.8.0.... 8.. .... 8.. .... 8.. 8. .8.. .8. .8.. .8. .8.. 8...... 8. .8: .... .8. .8. .8. .818. .818. .8. .8. .818. .8. .8. 8.. 888-8 .88.... 8.88.... 8.88.... 8.88.8 8.88.... .88.... .... .... 88-8-8 885...... 0.... ...... .88.. 08.. 2828 8.8.8 888.. 2:2 0.2:. 89.55.... 00.80.85 .8000... 0.8 8.5.8.. 05.3.0... 8..... .000 0882 .0.. S... 8:08.... .o 3580... .. 0.0:. E. Table 1. (continued) m... . .0.». .N0. - - - - 8 8 «8-8-8 8.8-8.00 vo...o. 8...... 8. ...: .8. .8 .8. .... 2-8-8 - - - - .8. ... .88 .8 ... .8.. .8. .8.-8-8 88-8.00 .... .... .8. - - - - .8. .8. .8. .8. .818. 2-8-8 8..-.800 8. ..8. 8. ..8. .88. .8. .88. .8. 8. - - - - .8. ..8 .8. ..8 .... .8. 88.8.8.8 8.8-89.0: 8. .... .8.-8-8 - - - - 8 .8. .8. 8 88.88 88.88.00 .882. 8. .8. .88. 8. 8. .8. .8. .8. .8.: 8. .8. 8. .8. - - .8. .8. 8..-8-8 82.. .8. .8. 8... ...; ...-8-8 - - - - .... .8. .8. .8.. 8 .8.-8-8 .8.)... - - - - 8.. .88. 8.. 88. 8 88-8-... 8285 8. 8. 88. .8. .8. .8. .88. .8. 88. .8.. .8. .8.. .8.. .8.. .8. .8.. - - 8.. 8-8-8 898...... ...-<08. 8. 8. .8.. .8. 8. .8.. .8. .8. .8. 8. .8. 8. .8. .8. .8. 8. .8. .8. .8. .... .8.. .8. .... .8. .... .... .8.. .8. .... .... .8.. .8. .8. 8.-. .-.o 8.68.. .8... ...... .88.. 08.. 2888 :28... .88.. 8a.. 885. ...... be... 45 E. Table 1. (continued) .... .... .... .8. .... .8. .... - - - - .8. .8 .8. .8 .... .8. 888-8 .3885 .8. .8. .8. .8. .8. .88. - - - - .8.... 8..... .81.... ..-8-8 ...-8.00 ... 8.. - - - - .8. .8.. ... .... .8. 88.88 8.18.00 8.. 8. .8. 8. 8.. .8. .8. .8. 88.8.8 .8. ...: 8. .8. .8. ...: .8. .8. .8. .8. .8. .8. 8. ...: .88-8-8 8..-8.00 .8.... ...... .88.. 08.. 8828 8.8.8.. .23.. 25: 882:. ...... as... 46 APPENDIX F Table 1. Presence or absence for B-AAA/M-CGT265 fiom all the parents and progeny evaluated. L Balaton I absent I Surefire I presetfl 1(66) absent 2(55) absent 3(34) present 4(13) absent 2(04) absent 2(56) present 3(36) present 4(14) absent 2(05) absent 2(57) present 3(37) absent 4(15) present 2(06) absent 2(58) present 3(38) present 4(16) ‘ present 2(07) present 2(59) absent 3(3 9) absent 4(18) absent 2(09) present 2(60) absent 3(40) absent 4(19) absent 2(10) absent 2(61) absent 3(41) Jresent 4(21) absent 2(1 1) absent 2(62) absent 3(42) absent 4(22) present 2(12) present 2(63) absent 3(43) absent 4(24) present 203) absent 2(64) absent 3(45) present 4(25) absent 2(15) absent 2(65) absent 3(46) absent 4(26) absent 2(16) absent 2(66) absent 3(47) present 4(27) present 2(17) absent 3(01) absent 3(48) absent 4(28) present 2(18) present 3(02) absent 3(49) absent 4(29) absent 2(19) present 3(03) absent 3(50) absent 4(30) absent 2(22) absent 3(05) present 3(51) present 4(3 1) absent 2(24) absent 3(07) present 3(52) absent 4(33) Jresent 2(28) absent 3(08) present 3(53) absent 4(34) absent 2(29) absent 3(10) absent 3(54) absent 4(35) absent 2(30) absent 3(12) absent 3(55) absent 4(37) absent 2(3 1) absent 3(14) absent 3(57) absent 4(40) absent 2(32) present 3(15) present 3(58) present 4(41) absent 2(35) present 3(16) absent 3(59) present 4(42) absent 2(36) present 3(17) absent 3(60) absent 4(43) absent 2(37) present 3(19) absent 3(61) absent 4(44) absent 2(38) absent 3(20) present 3(62) absent 4(45) absent 2(40) absent 3(21) present 3(63) absent 4(46) absent 2(41) absent 3(22) absent 3(64) absent 4(49) present 2(42) absent 3(24) absent 3(65) present 4(50) absent 2(43) absent 3(25) absent 4(01) absent 4(51) absent 2(45) absent 3(26) present 4(02) absent 4(52) absent 2(46) present 3(27) present 4(03) absent 4(54) present 2(47) present 3(28) present 4(04) absent 4(55) absent 2(48) absent 3(29) absent 4(07) absent 4(56) absent 2(50) absent 3(30) absent 4(09) absent 4(58) present 2(51) absent 3( 3 1) absent 4(10) absent 4(59) absent 2(52) present 3(32) absent 4(11) present 4(61) absent 2(53) absent 3(33) absent 4(12) absent Balaton = absent Surefire = present Progeny = 106 absent and 45 present (Total = 151) 47 APPENDIX G Table 1. SAS output for “group” comparison where group a is absent for E-AAA/M- CGT265 and b is present (therefore representing blm3 and BIm3 respectively, includes 2002 and 2003). Least Sguares Means Effect Group Estimate Std Error DF t Value Pr >j t 1 group a 220.08 1.9483 141 112.96 < .0001 group b 235.21 3.0052 137 78.27 ' < .0001 48 VI. REFERENCES Agricultural Research Service (ARS) Web page (last updated 5/13/2004) http://www.arsusdggov/is/Aanichive/apr99/cher0499.htm Aranzana, M.J., Garcia-Mas, J ., Carbo, J., Arus, P. 2001. Development and variability analysis of microsatellite markers in peach. Plant Breeding. 120: 1-6. Aranzana, M.J., Pineda, A., Cosson, P., Dirlewanger, E., Ascasibar, J., Testolin, R., Abbott, A., King, 6., Iezzoni, A., Arus, P. 2003. A set of simple-sequence repeat (SSR) markers covering the Prunus genome. Theoretical and Applied Genetics. 106(5):819-825. Arus, P., Dirlewanger, E., Quarta, R., Tobutt, K.R., Ballester, J., Boskovic, R., Dettori, M.T., de Vicente, M.C., Jauregui, B., Joobeur, T., Russell, K., Verde, 1., and Viruel, M.A. Location of 20 major genes of peach, almond and cherry on the Prunus linkage map. Abstract. Plant and Animal Genome Meeting, San Diego, CA. Bailey, CH, and LE Hough. 1975. Apricots, p. 375-380. In: J. Janick and 1N. Moore (eds). Advances in fruit breeding, West Lafayette, Indiana: Purdue University Press. Ballard, J.K., Probesting, EL, and Tukey, RB. 1971. Critical temperatures for blossom buds: Cherries. Wash. State Univ. Coop. Ext. Circl no. 371. Ballester, J ., Boskovic, R., Batlle, I., Arus, P., Vargas, F., de Vincente, M.C.1998. Location of the self-incompatibility gene on the almond linkage map. Plant Breeding. 117169-72. Baskerville GL and P Emin. 1969. Rapid estimation of heat unit accumulation from maximum and minimum temperatures. Ecology 50:514-517. Bliss, F.A., Arulsekar, S., Foolad, M.R., Becerra, V., Gillen, A.M., Warburton, M.L., Dandekar, A.M., Kocsisne, G.M., Mudin, K.K. 2002. An expanded linkage map of Prunus based on an interspecific cross between almond and peach. Genome. 45:520-529. Cantini, C., Iezzoni, A.F., Lamboy, W.F., Boritzki, M., Struss, D. 2001. DNA Fingerprinting of Tetraploid Cherry Germplasm Using Simple Sequence Repeats. Journal of American Society of Horticultural Science. 126(2):205-209. 49 Cipriani, G., Lot, 6., Huang, W.G., Marrazzo, M.T., Peterlunger, E., Testolin, R. 1999. AC/GT and AG/CT microsatellite repeats in peach [Prunus persica (L) Batsch]: isolation, characterization and cross-species amplification in Prunus. Theoratical Applied Genetics. 99: 65-72. Dennis, PG. and Howell, GS. 1974. Cold hardiness of tart cherry bark and flower buds. Res Rep MSU Ag Ext Station 220. Dirlewanger, E., Cosson, P., Tavaud, M., Aranzana, M.J., Poizat, C., Zanetto, A., Arus, P., Laigret, F., 2002. Development of microsatellite markers in peach [Prunus persica (L.) Batsch] and their use in genetic diversity analysis in peach and sweet cherry (Prunus avium L.) Theoretical and Applied Genetics. 105:127- 138. Downey, S. L., and Iezzoni, AF. 2000. Polymorphic DNA Markers in Black Cheny (Prunus serotina) Are Identified Using Sequences from Sweet Cherry, Peach, and Sour Cheny. Journal of American Society of Horticultural Science. 125(1):76-80. Etienne, C., Rothan, C., Moing, A., Plomion, C., Bodenes, C., Svanella-Dumas, L., Cosson, P., Pronier, V., Monet, R Dirlewanger, E. 2002. Candidate genes and QTLs for sugar and organic acid content in peach [Prunus persica (L.) Batsch]. Theoretical and Applied Genetics. 105: 145-159. Foulongne, M., Pascal, T., Arus, P., Kervella, J. 2003. The potential of Prunus davidiana for introgression into peach [Prunus persica (L.) Batsch] assessed by comparative mapping. Theoretical and Applied Genetics. 107:227-23 8. Hauck, N., Iezzoni, A., Yamane, H., Tao, R. 2001. Revisiting the S-allele Nomenclature in Sweet Cheny (Prunus avium) Using RFLP Profiles. Journal of American Society of Horticultural Sciences. 126(6):654-660. Hormaza, J .I., 2001. Identification of apricot (Prunus armeniaca L.) genotypes using microsatellite and RAPD markers. Acta Horticulturae. 546: 209-215. Hurtado, M.A., Romero, C., Vilanova, S., Abbott, A.G., Llacer, G., Badenes, ML. 2002. Genetic linkage maps of two apricot cultivars (Prunus armeniaca L.) and mapping of PPV sharka resistance. Theoreical and Applied Genetics. 105: 182-191. 50 Iezzoni, A., Schmidt, H., and Albertini, A. 1990. Cherries (Prunus spp.) pp. 110-173, In: JN Moore (ed.) Genetic Resources of Temperate Fruit and Nut Crops. ISHS, Wageningen, Netherlands. Iezzoni, AF. and Hamilton, KL. 1985. Differences in spring floral bud development among sour cherry cultivars. HortScience 20:915-916. Iezzoni, AF. and Mulinix, CA. 1992. Variation in bloom time in a sour cherry germplasm collection. HortScience. 27:1 1 13-1114. Jaglo, K.R., Klefi‘, S., Amundsem, K.L., Zhang, X., Haake, V., Zhang, J.Z., Deits, T., and Thomashow, M.F. 2001. Components of the Arabidopsis C-repeat/dehydration responsive element binding factor cold response pathway are conserved in Brassica napus and other plant species. Plant Physiology. 127:910-917. Joobeur, T., Periam, N., de Vicente, M.C., King, G.J., Arus, P. 2000. Development of a second generation linkage map for almond using RAPD and SSR markers. Genome. 43: 649-655. ' f . Jones, S. and Luchsinger, A. Plant systematics. McGraw-Hill, Inc: New York, New York. 1986. Kleweno, DD. and Matthews, V. 2003. Michigan Agricultural Statistics 2002-2003. Mich. Dept. Agric. 2002 Annual Report, 86 pp. Li, G., Quiros, CF. 2001. Sequence-related amplified polymorphism (SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica. Theoretical & Applied Genetics. 103:455-461. Lopes, M.S., Sefc, K.M., Laimer, M., Machado, A. 2002. Identification of microsatellite loci in apricot. Molecular Ecology Notes. 2:24-26. Lu, Z.X., Sosinski, B., Reighard, G.L., Baird, W.V., and Abbott, AG. 1998. Construction of a genetic linkage map and identification of AFLP markers for resistance to root-knot nematodes in peach rootstocks. Genome 41:199-207. 51 Michelmore, R.W., Paran, I., and Kessei, RV. 1991. Determining the linkage of quantitative trait loci to RFLP markers using extreme phenotypes of recombinant inbreds of soybean (Glycine max L. Merr). Theoretical and Applied Genetics. 86: 914-918. _ Michigan Balaton Cherries official Website (last updated 4/9/2004) http://www.michig2_1nbalatoncherries.com/ Michigan State University Oflicial Balaton Website (last updated 5/13/2004) http://wwwhrtmsu.edu/Balaton.html Miklas, P.N., Johnson, E., Stone, V., Beaver, J .S., Montoya, C., and Zapata, M. 1996. Selective mapping of QTL conditioning disease resistance in common bean. Crop Science 36: 1344-1351. Olden, E.J., and Nybom, N. 1968. On the origin of Prunus cerasus L. Hereditas 59: 327- 345. Owens, CL. 2001. The role of CBF orthologs in the cold acclimationo f sour cherry (Prunus cerasus) and strawberry (F ragaria x ananassa). PhD Thesis. Mich. State Univ., 7 5 pp. Owens, C.L., Thomashow, M.F., Hancock, IF, and Iezzoni, AF. 2002. CBF] orthologs in sour cherry and strawberry and the heterologous expression of CBF] in strawberry. Journal of American Society of Horticultural Science. 127(4):489-494. Rajapakse, S., Byme, D.H., Zhang, L., Anderson, N., Arumuganathan, K., Ballard, RE. 2001. Two genetic linkage maps of tetraploid roses. Theroretical and Applied Genetics. 103:575-583. Raman, H., Moroni, J.S., Sato, K., Read, B.J., Scott, B]. 2002. Identification of AFLP and microsatellite markers linked with an aluminum tolerance gene in barley (Hordeum vulgare L.). Theoretical and Applied Genetics. 1051458-464. Ricks, D. 1992. Supply fluctuations and long term acreage cycles in tart cherries and other tree crop industries. Great Lakes Fruit Growers News 31:22,24. 52 SAS Institute, Inc. 1999. Version 8.00. Cary, NC. Schueler, S., Tusch, A., Schuster, M., and Ziegenhagen, B. 2003. Characterization of microsatellites in wild and sweet cherry (Prunus avium L.) — markers for individual identification and reproductive processes. Genome. 46:95-102. Sosinski, B., Gannavarapu, M., Hager, L.D., Beck, L.E., King, G.J., Ryder, C.D., Rajapakse, S., Baird, W.V., Ballard, R.E., Abbott, AG. 2000. Characterization of microsatellite markers in peach [Prunus persica (L.) Batsch]. Theoretical and Applied Genetics. 101:421-428. Staub, 1E. and Serquen, F. 1996. Genetic Markers, Map Construction, and their applications in Plant Breeding. HortScience. 31(5):729-74l. Stockinger, E.J., Mulinix, CA, Long, C.M., Brettin, TS. and Iezzoni. AF. 1996. A linkage map of sweet cherry based on RAPD analysis of a microspore-derived callus culture population. Journal of Heredity. 87:214-218. Stniss, D., Boritzki, M., Karle, R., and Iezzoni, AF. 2002. Microsatellite markers differentiate eight Giessen cherry rootstocks. HortScience. 37: 191-193. Testolin, R., Marrazzo, T., Cipriani, G., Quarta, R.,Verde, I., Dettori, M.T., Pancaldi, M., Sansavini, S. 2000. Microsatellite DNA in peach (Prunus persica L. Batsch) and its use in fingerprinting and testing the genetic origin of cultivars. Genome. 43 :5 12-520. Vilanova, S., Romero, C., Abbott, A.G., Llacer, G., Badenes, ML. 2003. An apricot (Prunus armeniaca L.) F2 progeny linkage map based on SSR and AFLP markers, mapping plum pox virus resistance and self-incompatibility traits. Theoretical and Applied Genetics. 107:239-247. Wang, D., Karle, R., Brettin, TS, and Iezzoni, AF. 1998. Genetic linkage map in sour cherry using RFLP markers. Theoretical and Applied Genetics. 97: 1217-1224. Wang, D., Karle, R., Iezzoni, AF, 2000. QTL analysis of flower and fi'uit traits in sour cherry. Theoretical and Applied Genetics. 1002535-544. 53 Watkins R 1976. Chery, plum, peach, apricolt and almond, p. 242-247. In: NW Simmonds (eds). Evolution of crop plants. Longrnan, NY. Website for projects fiinded under Area 4 of the FAIR programme (Europa) http ://europa. eu . int/comm/research/ago/fg’r/en/fi4 1 3 9. htm] Wunsch, A., and Hormaza, J .I. 2002. Molecular cheracterisation of sweet cherry (Prunus avium L.) genotypes using peach [Prunus persica (L.) Batsch] SSR sequences. Heredity. 89:56-63. . Yamamoto, T., Mochida, K., Imai, T., Shi, Y.Z., Ogiwara, I., and Hayashi, T., 2002. Microsatellite markers in peach [Prunus persica (L.) Batsch] derived from an enriched genomic and cDNA libraries. Molecular Ecology Notes. 2:298-301. Yamamoto, T., Shimada, T., Imai, T., Yaegaki, H., Haji, T., Matsuta, N., Yamaguchi, M., Hayashi, T., 2001. Characterization of Morphologival Traits Based on Genetic Linkage Map in Peach. Breeding Science. 51:271-278. Yamane, H., Tao, R., Sugiura, A., Hauck, N., Iezzoni, AF. 2001. Identification and Characterization of S-RNases in Tetraploid Sour Cherry (Prunus cerasus). Journal of American Society of Horticultural Sciences. 126(6):661-667. 54 A“ “j“,