Iowa State University Capstones, Theses and Graduate Theses and Dissertations Dissertations

2010 Soybean improvement strategies: Insect-pollinator attraction and genetic resistance to whitefly and to brown stem rot Paola T. Perez Iowa State University

Follow this and additional works at: https://lib.dr.iastate.edu/etd Part of the Agronomy and Crop Sciences Commons

Recommended Citation Perez, Paola T., "Soybean improvement strategies: Insect-pollinator attraction and genetic resistance to whitefly nda to brown stem rot" (2010). Graduate Theses and Dissertations. 11402. https://lib.dr.iastate.edu/etd/11402

This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Graduate Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Soybean improvement strategies: Insect-pollinator attraction and genetic resistance to whitefly and to brown stem rot

by

Paola Tatiana Perez

A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY

Major: Breeding

Program of Study Committee: Reid G. Palmer, Co-Major Professor Silvia R. Cianzio, Co-Major Professor Dianne Cook John Nason Randy Shoemaker

Iowa State University

Ames, Iowa

2010 ii

TABLE OF CONTENTS

LIST OF TABLES ...... iv LIST OF FIGURES ...... vi ACKNOWLEDGEMENTS ...... viii CHAPTER 1: GENERAL INTRODUCTION...... 1 INTRODUCTION ...... 1 LITERATURE REVIEW ...... 3 Brown stem rot (BSR) ...... 3 Whitefly ...... 4 Soybean hybrids...... 5 Wild perennial species...... 7 DISSERTATION ORGANIZATION ...... 8 REFERENCES ...... 9 CHAPTER 2: GENETIC ANALYSIS OF NEW SOURCES OF SOYBEAN RESISTANCE TO BROWN STEM ROT...... 16 ABSTRACT...... 16 INTRODUCTION ...... 17 MATERIALS AND METHODS...... 20 Plant material ...... 20 Inoculation protocol...... 21 Severity assays...... 21 DNA marker analysis...... 23 Data analysis ...... 23 RESULTS ...... 24 Phenotypic analysis...... 24 Single marker analysis ...... 26 DISCUSSION...... 27 Acknowledgements...... 29 REFERENCES ...... 30 CHAPTER 3: QTL MAPPING OF WHITEFLY RESISTANCE IN SOYBEAN ...... 43 ABSTRACT...... 43 INTRODUCTION ...... 44 MATERIALS AND METHODS...... 47 Population development...... 47 Phenotypic evaluation...... 48 DNA marker analysis...... 49 Statistical analysis...... 49 RESULTS ...... 51 Phenotypic analysis...... 51 Genotypic analysis...... 52 QTL analysis...... 53 iii

DISCUSSION...... 54 Acknowledgements...... 57 REFERENCES ...... 58 CHAPTER 4: NECTAR COMPOSITION AND FLOWER CHARACTERISTICS OF WILD PERENNIAL GLYCINE SPECIES ...... 71 ABSTRACT...... 72 INTRODUCTION ...... 73 MATERIALS AND METHODS...... 76 Study ...... 76 Study site and methods ...... 77 Nectar analysis and flower evaluation ...... 77 Data analysis ...... 78 RESULTS ...... 79 Summary statistics ...... 79 Nectar and flower characteristics of perennial Glycine species...... 80 Analysis of variance...... 81 DISCUSSION...... 82 Acknowledgements...... 85 REFERENCES ...... 85 CHAPTER 5: GENERAL CONCLUSIONS ...... 106

iv

LIST OF TABLES

CHAPTER 2 Table 1. Visual assessment scale of brown stem rot (BSR) severity in soybean plants five weeks after inoculation...... 36 Table 2. Brown stem rot (BSR) resistance ratings for plate susceptibility and vigor in four F2 populations and their parental lines...... 37

Table 3. Number of F2:3 lines with higher brown stem rot (BSR) resistance than the most resistant parent or lower than the most susceptible parent by two standard deviations...... 38 Table 4. Means of genotypic classes and R2 values of molecular markers on chromosome 16 (formerly molecular linkage group J) associated with brown stem rot (BSR) resistance...... 39

CHAPTER 3 Table 1. Population mean, standard error, and tests of significance of the difference between genotypic classes, for whitefly infestation in populations Williams 79 x Corsoy 79, and Williams 79 x Cajeme, in 2003, 2004, and the combined analysis of both years...... 64 Table 2. Summary of significant and suggestive whitefly resistance quantitative trait loci (QTL) for populations Williams 79 x Corsoy 79, and Williams 79 x Cajeme, in 2003, 2004, and in the combined analysis across years...... 65 CHAPTER 4 Table 1. Species evaluated of the genus Glycine Willd., subgenus Glycine, and the interspecific hybrid G. max x G. tomentella, ploidy level, genome type, and geographic distribution...... 91 Table 2. Nectar sugar content and correlations across Glycine species and the interspecific hybrid G. max x G. tomentella...... 92 Table 3. Sugar ratio, sugar percentages, and correlations across Glycine species and the interspecific hybrid G. max x G. tomentella...... 93 Table 4. Flower size across species and correlations across Glycine species and the interspecific hybrid G. max x G. tomentella...... 94 Table 5. Nectar sugar content [mean (s.d.)] of Glycine species and the interspecific hybrid G. max x G. tomentella...... 95 v

Table 6. Sugar ratio and sugar percentages [mean (s.d.)] of Glycine species and the interspecific hybrid G. max x G. tomentella...... 96 Table 7. Flower size and petal perimeter [mean (s.d.)] of Glycine species and the interspecific hybrid G. max x G. tomentella...... 97 Table 8. Multivariate analysis of variance (MANOVA) for nectar sugar content, using species with a sample size larger than ten. Wilk's lambda test and p-value (in parentheses) are reported. Wilk's lambda test and p-value of a pair of species that are non-significantly different in nectar sugar content are reported in bold...... 98 Table 9. Multivariate analysis of variance (MANOVA) for nectar sugar ratio and percentages, using species with a sample size larger than ten. Wilk's lambda test and p-value (in parentheses) are reported. Wilk's lambda test and p-value of a pair of species with non-significant differences in nectar sugar proportions are reported in bold...... 99 Table 10. Pair wise comparisons of flower size, using species with sample size three or larger. Student t-test and p-values (in parentheses) are reported. t-test and p-value of a pair of species with non-significant differences in flower size are reported in bold...... 100

vi

LIST OF FIGURES

CHAPTER 2 Figure 1. Plate susceptibility. A) Susceptible plant with infection of P. gregata in stem fragments from the inoculation point and from the plant top. B) Resistant plant with infection only at the inoculation point...... 40 Figure 2. Distribution of plate susceptibility (1 = infection at inoculation point and at plant top; 0 = infection only at inoculation point) for four segregating populations and their parental lines. Arrows indicate parental mean values (continuous arrows for Century 84 and dashed arrows for PI 594637, PI 594638B, PI 594650A, and PI 594658B)...... 41 Figure 3. Distribution of vigor (measured on a scale from 1 to 7; 1= dead plant; 7 = healthy plant) for four segregating populations and their parental lines. Arrows indicate parental mean values (continuous arrows for Century 84 and dashed arrows for PI 594637, PI 594638B, PI 594650A, and PI 594658B)...... 42 CHAPTER 3

Figure 1. Distribution of whitefly infestation among 150 F2:3 lines from population Williams 79 x Corsoy 79. Histograms show distribution of mean nymph density (number of nymphs in an area of 6.5 cm2) in 2003, 2004, and in the combined analysis of both years. Arrows indicate parental mean values (continuous arrows for Williams 79 and dashed arrows for Corsoy 79)...... 66

Figure 2. Distribution of whitefly infestation among 90 F2:3 lines from the population Williams 79 x Cajeme. Histograms show distribution of mean nymph density (number of nymphs in an area of 6.5 cm2) in 2003, 2004, and in the combined analysis of both years. Arrows indicate parental mean values (continuous arrows for Williams 79 and dashed arrows for Cajeme)...... 67 Figure 3. Logarithm of Odds (LOD) curves from interval mapping analysis in the combined analysis of two years for Williams 79 x Corsoy 79, and in 2003 for Williams 79 x Cajeme. The dashed horizontal lines indicate the genome-wide 95% LOD threshold...... 68 Figure 4. Nymph density (number of nymphs in an area of 6.5 cm2) for different genotypes at the two SSR markers flanking a QTL on chromosome 18 (formerly molecular linkage group G) in the combined analysis. Error bars represent 95% confidence intervals. A: alleles from Williams 79, and B: alleles from Corsoy 79...... 69 Figure 5. Nymph density (number of nymphs in an area of 6.5 cm2) for different genotypes at the two SSR markers flanking a QTL on chromosome 12 (formerly vii

molecular linkage group H) in 2003. Error bars represent 95% confidence intervals. A: alleles from Williams 79 and B: alleles from Cajeme...... 70 CHAPTER 4 Figure 1. Typical soybean flower. A) Standard petal. B) Two wing petals. C) Fused keel petals...... 101 Figure 2. Distribution of nectar sugar content of perennial Glycine species with a sample size larger than ten. Mean value represented by X symbol...... 102 Figure 3. Distribution of nectar sugar ratio and sugar percentages of perennial Glycine species with a sample size larger than ten. Mean value represented by X symbol...... 103 Figure 4. Distribution of petal perimeter and flower size of perennial Glycine species with a sample size larger than ten. Mean value represented by X symbol...... 104 Figure 5. Geographic distribution of G. clandestina, G. cyrtoloba, G. latifolia, G. microphylla, G. tomentella, G. tabacina, G. canescence, G. falcata, and G. lactovirens...... 105

viii

ACKNOWLEDGEMENTS

First of all, I want to thank my major professors, Dr. Reid Palmer and Dr. Silvia

Cianzio. They truly have been mentors to me during this stage of my life, both academically and personally. They have kept me motivated with their enthusiasm and willingness to share all their knowledge and experience. Thanks to them I had the opportunity to continue with my Ph.D. program at Iowa State Univeristy and to attend several national and international meetings, which improved my professional life (and my personal life too, thanks Dr. Palmer for your suggestion to prepare business cards!). Many thanks also to the professors who served as POS committee members: Dr. Dianne Cook, Dr. John Nason, and Dr. Randy

Shoemaker. Their comments have helped me to reshape my research projects.

I would like to thank my family for always believing in me (special thanks to my grandma who strongly believes that all her prayers made my journey possible!). Thanks to my friends in Colombia and in the U.S. for always being there, no matter what. Last, but not least, thanks to my husband Adam, for his love, patience, and encouragement.

1

CHAPTER 1: GENERAL INTRODUCTION

INTRODUCTION

Soybean [Glycine max (L.) Merr.] is one of the most important crop plants cultivated worldwide. It is a source of oil and protein and is used for livestock feed, human food, and also for industry. The major producers are the United States, which accounts for 34% of world production, Brazil with 27%, Argentina with 22%, and China with 6%, (FAO

Statistics, 2010). Research goals are related to increases in grain yield along with desirable quality characteristics such as high seed protein and seed oil content. The majority of farmers in the U.S. derive their income on the basis of seed yield, which is the final plant expression of all metabolic processes and interactions with environment (biotic and abiotic factors). Because of this, yield is the main trait in breeding improved cultivars, both in the public and private sector.

Over the years, yield increases have been accomplished due to genetic improvement and enhancement of agronomic practices. These increases have been cited extensively in the literature. Boerma (1979) reported average yield increases of 0.7% per year for cultivars that were released from 1914 to 1973 in the southern U.S. In a similar study, Luedders (1977) reported a 26% yield increase of cultivars released from 1940 to 1950 compared to plant introductions, and a 16% yield increase from cultivars released from 1960 to 1970 compared with those released from 1940 to 1950. Specht and Williams (1984) reported a soybean yield increase of 21 kg ha-1 annually in the period from 1924 to 1980, with the majority of the 2

yield increase attributed to breeding. Wilcox (2001) estimated that public soybean breeders in the northern U.S. soybean production region have increased seed yield approximately 60% over the past 60 years. Specht et al. (1999) estimated that similar progress in breeding also has been made in the private sector, although these results were obtained during the last two to three decades. In general, yield improvement may be achieved by breeding directly for the trait, i.e. crossing high-yielding parents to bring together yield genes and selecting at every stage of testing the highest producing individuals. Additionally, improved yield can be obtained indirectly by improving other aspects of plant production, i.e. physiological aspects and plant health. Soybean breeding efforts for yield have been made by applying different strategies related to selection methods (Piper and Fehr, 1987; Sumarno and Fehr, 1982;

Uphoff et al., 1997; Voigt and Weber, 1960) and to population development (Cerna et al.,

1997; Thompson and Nelson, 1998).

Soybean breeding programs have been focused on the development of inbred lines, which represent the totality of cultivars available in the market. However, increases in soybean production may be possible through the development of hybrids (Brim and

Cockerham, 1961; Burton and Brownie, 2006; Cerna et al., 1997; Lewers, 1996; Manjarrez-

Sandoval et al., 1997; Nelson and Bernard, 1984; Ortiz-Perez, et al., 2007; Palmer et al.,

2001, 2003; Pandini et al. 2002; Perez, et al., 2009a; Sun et al., 1999). Breeding for resistance to diseases and pests is an extremely important objective to improve soybean yield.

Soybean is threatened by insect pests (i.e. whitefly) that impact soybean yield in commercial plantings, or in research plantings under confined conditions, such as greenhouses. In 3

addition, several diseases [i.e. brown stem rot (BSR)] lowers yield, which decreases production and results in less income (Wrather et al., 2001).

The objectives of this study were: 1) to determine if four plant introductions (PI) from south-central China, identified as resistant to the disease brown stem rot (BSR), represent new sources of BSR resistance genes, 2) to identify quantitative trait loci (QTL) for whitefly resistance in soybean, and 3) to determine nectar composition (sucrose, glucose, fructose, and total carbohydrates) and flower size of greenhouse-grown wild perennial plants.

LITERATURE REVIEW

Brown stem rot (BSR)

BSR is an economically important disease prevalent in soybean producing regions in the north-central U.S. and Canada (Gray and Grau, 1999). Under conditions that favor disease development, yield losses of up to 38% have been reported (Gray, 1972).

Phialophora gregata f. sp. sojae is the causal organism, infecting plants through the roots, colonizing the pith and vascular system and moving up the stem into the leaves of susceptible genotypes (Allington and Chamberlain, 1948; Schneider et al., 1972). Two types of

Phialophora gregata f. sp. sojae have been identified. Type I, causes leaf necrosis and defoliation, and stem vascular discoloration, and Type II, causes vascular discoloration only

(Gray, 1972; Harrington et al., 2003).

The most effective way to manage the disease is by planting BSR-resistant soybean cultivars (Bachman et al., 1997; Mengistu et al., 1986). To date, three BSR resistance genes 4

(Rbs1, Rbs2, and Rbs3) have been found and these genes map to the same general region on chromosome 16 (formerly molecular linkage group J) (Bachman et al. 2001; Klos et al.,

2000; Lewers et al., 1999; Patzoldt et al, 2005; Webb, 1997). The identification of new sources of resistance to BSR is important for soybean breeding programs. Plant introductions

(PI) have been valuable sources of genetic diversity in soybean improvement because they can carry genes different than those found in commercial cultivars (Carter et al., 2004).

Whitefly

Whitefly, Bemisia tabaci (Gennadius) (order: Homoptera, family: Aleyrodidae), is a common economic pest on many crops throughout the world. Historically, whiteflies have been major pests in vegetable crops, cotton, and greenhouse plants. In addition, they have a wide range of hosts, including wild plant species, which makes it difficult to control once they have infested a crop (Perring, 2001). The whitefly species complex is formed by highly cryptic sibling species (Campbell et al., 1996; Perring, 2001). To date, 41 distinct populations of B. tabaci have been studied; 24 of these populations have been given a specific biotype designation (Perring, 2001). Adults and nymphs have piercing-sucking mouthparts and feed on the lower leaf surfaces where they lay their eggs (Byrne and Bellows, 1991; Perring,

2001). Whiteflies damage crops by extracting large quantities of phloem sap (Jones, 2003), and heavy infestations may result in development of chlorotic spots on leaves, wilting, and stunting of plants. In addition, these insects excrete a sticky material called honeydew, which in high concentrations promotes the growth of sooty mold fungi, which interferes with photosynthesis. Some species of the whitefly complex serve as vectors of several economically important viral plant pathogens (Byrne and Bellows, 1991). 5

Host plant resistance to whitefly can be due to physical characteristics of the leaf surface, e.g. hairiness vs. glabrousness, sticky glandular trichomes, leaf shape, microclimate due to foliage density; and chemical characteristics, such as pH of leaf sap (Berlinger, 1986).

However, the main reason for the resistance to whitefly is unknown.

In soybean, whitefly has become an important pest problem in tropical locations and in greenhouse conditions (Costa, 1976; Vaishampayan et al. 1975). The sweetpotato whitefly, Bemisia tabaci (Gennadius), and the bandedwinged whitefly, Trialeurodes abutiloneus (Haldeman) have been recorded colonizing field-grown soybean (Vaishampayan et al. 1975). In greenhouse conditions, soybean is colonized by Trialeurodes vaporariorum

(Westwood). Viruses such as soybean crinkle mosaic, and soybean dwarf mosaic can be transmitted by whiteflies (Costa, 1976). Breeding of soybean cultivars resistant to this pest is the most effective management alternative (Arioglu et al., 1989).

Soybean hybrids

The interest in hybrid soybean developed after the identification of the first male- sterile, female-fertile mutant (Brim and Young, 1971). Its use in recurrent selection breeding programs (Brim and Stuber, 1973), increased the awareness of the potential to produce commercial hybrid soybean. According to Palmer et al. (2001) there are five components that are crucial for the successful development of commercial hybrid soybean:

1. Parental combinations that produce heterosis levels superior to the best pure-line cultivars.

2. A stable male-sterile, female-fertile sterility system.

3. A selection system to obtain 100% female (pod parent) plants that set seed normally and can be harvested mechanically. 6

4. An efficient pollen transfer mechanism from pollen parent to pod parent.

5. An economical level of seed increase for the seedsman and growers that ultimately benefits the consumer.

Regarding the first requirement for hybrid soybean production, several heterosis studies have showed that heterosis levels, above the high parent, are possible (Brim and

Cockerham, 1961; Burton and Brownie, 2006; Cerna et al., 1997; Lewers, 1996; Manjarrez-

Sandoval et al., 1997; Nelson and Bernard, 1984; Ortiz-Perez et al., 2007; Palmer et al.,

2001; Pandini et al. 2002; Perez et al., 2009a, 2009b; Sun et al., 1999). In order to meet the second requirement for soybean hybrid production (a stable male-sterile, female-fertile sterility system), mutations affecting male cell and organ development have generated male- sterile, female-fertile lines that can be used as female parents for hybrid seed production

(Palmer, 2000). The third requirement is to have a selection system to obtain 100% female plants because in hybrid seed production fields, female rows will be segregating for the male sterility (ms) mutation. According to Palmer et al. (2003), any of the selection systems employed with the soybean nuclear male-sterile genes, seed size differential (Carter et al.,

1984), linkage between genes controlling the green cotyledon trait and the Ms5 locus (Burton and Carter, 1983), as well as the W1 flower color locus and the Ms6 locus (Lewers and

Palmer, 1997; Lewers, 1996; Lewers et al., 1998a, 1998b) are suitable to identify male- sterile, female-fertile plants for hybrid seed production.

Upon obtaining a stable male sterility system, it is necessary to transfer the pollen from the male parent to the female parent. An efficient pollen transfer mechanism is the most limiting factor for soybean hybrid production (Brim, 1975; Davis, 1975; Palmer et al., 2001,

2003). In soybean, manual cross-pollination to produce large quantities of hybrid seed is 7

difficult and time consuming. The small size of the soybean flowers, the low success rate and the few seeds obtained per hybrid pod contribute to the difficulty of manually producing large quantities of hybrid seed (Fehr, 1991). Insect cross-pollination of male-sterile soybean plants facilitates the production of hybrid seed (Lewers, 1996; Nelson and Bernard, 1984;

Ortiz-Perez et al., 2007). Pollinator insects such as honeybees Aphis melliphera and alfalfa leaf cutter bee, Megachile rotundata F. are used in hybrid soybean production. Also, some wild native bees primarily from families Megachilidae, Halictidae, Anthophoridae, and

Andrenidae could be efficient pollinators of soybean flowers (Ortiz-Perez et al., 2007).

Soybean is an obligate self-pollinated crop with less than 1% natural cross-pollination

(Palmer et al., 2001). Before soybean hybrids can become a reality, it is necessary to understand some of the mechanisms for insect attraction and reward. Wild relatives of the cultivated soybean represent a germplasm reservoir for agronomic improvement of the cultivated species (Singh and Hymovitz, 1999).

Wild perennial Glycine species

The genus Glycine (Willd.) is composed of two subgenera: Soja (Moench) F. J.

Herm. and Glycine. The subgenus Soja contains the cultivated soybean Glycine max (L.)

Merr. and the wild annual Glycine soja (Sieb. and Zucc.). The subgenus Glycine has approximately 25 defined wild perennial species distributed throughout Australia and islands of the Pacific Ocean, extending to Taiwan. The unexploited valuable pool of genetic diversity of the wild perennial species is largely unexplored by soybean breeders (Singh and

Hymowitz, 1999). Traits of interest that could be incorporated into the cultivated soybean 8

include resistance to soybean rust, soybean brown spot, powdery mildew, phytophthora root rot, white mold, sudden death syndrome, tobacco ring spot virus, yellow mosaic virus, soybean cyst nematode, and tolerance to certain herbicides and to salt (Palmer and

Hymowitz, 2004).

Glycine species are predominantly self-pollinated, although they possess chasmogamous flowers capable of outcrossing. Within the wild perennial species Glycine argyrea (Tind.), Brown et al. (1986) recorded outcrossing rates among chasmogamous flowers from zero to complete outcrossing, with an average of 40%. The identification and characterization of sources of insect attraction among the wild perennial species of the genus

Glycine could be useful in soybean breeding programs for hybrid production, because in contrast with the low natural out-crossing in the cultivated soybean, out-crossing in the wild perennial soybean species can exceed 50%.

DISSERTATION ORGANIZATION

This dissertation is organized into five chapters. Chapter one is the general introduction. Chapter two is a manuscript submitted to Crop Science entitled "Genetic analysis of new sources of soybean resistance to brown stem rot". Chapter three is a manuscript to be submitted to Journal of Crop Improvement entitled "QTL mapping of whitefly resistance in soybean". Chapter four is a manuscript to be submitted to Journal of

American Botany entitled "Nectar composition and flower characteristics of wild perennial

Glycine species". Chapter five is the general conclusion of this study. Different objectives will be addressed in chapters two, three, and four, each of them with the ultimate purpose to 9

contribute to improve yield in soybean. Chapters two and three each have objectives that relate to plant health, as a means to improve yield potential through genetic resistance to a fungal disease and an insect pest. Chapter two describes the search for new sources of resistance to brown stem rot in four plant introductions from south-central China. Chapter three describes quantitative trait loci mapping of whitefly resistance in two segregating populations. Chapter four evaluates factors that affect insect-pollinator attraction and reward in the wild perennial Glycine species, which could be used to improve pollen transfer in the cultivated soybean for hybrid production.

REFERENCES

Arioglu, H.H., A.F. Ozgur, and N. Isler. 1989. Influence of soybean pubescence type and

density on whitefly (Bemisia tabaci Genn.) resistance. In: A.J. Pascale (ed.), Proc.

World Soybean Res. Conf. IV, Argentina, 5-9 March 1989, pp. 1235.

Allington, W.B., and D.W. Chamberlain. 1948. Brown stem rot of soybean. Phytopathology

38: 793–802.

Bachman, M.S., C.D. Nickell, P.A. Stephens, and A.D. Nickell. 1997. Brown stem rot

resistance in soybean germplasm from central China. Plant Dis. 81:953–956.

Bachman, M.S., J.P. Tamulonis, C.D. Nickell, and A.F. Bent. 2001. Molecular markers

linked to brown stem rot resistance genes Rbs1 and Rbs2 in soybean. Crop Sci.

41:527–535.

Berlinger, M.J. 1986. Host plant resistance to Bemisia tabaci. Agriculture, Ecosystems and

Environment 17: 69-82.

Boerma, H.R. 1979. Comparison of past and recently developed soybean cultivars in 10

maturity groups VI, VII, and VIII. Crop Sci. 19:611-613.

Brim, C.A. and M.F. Young. 1971. Inheritance of a male-sterile character in soybeans. Crop

Sci. 11: 564-566.

Brim, C.A. and C.W. Stuber. 1973. Application of genetic male sterility to recurrent selection

schemes in soybeans. Crop Sci. 13: 528-530.

Brim, C.A. 1975. Implications of male sterility in soybean. Proc. Sixth Soybean Seed Res.

Conf., Am. Seed Trade Assoc. 6: 67-71.

Brim, C.A., and C.C. Cockerham. 1961. Inheritance of quantitative characters in soybeans.

Crop Sci. 1: 187-190.

Brown, A. H. D., J. E. Grant, and R. Pullen. 1986. Outcrossing and paternity in Glycine

argyrea by paired fruit analysis. Biol. J. Linn. Soc. 29: 283-294.

Burton, J.W., and C. Brownie. 2006. Heterosis and inbreeding depression in two soybean

single crosses. Crop Sci. 46: 2643-2648.

Burton, J.W., and T.E., Jr. Carter. 1983. A method for production of experimental quantities

of hybrid soybean seed. Crop Sci. 23: 388-390.

Byrne, D.N. and T.S. Bellows. 1991. Whitefly biology. Ann. Rev. Entomol. 36: 431-457.

Campbell, B.C., J.D. Seffen-Campbell, and R.J. Gill. 1996. Origin and radiation of

whiteflies: an initial molecular phylogenetic assessment. In: D. Gerlin, and R.T. Mayer

(eds.), Bemisia: 1995 , biology, damage, control, and management. Intercept

Ltd., Andover, Hants, UK, pp. 29-51.

Carter, T.E., Jr., J.W. Burton, and E.B. Jr., Huie. 1984. Mechanical separation of seed from

male-sterile and fertile plants by seed size. Soybean Genet. Newsl. 11: 146-149. 11

Carter Jr, T.E., R.L. Nelson, C. Sneller, and C. Zhanglin. 2004. Genetic diversity in

soybean. H.R. Boerma and J.E. Specht (eds), Soybean Monograph, 3rd Edition, Am.

Soc. of Agron., Madison, WI. pp. 114.

Cerna F.J., S.R. Cianzio, A. Rafalski, S. Tingey, and D. Dyer. 1997. Relationship between

seed yield heterosis and molecular heterozygosity in soybean. Theor. Appl. Genet. 95:

460-467.

Costa, A.S. 1976. Whitefly-transmitted plant diseases, Ann. Rev. Phytopathol. 16: 429-

449.

Davis, W.H. 1975. Technology of hybrid soybeans. Proc. Sixth Soybean Seed Res. Conf.,

Am. Seed Trade Assoc. 6: 72-74.

FAOSTAT data. 2010. http://faostat.fao.org/faostat.

Fehr, W. R. 1991. Principles of cultivar development. Theory and technique. Macmillan

Publishing Company, Ames, IA.

Gray, L.E. 1972. Effect of Cephalosporium gregatum on soybean yield. Plant Dis. Rep. 56:

580-581.

Gray, L.E., and C.R. Grau. 1999. Brown stem rot. In: G.L. Hartman, J.B. Sinclair, and J.C.

Rupe (eds.), Compendium of soybean diseases. 4th ed. APS Press. St. Paul, MN.

Harrington, T.C., J. Steimel, F. Workneh, and X.B. Yang. 2003. Characterization and

distribution of two races of Phialophora gregata in the north-central United States.

Phytopathology 93: 901-912.

Jones, D. 2003. Plant viruses transmitted by whiteflies. Euro. J. Plant Pathol. 109: 195-219. 12

Klos, K.L.E., M.M. Paz, L.F. Marek, P.B. Cregan, and R.C. Shoemaker. 2000. Molecular

markers useful for detecting resistance to brown stem rot in soybean. Crop Sci. 40:

1445–1452.

Lewers, K. S. 1996. Production, evaluation, and utilization of hybrid soybean [Glycine max

(L.) Merr.], Ph.D. Diss., Iowa State University, Ames (Diss. Abstr. 96-26046).

Lewers, K.S., and R.G. Palmer. 1997. Recurrent selection in soybean. Plant Breed. Rev. 16:

275-313.

Lewers, K.S., S.K. St. Martin, B.R. Hedges, and R.G. Palmer. 1998a. Effects of the Dt2 and S

alleles on agronomic traits of F1 hybrid soybean. Crop Sci. 38: 1137-1142.

Lewers, K.S., S.K. St. Martin, B.R. Hedges, and R.G. Palmer. 1998b. Testcross evaluation of

soybean germplasm. Crop Sci. 38: 1143-1149.

Lewers, K.S., E.H. Crane, C.R. Bronson, J.M. Schupp, P. Keim, and R.C. Shoemaker. 1999.

Detection of linked QTL for soybean brown stem rot resistance in ‘BSR 101’ as

expressed in a growth chamber environment. Molec. Breed. 5: 33–42.

Luedders, V.D. 1977. Genetic improvement in yield of soybeans. Crop Sci. 17: 971-972.

Manjarrez-Sandoval, P., T.E. Carter, Jr., D.M. Webb, and J.W. Burton. 1997. Heterosis in

soybean and its prediction by genetic similarity measures. Crop Sci. 37: 1443-1452.

Mengistu, A., C.R. Grau, and E.T. Gritton. 1986. Comparison of soybean genotypes for

resistance to and agronomic performance in the presence of Brown Stem Rot. Plant Dis.

70: 1095-1098.

Nelson, R.L. and R.L. Bernard. 1984. Production and performance of hybrid soybeans. Crop

Science 24: 549-553.

Ortiz-Perez, E., S.R. Cianzio, H. Wiley, H.T. Horner, W.H. Davis, and R.G. Palmer. 2007. 13

Insect-mediated cross-pollination in soybean [Glycine max (L.) Merr.]: I. Agronomic

performance. Field Crops Res. 101: 259-268.

Palmer, R.G. 2000. Genetics of four male-sterile, female-fertile soybean mutants. Crop Sci.

40: 78-83.

Palmer, R.G., E. Ortiz-Perez, I.G. Cervantes-Martinez, H. Wiley, S.J. Hanlin, R.A. Healy, T.

Horner, and W.H. Davis. 2003. Hybrid soybean-current status and future outlook. 33rd

Soybean Seed Research Conference. American Seed Trade Association. Seed Expo 2003.

Available in CD-ROM.

Palmer, R.G., J. Gai, H. Sun, and J.W. Burton. 2001. Production and evaluation of hybrid

soybean. Plant Breed. Rev. 21: 263-307.

Palmer, R. G. and T. Hymowitz. 2004. Soybean Genetics and Breeding. In: C. Wrigley, H.

Corke, and C. Walker (eds.), Encylopedia of grain sciences. Academic Press, London,

England, pp. 136-146.

Pandini, F., A.V. Natal, A. Celis de Almeida Lopes. 2002. Heterosis in soybeans for seed

yield components and associated traits. Braz. Arch. Biol. Tech. 45: 401-412.

Patzoldt, M.E., C.R. Grau, P.A. Stephens, N.C. Kurtzweil, S.R. Carlson, and B.W. Diers.

2005. Localization of a quantitative trait locus providing brown stem rot resistance in the

soybean cultivar Bell. Crop Sci. 45: 1241-1248.

Perez. P.T., S. R. Cianzio, and R.G. Palmer. 2009a. Evaluation of soybean [Glycine max (l.)

Merr.] F1 hybrids. J. Crop Improve. 23: 1-18.

Perez. P.T., S.R. Cianzio, E. Ortiz-Perez, and R.G. Palmer. 2009b. Agronomic performance

of soybean hybrids from single, three-way, four-way, and five-way crosses, and

backcross populations. J. Crop Improve. 23: 95-118. 14

Perring, T.M. 2001. The Bemisia tabaci species complex. Crop Prot. 20: 725-737.

Piper, T.E., and W.R. Fehr. 1987. Yield improvement in a soybean population by utilizing

alternative strategies of recurrent selection. Crop Sci. 27: 172-178.

Schneider, R.W., J.B. Sinclair, and L.E. Gray. 1972. Etiology of Cephalosporium gragatum

in soybean. Phytopathology 62: 345-349.

Singh, R. J., and T. Hymowitz. 1999. Soybean genetic resources and crop improvement.

Genome 42: 605-616.

Specht, J.E., D.J. Hume, and S.V. Kumundini. 1999. Soybean yield potential-A genetic and

physiological perspective. Crop Sci. 39:1560-1570.

Specht, J.E., and J.H. Williams. 1984. Contribution of genetic technology to soybean

productivity-retrospect and prospect. In: W.R. Fehr (ed.). Genetic contribution to yield

gains of five major crops plants. CSSA-ASA. Special Publ. 7, Madison, WI. pp. 49-74.

Sumarno, and W.R. Fehr. 1982. Response to recurrent selection for yield in soybeans. Crop

Sci. 22: 295-299.

Sun, H., L. Zhao, J. Li, and S. Wang. 1999. The investigation of heterosis and pollen transfer

in soybean. In: H.E. Kauffman (ed.), World Soybean Res. Conf. VI. Superior Printing,

Champaign, IL. p. 489.

Thompson, J.A., and R.L. Nelson. 1998. Utilization of diverse germplasm for improving

soybean yield. Crop Sci. 38: 1362-1368.

Uphoff, M.D., W.R. Fehr, and S.R. Cianzio. 1997. Genetic gain for soybean seed yield by

three recurrent selection methods. Crop Sci. 37: 1155-1158. 15

Vaishampayan S.M., M. Kogan, G.P. Waldbauer, and J.T. Woolley. 1975. Spectral specific

responses in the visual behavior of the greenhouse whitefly, Trialeurodes vaporariorum

(Homoptera: Aleyrodidae). Ent. Exp. Appl. 18:344–356.

Voigt, R. L., and C.R. Weber. 1960. Effectiveness of selection methods for yield in soybean

crosses. Agron J. 52: 527-530.

Webb, D.M. 1997. Brown stem rot resistance in soybeans. Pioneer Hi-Bred International,

Inc., Johnston, Iowa, USA:534,091 (U.S. Patent No. 5,689,035), 09-26-1995.

Wilcox, J.R. 2001. Sixty years of improvement in publicly developed elite soybean lines.

Crop Sci. 41: 1711-1716.

Wrather, J.A, W.C. Stienstra, and S.R. Koenning. 2001. Soybean disease loss estimates for

the United States from 1996 to 1998. Can. J. Plant Pathol. 23: 122–131. 16

CHAPTER 2: GENETIC ANALYSIS OF NEW SOURCES OF SOYBEAN

RESISTANCE TO BROWN STEM ROT

A paper submitted to Crop Science

Paola T. Perez, Brian W. Diers, Peter Lundeen, Girma M. Tabor, and Silvia R. Cianzio*

P.T. Perez, P. Lundeen, and S.R. Cianzio, Dep. of Agronomy, Iowa State University, Ames, IA 50011-1010; B.W. Diers, Dep. of Crop Sciences, Univ. of Illinois, Urbana, IL 61801; G.M. Tabor, Pioneer Hi-Bred International, Inc. Johnston, IA 50131. *Corresponding author ([email protected])

ABSTRACT

Brown stem rot (BSR) of soybean [Glycine max (L.) Merr.], caused by Phialophora gregata, is an economically important disease prevalent in soybean producing regions of the north- central U.S. and Canada. To date, all BSR resistant genes identified are located on chromosome 16 (formerly molecular linkage group J). The objective of this study was to determine if four plant introductions (PI) from south-central China identified as BSR- resistant have resistance genes mapping to the same location on chromosome 16 as previously mapped BSR resistance genes. The four plant introductions, PI 594637, PI

594638B, PI 594650A, and PI 594858B were crossed to the BSR-susceptible cultivar

‘Century 84’ to develop four F2 populations. Each segregating population and the parental lines were screened for BSR resistance in growth chamber conditions. The F2:3 individual plants of each population were tested with the simple sequence repeat (SSR) markers Satt431 17

or Satt547, which map closely to BSR resistance quantitative trait loci (QTL) on chromosome 16. Associations between molecular data and phenotypic data used to validate

QTL were analyzed using single factor analysis of variance. Three of the four populations had markers on chromosome 16 significantly associated with BSR resistance with R2 values from 24 to 48%. However, when marker Satt547 was regressed on BSR resistance in population PI 594637 X Century 84, no significant association was observed. This result suggests that PI 594637 could have a new BSR resistance gene. Transgressive segregation also was observed in this population, and highly BSR resistant progeny could be used in the development of BSR resistant cultivars. Additional research and testing in this population will be conducted to identify resistance QTL(s) from this source.

Abreviations: BSR, brown stem rot • LSMEANS, Least square means • MAS, marker assisted selection • MG, maturity group • PCR, polymerase chain reaction • PI, plant introduction • QTL, quantitative trait loci • SSR, simple sequence repeat.

INTRODUCTION

Brown stem rot (BSR) of soybean is an economically important disease prevalent in soybean producing regions in the north-central U.S. and Canada (Hartman et al., 1999; Sinclair and

Backman, 1989). Under conditions that favor disease development, yield loses of up to 38% have been reported (Bachman et al., 1997; Gray, 1972). Phialophora gregata f. sp. sojae is the causal organism, infecting plants through the roots, colonizing the pith and vascular 18

system and moving up the stem and into the leaves of susceptible genotypes (Allington and

Chamberlain, 1948; Schneider et al., 1972). Two types of Phialophora gregata f. sp. sojae have been identified. Type I, causes leaf necrosis and defoliation, and stem vascular discoloration, and Type II, only causes vascular discoloration (Gray, 1971; Harrington et al.,

2003). The most effective way to manage the disease is by planting BSR-resistant soybean cultivars (Bachman et al., 1997; Mengistu et al., 1986).

Resistant genes have been identified and utilized in cultivar development and germplasm enhancement (Bachman et al., 2001; Chamberlain and Bernard, 1968; Tachibana and Card, 1972). Genetic analyses indicated that BSR resistance is conferred by dominant alleles at three independent loci, Rbs1, Rbs2 and Rbs3. The Rbs1 gene was identified in PI

84946-2 and L78-4094 (Bachman et al., 2001; Hanson et al., 1988; Klos et al., 2000; Lewers, et al., 1999; Sebastian et al., 1985), Rbs2 in PI 437833 (Bachman et al., 2001; Hanson et al.,

1988; Klos et al., 2000; Patzoldt et al., 2005b), and Rbs3 in PI 437970, PI 84946-2, 'BSR

101', 'BSR 201', 'BSR 301', 'BSR 302' and 'IA 2008' (Klos et al., 2000; Lewers, et al., 1999;

Willmot and Nickell, 1989). Other resistance sources have been identified, but their genetic control has yet to be determined (Lewers, et al., 1999; Nelson et al., 1989).

BSR disease expression exhibits a wide range of continuous variation in segregating populations. This may be an indication of environmental effects on the phenotypic disease expression and possible evidence of polygenic inheritance of BSR resistance (Bachman and

Nickell, 2000; Bachman et al., 1997), in addition to the major genes identified. These characteristics lead to low heritability of resistance, making development of high-yielding resistant cultivars difficult. Marker assisted selection (MAS) might be an alternative for 19

selection of traits such as BSR resistance. Once precise and accurate screening methods are achieved, environmental effects would less influence selection for BSR resistance.

To identify quantitative trait loci (QTL) that confer BSR resistance, mapping studies with various types of molecular markers have been conducted. Webb (1997) reported that the

Rbs3 gene was associated with BSR resistance on molecular linkage group J (chromosome

16, Schumtz et al., 2010) using restriction fragment length polymorphism (RFLP) markers on

328 recombinant inbred lines (RIL) from a cross between the resistant cultivar BSR 101 and the susceptible accession PI 437654. Lewers et al. (1999) identified one major and one minor

QTL on chromosome 16 responsible for BSR resistance in BSR101 using RFLP and amplified fragment length polymorphism (AFLP) markers. Klos et al. (2000) evaluated 17 resistant and 29 susceptible cultivars and plant introductions (PI) as well as RIL derived from the cross of BSR 101 with PI 437654. These authors used nine new DNA markers derived from RFLP, SSR, and bacterial artificial chromosome sequences, which mapped to chromosome 16. Bachman et al. (2001) observed in a population of F2:3 lines from the cross of L78-4094, a resistant line with the Rbs1 gene, to the susceptible cultivar Century, SSR markers Satt215 and Satt431 were linked to Rbs1. In the same study, the authors also identified SSR markers Satt244 and Satt431 linked to Rbs2 in a F2:3 population developed from crossing PI 437833 to Century. Patzoldt et al. (2005b) localized a BSR resistant QTL from ‘Bell’, a BSR resistant cultivar, close to molecular markers Satt431, Satt547,

21E22.sp1, 21E22.sp2, K375, and 35E22.sp1. These studies placed the BSR resistance genes

Rbs1, Rbs2, Rbs3 and a gene from Bell to the same region on chromosome 16.

New sources of resistance are important in breeding soybeans due to its narrow genetic 20

base (Gizlice et al., 1994). Plant introductions have been particularly valuable to soybean improvement because they can contribute different genes than those found in adapted genotypes (Carter et al., 2004; Patzoldt et al., 2005a; Thompson and Nelson, 1998; Thorne and Fehr, 1970; Vello et al., 1984). In an effort to expand the genetic base of resistance to

BSR, Patzoldt et al. (2005a) mapped BSR resistance QTL in five accessions from central

China. In all five accessions, they mapped a QTL for BSR resistance to the same chromosome 16 region as previously reported for the three BSR resistant genes. Our study continues the search in soybean for sources of new BSR resistance genes. Our objective was to determine if four BSR resistant plant introductions from south-central China had resistance genes different from those previously reported on chromosome 16.

MATERIALS AND METHODS

Plant material

Four plant introductions (PI) from south-central China previously identified as BSR- resistant (Patzoldt et al., 2003) were crossed to the maturity group (MG) II susceptible cultivar Century 84 (Walker et al., 1986). PIs were PI 594637, PI 594638B, PI 594650A of

MG IV, and PI 594858B from MG V.

The F1 seed was produced in Urbana, IL, and the F2 :3 mapping populations were developed at Isabela, PR. The F2:3 lines along with the parental lines were evaluated for BSR resistance in Ames, IA. One hundred and fifty five F2:3 individuals from each segregating population and their parental lines were grown in a growth chamber in a randomized 21

complete block design with three replications. One seed of each individual genotype was planted in a mix of soil, sand, and perlite in 4-by-21-cm plastic, cone-shaped containers.

Temperature in the growth chamber was constant at 19oC, with 16 hours of light. Two-week old plants were inoculated with P. gregata by introducing inoculum (described below), into the stems.

Inoculation protocol

Inoculum of P. gregata f. sp. sojae isolate Oh2-3, was prepared in cultures as described by Tabor et al. (2003). Cultures were started on green bean agar (GBE) medium

(ground frozen Phaseolus vulgaris L. green pods at 35 g/liter) supplemented with ampicillin at 50 mg/liter and were incubated for 44 to 64 days at room temperature (21 to 23oC) in the dark until abundant sporulation was observed. Mycelium of P. gregata was harvested from the agar plates and conidia were suspended in 0.8% water agar (2.7 x 107 conidia/ml). The conidial suspension was mixed into a paste, then, the bevel of an 18-gauge needle was filled with this inoculum paste. Stems of two-week-old plants were punctured approximately 2 cm above the soil line to introduce inocula. After inoculation, plants were grown under 16-h light-photoperiod, at a constant temperature of 19 oC, watered daily, and fertilized weekly.

Severity assays

Five weeks after inoculation, BSR severity was evaluated using both a BSR severity or symptom scale and a measurement of the recovery of P. gregata from the stem after plating stem pieces on green bean agar (Tabor et al., 2003). For the first method, disease severity was accessed by evaluating each plant for its vegetative stage (V-stage, Fehr et al.,

1971) and for its vigor (overall plant health condition) (Table 1). For each individual plant in 22

the three replications, vigor was measured on a scale from 1 to 7. In this scale, a score of 1 corresponds to a dead plant; 2 to a plant with a green stem and no leaves; 3 when chlorotic and necrotic leaves are prominent; 4 when some stunting, mosaic chlorosis and necrosis on leaves is observed; 5 when leaf area is normal except for some yellowing; 6 when leaf area is normal, and plants are small but healthy; and 7 when plants are completely healthy (P.

Lundeen, unpublished).

For the plating method, the procedure described by Tabor et al. (2003) was used.

After vigor scores were taken, plants were defoliated and stems were cut at the soil line, immersed for 3 min in 70% ethanol, followed by 5 min in 10% sodium hypoclorite and finally rinsed with sterile water. Incidence of BSR was estimated by cutting 1-cm stem segments at the inoculation point and at the top of the plants. The two stem segments were plated on green bean agar, supplemented with ampicillin. After plating stem segments, cultures were allowed to grow in the dark at 5 oC for two weeks. Following this period, plates were evaluated for P. gregata infection. A plant was considered colonized if P. gregata was recovered from the stem segment corresponding to the inoculation point. Plants were considered susceptible if P. gregata also was recovered from the top portion of the stem in addition to the inoculation point. Conversely, a plant was considered resistant if the fungus was recovered only from the inoculation point and was not present at the top segment (Figure

1). The variable plate susceptibility was recorded as 0 for a resistant plant, which showed infection only at the inoculation point; and as 1 for a susceptible plant, infected at the inoculation point and plant top. 23

DNA marker analysis

DNA extraction and genotyping was performed in Urbana, IL. For each mapping population, DNA was extracted from leaves of 10 plants in each F2:3 line. DNA was extracted according to Keim and Shoemaker (1988) with modifications outlined by Patzoldt et al.

(2005b). Each line was genotyped with SSR markers Satt431 or Satt547 from chromosome

16 (Song et al., 2004) according to Patzoldt, et al. (2005b). SSR markers were amplified through polymerase chain reaction (PCR). Amplified samples were run on 6% nondenaturing acrylamide (Wang et al., 2003) or 3% Metaphor agarose gels (Patzoldt et al., 2005b).

Ethidium bromide was used to stain the gels and bands were visualized under UV light.

Data analysis

Phenotypic data were analyzed using PROC MIXED of SAS v. 9.1 (SAS Institute,

2 2003). Broad-sense heritabilities were estimated for each F2 population with the formula H =

2 2 2 σ G / ((σ E/r) + σ G), where r is the number of replications. All effects were treated as random in the model and variance components were estimated with the COVTEST option.

Subsequently, F2:3 lines and parental lines were considered fixed effects. Least square means

(LSMEANS) were used to estimate average BSR resistance level for each F2:3 individual and parental lines. The ESTIMATE statement was used to test the difference between average of the F2:3 individuals and the average of each parental line. Transgressive segregation was estimated according to the method used by Shahid et al. (2008). The BSR resistance level of each F2:3 line was compared with the resistance exhibited by each parent. Transgressive segregants for resistance were defined as F2:3 lines which exceed the most resistant parent by two standard deviations. Similarly, the number of F2:3 lines with transgressive segregation for 24

BSR susceptibility was calculated as the number of F2:3 individuals with susceptibility higher than that of the most susceptible parent by two standard deviations. Single-factor analysis of variance using PROC GLM of SAS v.9.1 was performed to detect associations between SSR markers and BSR resistance in segregating populations. A significant association was declared if P<0.05.

RESULTS

Phenotypic analysis

The broad-sense heritability calculated for vigor ranged from 22% to 60%, and from

35% to 70% for plate susceptibility (data not shown). Each segregating population, along with the susceptible and resistant parents, was evaluated for plate susceptibility and vigor to

BSR infection (Table 2). The resistant parents, PI 594638B, PI 594650A, and PI 594658B, showed zero plate susceptibility, meaning that P. gregata was not recovered from the top of the plants. The PIs also showed high vigor. In contrast, PI 594637 exhibited greater BSR susceptibility and was not statistically different from Century 84. The susceptible parent,

Century 84, also showed susceptibility for plate susceptibility and vigor scores. In all F2 populations, segregation was observed for visual symptoms and plate susceptibility (Figures

2-3).

The average of F2:3 lines in the PI 594637 X Century 84 population showed intermediate values of BSR susceptibility, the average for plate susceptibility was 0.48 and

3.63 for vigor (Table 2). Average value of the F2:3 lines for plate susceptibility was not statistically different from both parental lines. However, the average of F2:3 lines showed 25

significantly higher average vigor than PI 594637 (Table 2; Figures 2-3). Although PI

594637 was used in these tests as a resistant parent based on previous screening (parents alone), when the screening was repeated for the current tests, it showed higher susceptibility to BSR symptoms than those of the F2 population. This population showed transgressive segregants in both directions (Table 3). There were 8 transgressive segregants two standard deviations more resistant than the most resistant parent, and 24 transgressive segregants two standard deviations more susceptible than the most susceptible parent.

Population PI 594638B X Century 84 showed segregation for the two resistance traits

(Table 2; Figures 2-3). F2:3 individuals exhibited averages of 0.34 and 5.54 for plate susceptibility and vigor, respectively (Table 2). The average of the F2:3 lines was significantly (P-value < 0.5) more resistant than Century 84 for plate susceptibility. No statistical differences were detected between the average of the F2:3 lines and PI 594638B for plate susceptibility and vigor. Only three transgressive segregants for BSR resistance were observed (Table 3).

For Population PI 594650A X Century 84, averages of F2:3 individuals were 0.26 and

5.15 for plate susceptibility and vigor, respectively (Table 2; Figure 2-3). The population average was significantly (P-value < 0.01) more BSR resistant than the susceptible parent

Century 84 for plate susceptibility and vigor (Table 2). Conversely, no significant differences were detected between the average of the F2:3 individuals and the resistant parent PI

594650A. Only one transgressive segregant was observed and it was for BSR susceptibility

(Table 3). 26

Population PI 594658B X Century 84 showed segregation for the two variables

(Table 2; Figure 2-3). Averages of plate susceptibility and vigor were 0.29 and 5.94, respectively. The average of the F2:3 lines was significantly (P-value <0.01) different from the susceptible parent, with lower values for plate susceptibility and higher values of vigor.

No statistical differences were detected between the mean of the F2:3 lines and PI 594638B for vigor. This population had 20 transgressive segregants more susceptible than the susceptible parent (Table 3).

Single marker analysis

Marker Satt431 showed polymorphism for population PI 594650A x Century 84, and

Satt547 was polymorphic for populations PI 594637 X Century 84, PI 594638B X Century

84, and PI 594658B X Century 84. Phenotypic data from each segregating population was used to test for associations between BSR resistance and molecular markers Satt431 and

Satt547, which are in the region where BSR resistance was previously mapped, using single marker analysis. In three of the four populations, F2:3 lines homozygous for the marker allele from the PI showed significantly greater BSR resistance than lines homozygous for alleles from Century 84 (Table 4). For the three populations in which significant associations were detected between the markers and the phenotypes, lines segregating for the marker exhibited average resistance levels that were intermediate between the two classes of homozygous lines.

Single marker regression revealed a significant association between BSR resistance and the SSR markers on chromosome 16 for three of the four populations (Table 4).

Molecular marker Satt547 was significantly associated with BSR resistance in population PI 27

594638B X Century 84, with R2 values of 44% for plate susceptibility and 33% for vigor. In population PI 594650A X Century 84, Satt431 was significantly associated with BSR resistance, with R2 values of 34 and 37% for plate susceptibility and vigor, respectively.

Satt547 was significantly associated to BSR resistance in population PI 594658B X Century

84, with R2 values of 48% for plate susceptibility and 25% for vigor. Conversely, population

PI 594637 X Century 84 did not show significant association with molecular marker Satt547.

These results are suggestive of a possible new source of resistance in PI 594637, mapping to a different region from those previously identified on chromosome 16, where Rsb1, Rsb2, and Rsb3 map.

DISCUSSION

The study was conducted to determine if BSR resistance from PIs from south-central

China is governed by genes mapping to genetic regions different from the previously reported resistance genes, Rbs1, Rbs2, and Rbs3. We hypothesized that association between

BSR resistance and SSR molecular makers Satt431 or Satt547, in F2 segregating populations, could indicate the presence of any of the previously identified resistance genes, which map close to these SSR markers on chromosome 16.

Given the high association between the SSR markers used and BSR resistance in populations PI 594638B X Century 84, PI 594650A X Century 84, and PI 594658B X

Century 84, the BSR resistance genes present in these PIs could be allelic to previously identified BSR resistance genes on chromosome 16. This is because the three BSR resistance genes, Rbs1, Rbs2, and Rbs3 were confirmed to map to a region on chromosome 28

16 that the SSR markers Satt431 and Satt537 map (Bachman et al., 2001; Klos et al., 2000;

Lewers et al., 1999; Patzold et al., 2005b; Webb, 1997). BSR resistance in population PI

594637 X Century 84 showed non-association with SSR Satt547. This could be indicative of the presence of a new resistance gene not associated with the molecular marker Satt547 on chromosome 16. However, in order to test this assumption, allelism tests between PI 594637 and genotypes possessing Rsb1, Rbs2, and Rbs3 genes need to be performed to confirm if PI

594637 carries a gene non allelic to any of the three genes. This research is currently in progress.

In contrast to the other three PIs used, which were highly resistant to BSR, PI 594637 showed intermediate resistance. Highly resistant individuals, however, were present in the

F2:3 population derived from this PI crossed with susceptible cultivar Century 84 (Table 3;

Figure 2-3), and transgressive segregant individuals for BSR resistance also were detected.

Transgressive segregation of traits is a common phenomenon in segregating hybrid populations, especially in plants, where transgressive segregants might be the rule rather than the exception (Rieseberg et al., 1999; Rieseberg et al., 2003). A number of different mechanisms that could be responsible for transgressive segregation in hybrids are: elevated mutation rate, reduced developmental stability, epistasis, overdominance, unmasking of rare recessive alleles, chromosome number variation, and complementary gene action (Lynch and

Walsh, 1998; Rieseberg et al., 1999; Rieseberg et al., 2003). In plants, results from classical and QTL studies have provided evidence to conclude that complementary gene action is the most likely explanation for transgressive segregation (Li et al., 1995; Vega and Frey, 1980; de Vicente and Tanksley, 1993). Vega and Frey (1980) working in oats concluded parents 29

might have contributed one or more useful genes to traits probably with a cumulative effect on phenotype. de Vicente and Tanksley (1993) in tomatoes, reported transgressive segregation was directly linked to the presence of complementary QTL alleles in the parental cultivars. Transgressive segregants in soybean also have been observed (Fehr et al., 1991;

Alt et al., 2005).

For BSR resistance in soybean, Bachman and Nickell (2000) proposed a genetic model in which resistance is the result of epistatic interaction between pairs of loci. It is plausible that the high levels of resistant observed in some F2:3 progeny of population PI

594637 X Century 84 could be due to epistatic or additive interaction between resistance loci from PI 594637 and Century 84. Although Century 84 was considered as the susceptible parent, some degree of resistance was also observed in it for the variables evaluated (Table

2). In the cross PI 594637 X Century 84, the highly BSR resistant progeny identified can be used for development of BSR resistant cultivars. This alternative could even be more attractive if the PI is confirmed as possessing a new source of resistance for BSR.

Acknowledgements

The research was in part financed by funds provided by the Iowa Soybean Association (ISA), and the North Central Soybean Research Program (NCSRP). We thank Dr. Charlotte

Bronson for the use of her laboratory facilities, and to Mr. Allen Vogl for his assistance in the use of growth chambers at the Dept. of Agronomy, Iowa State University.

30

REFERENCES

Allington, W.B., and D.W. Chamberlain. 1948. Brown stem rot of soybean. Phytopathology

38:793–802.

Alt, J.L., W.R. Fehr, G.A. Welke, and J.G. Shannon. 2005. Transgressive segregation for

oleate content in three soybean populations. Crop Sci. 45:2005-2007.

Bachman, M.S., C.D. Nickell, P.A. Stephens, and A.D. Nickell. 1997. Brown stem rot

resistance in soybean germplasm from central China. Plant Dis. 81:953–956.

Bachman, M.S., and C.D. Nickell. 2000. Investigating the genetic model for brown stem rot

resistance in soybean. J. Hered. 91:316–321.

Bachman, M.S., J.P. Tamulonis, C.D. Nickell, and A.F. Bent. 2001. Molecular markers

linked to brown stem rot resistance genes Rbs1 and Rbs2 in soybean. Crop Sci. 41:527–

535.

Carter Jr, T.E., R.L. Nelson, C. Sneller, and C. Zhanglin. 2004. Genetic diversity in

soybean. H.R. Boerma and J.E. Specht (eds), Soybean Monograph, 3rd Edition, Am.

Soc. of Agron., Madison, WI. 114 pp.

Chamberlain, D.W., and R.L. Bernard. 1968. Resistance to brown stem rot in soybeans.

Crop Sci. 8:728–729.

Fehr, W.R., C.E. Caviness, D.T. Burmood, and J.S. Pennington. 1971. Stage of development

descriptions for soybeans, Glycine max (L.) Merrill. Crop Sci. 11:929– 931.

Fehr, W.R., G.A. Welke, E.G. Hammond, D.N. Duvick, and S.R. Cianzio. 1991. Inheritance

of elevated palmitic acid content in soybean seed oil. Crop Sci. 31:1522–1524. 31

Gizlice Z., T.E. Carter, Jr., and J.W. Burton. 1994. Genetic base for North American soybean

cultivars released between 1947 and 1988. Crop Sci. 34:1143-1151.

Gray, L.E. 1971. Variation in pathogenicity of Cephalosporium gregatum isolates.

Phytopathology 61:1410-1411.

Gray, L.E. 1972. Effect of Cephalosporium gregatum on soybean yield. Plant Dis. Rep.

56:580-581.

Hanson, P.M., C.D. Nickell, L.E. Gray, and S.A. Sebastian. 1988. Identification of two

dominant genes conditioning brown stem rot resistance in soybean. Crop Sci. 28:41–43.

Hartman, G.L., J.B. Sinclair, and J.C. Rupe. 1999. Compendium of soybean diseases (4th

edition). The American Phytopathological Society, St. Paul, Minnesota, USA.

Harrington, T.C., J. Steimel, F. Workneh, and X.B. Yang. 2003. Characterization and

distribution of two races of Phialophora gregata in the north-central United States.

Phytopathology 93:901-912.

Keim, P., and R.C. Shoemaker. 1988. Construction of a random recombinant DNA library

that is primarily single copy sequence. Soybean Genet. Newsl. 15:147-148.

Klos, K.L.E., M.M. Paz, L.F. Marek, P.B. Cregan, and R.C. Shoemaker. 2000. Molecular

markers useful for detecting resistance to brown stem rot in soybean. Crop Sci.

40:1445–1452.

Lewers, K.S., E.H. Crane, C.R. Bronson, J.M. Schupp, P. Keim, and R.C. Shoemaker.

1999. Detection of linked QTL for soybean brown stem rot resistance in ‘BSR 101’ as

expressed in a growth chamber environment. Molec. Breed. 5:33–42. 32

Li, Z., S.R.M. Pinson, J.W. Stansel, and W.D. Park. 1995. Identification of quantitative trait

loci (QTLs) for heading date and plant height in cultivated rice (Oryza sativa L.).

Theor. Appl. Genet. 91:374-381.

Lynch M., and B. Walsh. 1998. Genetics and analysis of quantitative traits. Sinauer

Associates, Sunderland, MA.

Mengistu, A., C.R. Grau, and E.T. Gritton. 1986. Comparison of soybean genotypes for

resistance to and agronomic performance in the presence of Brow Stem Rot. Plant Dis.

70: 1095-1098.

Nelson, R.L., C.D. Nickell, J.H. Orf, H. Tachibana, E.T. Gritton, C.R. Grau, and B.W.

Kennedy.1989. Evaluating soybean germplasm for brown stem rot resistance. Plant

Dis. 73:110-114.

Patzoldt, M.E., W. Chen, and B.W. Diers. 2003. Evaluation of soybean plant introductions

from China for resistance to brown stem rot [Online]. Available at

http://www.plantmanagementnetwork.org/sub/php/research/2003/bsr. Plant Health

Prog.

Patzoldt, M. E., S. R. Carlson, and B.W. Diers. 2005a. Characterization of resistance to

brown stem rot of soybean in five accessions from central China. Crop Sci. 45:1092-

1095.

Patzoldt, M.E., C.R. Grau, P.A. Stephens, N.C. Kurtzweil, S.R. Carlson, and B.W. Diers.

2005b. Localization of a quantitative trait locus providing brown stem rot resistance in

the soybean cultivar Bell. Crop Sci. 45:1241-1248. 33

Rieseberg, L.H., M.A. Archer, and R.K. Wayne. 1999. Transgressive segregation,

adaptation and speciation. Heredity 83:363-372.

Rieseberg L.H., A. Widmer, A.M. Arntz, and J.M. Burke. 2003. The genetic architecture

necessary for transgressive segregation is common in both natural and domesticated

populations. Philos. T. R. Soc. B. 358:1141-1147.

SAS System. Version 9.1. 2003. Cary, NC: SAS Institute Inc.

Schumtz, J., S.B. Cannon, J. Schlueter, J. Ma, T. Mitros, W. Nelson, D.L. Hyten, Q. Song,

J.J. Thelen, J. Cheng, D. Xu, U. Hellsten, G.D. May, Y. Yu, T. Sakurai, T. Umezawa,

M.K. Bhattacharyya, D. Sandhu, B. Valliyodan, E. Lindquist, M. Peto, D. Grant, S.

Shu, D. Goodstein, K. Barry, M. Futrell-Griggs, B. Abernathy, J. Du, Z. Tian, L Zhu, N.

Gill, T. Joshi, M. Libault, A. Sethuraman, X.-C. Zhang, K. Shinozaki, H.T. Nguyen, R.A.

Wing, P. Cregan, J. Specht, J. Grimwood, D. Rokhsar, G. Stacey, R.C.

Shoemaker, and S.A. Jackson. 2010. Genome sequence of the paleopolyploid soybean.

Nature 463: 178–183.

Schneider, R.W., J.B. Sinclair, and L.E. Gray. 1972. Etiology of Cephalosporium gragatum

in soybean. Phytopathology 62:345-349.

Sebastian, S.A., C.D. Nickell, and L.E. Gray. 1985. Efficient selection for brown stem rot

resistance in soybeans under greenhouse screening conditions. Crop Sci. 25:753–757.

Shahid, M., S. Han, H. Yoell, and J. Xu. 2008. Fitness distribution and transgressive

segregation across 40 environments in a hybrid progeny population of the human-

pathogenic yeast Cryptococcus neoformans. Genome 51:272-281. 34

Sinclair, J.B., and P.A. Backman. 1989. A compendium of soybean diseases (3rd edition).

American Phytopathological Society, St. Paul, Minnesota, USA.

Song Q.J., L.F. Marek, R.C. Shoemaker, K.G. Lark, V.C. Concibido, X. Delannay, J.E.

Specht, and P.B. Cregan. 2004. A new integrated genetic linkage map of the soybean.

Theor. Appl. Genet. 109:122-128.

Tabor, G.M., G.L. Tylka, S.R. Cianzio, and C.R. Bronson. 2003. Resistance to Phialophora

gregata is expressed in the stems of resistant soybeans. Plant Disease 87:970-976.

Tachibana, H., and L.C. Card. 1972. Brown stem rot resistance and its modification by

soybean mosaic virus in soybeans. Phytopathology 62:1314-1317.

Thompson J.A., and R.L. Nelson. 1998. Utilization of diverse germplasm for soybean yield

improvement. Crop Sci. 38:1348-1355.

Thorne J.C., and W.R. Fehr. 1970. Incorporation of high-protein, exotic germplasm into

soybean populations by 2- and 3-way crosses. Crop Sci. 10:652-655.

Vega, U., and K.J. Frey. 1980. Transgresssive segregation in inter and intraspecific crosses

of barley. Euphytica 29:585-594.

Vello, N.A., W.R. Fehr, and J.B. Bahrenfus. 1984. Genetic variability and agronomic

performance of soybean populations developed from plant introductions. Crop Sci.

24:511-514. de Vicente, M. C. and S.D. Tanksley. 1993. QTL analysis of transgressive segregation in an

interspecific tomato cross. Genetics 134:585-596. 35

Walker, A.K., A.F. Schmitthenner, R.J. Fioritto, S.K. St. Martin, R.L. Cooper, and R.J.

Martin. 1986. Registration of ‘Century 84’ soybean. Crop Sci. 26:199–200.

Wang, D., J. Shi, S.R. Carlson, P.B. Cregan, R.W. Ward, and B.W. Diers. 2003. A low-

cost, high-throughput polyacrylamide gel electrophoresis system for genotyping with

microsatellite DNA markers. Crop Sci. 43:1828–1832.

Webb, D.M. 1997. Brown stem rot resistance in soybeans. Pioneer Hi-Bred International,

Inc., Johnston, Iowa, USA:534,091 (U.S. Patent No. 5,689,035), 09-26-1995.

Willmot, D.B., and C.D. Nickell. 1989. Genetic analysis of brown stem rot resistance in

soybean. Crop Sci. 29:672–674.

36

Table 1. Visual assessment scale of brown stem rot (BSR) severity in soybean plants five weeks after inoculation.

Vigor† Plant Description 1 Dead plant 2 Green stem and no leaves 3 Chlorotic and necrotic leaves are prominent 4 Some stunting, mosaic chlorosis and necrosis on leaves 5 Leaf area is normal, with some yellowing 6 Leaf area is normal, plants are small but healthy 7 Healthy plants. † Vigor is a measure of overall plant health.

37

Table 2. Brown stem rot (BSR) resistance ratings for plate susceptibility† and vigor‡ in four F2 populations and their parental lines. P-value of the difference

between F population Mean values 2 Variable and parental lines F Century Mean F vs. Mean F 2 PI 2 2 population 84 Century 84 vs. PI PI 594637 X Century 84 Plate susceptibility 0.48 0.45 0.47 0.945 0.959 Vigor 3.63 2.95 3.47 0.527 0.072

PI 594638B X Century 84 Plate susceptibility 0.34 0.00 0.80 0.049 0.591 Vigor 5.54 5.71 4.32 0.137 0.846

PI 594650A X Century 84 Plate susceptibility 0.26 0.00 0.83 <0.0001 0.365 Vigor 5.15 5.87 3.44 <0.0001 0.172

PI 594658B X Century 84 Plate susceptibility 0.29 0.00 1.00 0.0001 0.097 Vigor 5.94 6.57 4.25 0.003 0.193 † Plate susceptibility: 1 = infection at inoculation point and at plant top; 0 = infection only at inoculation point. ‡ Vigor: Overall plant health, measured on a scale from 1 to 7; 1= dead plant; 7 = healthy plant. 38

Table 3. Number of F2:3 lines with higher brown stem rot (BSR) resistance than the most resistant parent or lower than the most susceptible parent by two standard deviations.

Number of F2:3 lines showing transgressive segregation Transgressive segregtion Plate Vigor Total Susceptibility PI 594637 X Century 84 For resistance† 2 6 8 For susceptibility‡ 23 1 24

PI 594638B X Century 84 For resistance 0 3 3 For susceptibility 0 0 0

PI 594650A X Century 84 For resistance 0 0 0 For susceptibility 1 0 1

PI 594658B X Century 84 For resistance 0 0 0 For susceptibility 20 0 20 † Transgressive segregants for BSR resistance. Lines with resistance higher than that of the most resistant parent by 2σ. ‡ Transgressive segregants for BSR susceptibility. Lines with susceptibility higher than that of the most susceptible parent by 2σ. 39

Table 4. Means of genotypic classes and R2 values of molecular markers on chromosome 16 (formerly molecular linkage group J) associated with brown stem rot (BSR) resistance.

Mean values for genotypic classes SSR Homozygous Homozygous Segregating R2 Variable marker for alleles for alleles for alleles (%) from from from Century from both PI 84 parents PI 594637 X Century 84 Satt547 Plate susceptibility 0.50 0.47 0.46 2.1 Vigor 3.61 3.67 3.62 0.1

PI 594638B X Century 84 Satt547 Plate susceptibility 0.05 0.70 0.31 44.4*** Vigor 6.51 4.56 5.63 33.0***

PI 594650A X Century 84 Satt431 Plate susceptibility 0.06 0.67 0.27 34.3*** Vigor 5.89 3.80 5.10 37.0***

PI 594658B X Century 84 Satt547 Plate susceptibility 0.03 0.66 0.21 48.1*** Vigor 6.47 5.18 6.09 25.0*** *** Significant at the 0.001 probability level. 40

A B

Figure 1. Plate susceptibility. A) Susceptible plant with infection of P. gregata in stem fragments from the inoculation point and from the plant top. B) Resistant plant with infection only at the inoculation point.

41

Figure 2. Distribution of plate susceptibility (1 = infection at inoculation point and at plant top; 0 = infection only at inoculation point) for four segregating populations and their parental lines. Arrows indicate parental mean values (continuous arrows for Century 84 and dashed arrows for PI 594637, PI 594638B, PI 594650A, and PI 594658B). 42

Figure 3. Distribution of vigor (measured on a scale from 1 to 7; 1= dead plant; 7 = healthy plant) for four segregating populations and their parental lines. Arrows indicate parental mean values (continuous arrows for Century 84 and dashed arrows for PI 594637, PI 594638B, PI 594650A, and PI 594658B).

43

CHAPTER 3: QTL MAPPING OF WHITEFLY RESISTANCE IN SOYBEAN

A paper to be submitted to Journal of Crop Improvement

Paola T. Perez,*1 Silvia R. Cianzio1, Pinar C. Kara2, Mayra Aviles3, and Reid G. Palmer4

1P.T. Perez and S.R. Cianzio, Dep. of Agronomy, Iowa State University, Ames, IA 50011-

1010 USA; 2 P.C. Kara, Cukurova Agricultural Research Institute, Adana, Turkey; 3 Mayra

Aviles, International Maize and Wheat Improvement Center (CYMMIT) - Ciudad Obregon,

Sonora, Mexico; 4 R.G..Palmer, USDA-ARS-CICGR, and Dep. of Agronomy, Iowa State

University, Ames, IA 50011-1010 USA. *Corresponding author ([email protected])

ABSTRACT

Whitefly Bemisia tabaci (Gennadius) is one of the most damaging insects attacking crops in the world. In soybean [Glycine max (L.) Merr.], in addition to causing plant stress and reduced seed size, they also can be vectors of viruses, e.g. soybean crinkle mosaic, and soybean dwarf mosaic. Resistance to whitefly in soybean has been reported; however, whitefly resistance genes have not been identified. The objective of the study was to identify quantitative trait loci (QTL) for whitefly resistance in soybean. Two F2 segregating populations, derived from crosses between whitefly resistant lines ‘Cajeme’ and 'Corsoy 79' and the susceptible line ‘Williams 79’, were developed to determine inheritance of resistance to B. tabaci. Parental lines were screened with simple sequence repeats (SSR). Segregating

F2 populations and parents were evaluated for whitefly infestation for two years at one 44

location in Northwest Mexico, where whitefly is a recurrent pest of soybean. QTL analyses were performed using standard interval mapping (SIM), and composite interval mapping

(CIM). Results indicate that whitefly resistance is a polygenic trait, controlled by QTL on different chromosomes. The significant QTL detected in this study were located on chromosomes 12, 18, and 19 (formerly molecular linkage groups H, G, and L, respectively).

In addition, there was evidence that other QTL could be involved in resistance, such as those detected on chromosomes 1, 7, 10, 16, and 17 (formerly molecular linkage groups D1a, M,

O, J, and D2, respectively). Fine-mapping of the large effect QTL detected in this study would be helpful to identify tightly linked markers for use in marker-assisted selection

(MAS) for whitefly resistance in soybean.

Abbreviations: • CIM, composite interval mapping • LOD, likelihood of odds • LSMEANS, Least square means • MAS, marker-assisted selection • PCR, polymerase chain reaction • QTL, quantitative trait loci • SIM, standard interval mapping • SSR, simple sequence repeat.

INTRODUCTION

Whitefly, Bemisia tabaci (Gennadius) (order Homoptera, family Aleyrodidae), is a common pest in many crops throughout the world. Historically, whiteflies have been major pests in vegetable crops, cotton, and greenhouse plants. In addition, they have a wide range of hosts, including wild plant species, which makes it difficult to control once they have infested a crop (Perring, 2001). In the United States, annual yield losses exceed US $200 million, and occur in cotton, peanuts, soybean, vegetables, and ornamentals (Faust, 1992).

Damage can be caused by whiteflies directly feeding on the crop and also by virus 45

transmission while feeding (Costa, 1976; Kogan and Herzog, 1980; Vaishampayan et al.

1975).

Infestations of whiteflies in soybeans have been recorded in Brazil, India, Japan,

Turkey, Australia, and Mexico (Arioglu, 1987; Costa, 1976; Vaishampayan et al. 1975). In certain regions of Mexico, soybean plantings have been banned due to heavy whitefly infestations that could also spread to other crop species (Morales and Anderson, 2001). In

Puerto Rico, whiteflies have caused losses to soybean planted in off-season nurseries. In the

U.S., losses have been reported in Florida and Georgia (McPherson and Douce, 1992).

Recently, the insect has been reported in soybean fields in Illinois (Gray, 1999) and in Iowa

(O’Neal and Rice, 2006; Rice, 2007).

Whitefly, B.tabaci, is a species complex formed by highly cryptic sibling species

(Campbell, 1993; Campbell et al., 1996; Perring, 2001). To date, 41 distinct populations of B. tabaci have been studied; 24 of these populations have been given a specific biotype designation (Perring, 2001). Adults and nymphs have piercing-sucking mouthparts and feed on the lower leaf surfaces where they lay their eggs (Byrne and Bellons, 1991; Perring,

2001). The life cycle averages 32 days, depending on temperature, i.e. temperatures of about

25 ºC favor development of the insect on soybean (University of Illinois - Integrated Pest

Management Bulletin, 2004). Whiteflies damage crops by extracting large quantities of phloem sap (Jones, 2003), and heavy infestations may result in development of chlorotic spots on leaves, wilting, and stunting of plants. In addition, the insects excrete a sticky material called honeydew, which in high concentrations promotes the growth of sooty mold fungi, which interferes with the host plant photosynthesis. Some species of the whitefly 46

complex also serve as vectors of several economically important viral plant pathogens (Byrne and Bellons, 1991). Their high reproduction rate and short generation interval during the growing season favor development of resistance to insecticides. Under this scenario, management of the pest is complex and chemical controls are usually ineffective. Breeding resistant cultivars is the most effective management alternative.

In soybean, viruses such as soybean crinkle mosaic, and soybean dwarf mosaic can be transmitted by whiteflies (Costa, 1976). Two biotypes of whiteflies have been recorded that colonize field-grown soybean: the sweetpotato whitefly, B. tabaci, and the bandedwinged whitefly, Trialeurodes abutiloneus (Haldeman) (Vaishampayan et al. 1975). In greenhouse conditions, soybean is colonized by a different species, Trialeurodes vaporariorum

(Westwood).

In the search for genetic resistance that could be used in breeding soybean, systematic screenings of germplasm collections have been done. In 1987, 109 soybean cultivars were screened in Turkey (Arioglu, 1987) for resistance to B. tabaci during the pod-filling stage.

Variation among genotypes was observed, with the indication that one factor that could be affecting whitefly resistance was pubescence type (long, medium, short, acute, prostrate, erect) and density (glabrous, sparse, dense). Soybean genotypes with medium to long pubescence have, in general, lower whitefly infestation. Conversely, cultivars with dense pubescence have higher whitefly infestation than glabrous cultivars (Arioglu et al., 1989;

Kilen and Lambert, 1993; Lambert et al., 1995; McAuslane, 1996; McAuslane et al., 1995).

Partial whitefly resistance of glabrous soybean cultivars is due in part, by reduced ovipositional preference related to lack of foliar pubescence (McAuslane, 1996). 47

The wide adaptation of the insect and its increase in spread into new regions are of utmost concern. This is especially true due to the limited information about whitefly on soybean in the Midwest, and the lack of information about the inheritance of resistance to whitefly in soybean. To date few inheritance studies of whitefly resistance in soybean have been reported. Xu et al. (2009) studied inheritance of whitefly resistance by analyzing different genetic models. The authors evaluated mean nymph whitefly infestation per leaflet of F1 and F2:3 plant populations from the cross between a highly susceptible variety 'Qihuang' and a highly resistant variety 'Huapidou'. They found that the model with a major gene + polygene(s) had the best fit to explain whitefly resistance in soybean. These findings indicate the complex inheritance of whitefly resistance in soybean. In the study, it was not determined, at the molecular level, how many and where the quantitative trait loci (QTL) were located. For a trait of such complexity, incorporation of molecular marker information in inheritance studies is an important requirement to identify and locate the resistance genes that could be linked to molecular markers. The objective of our study was therefore, to identify number and location of QTL for whitefly resistance in soybean.

MATERIALS AND METHODS

Population development

Two F2 segregating populations for whitefly resistance were developed. On the basis of Arioglu’s (1987) results, two cultivars 'Corsoy 79' (Bernard and Cremmens, 1988a) and

'Cajeme', a Mexican cultivar, were identified as whitefly resistant. The cultivars were crossed to the susceptible cultivar 'Williams 79' (Bernard and Cremmens, 1988b), which was used as female in both crosses. 48

The mapping populations Williams 79 x Corsoy 79, and Williams 79 x Cajeme consisted of 181 and 156 F2 plants, respectively. The F1 seed was obtained at Ames, IA, and was sent to Isabela, Puerto Rico, where the F2 plants and their F2:3 families were developed.

Leaves were collected from F2 plants for DNA extraction, which afterwards were grown to maturity. Each F2 plant was single-plant threshed to generate the F2:3 families. From the

Williams 79 x Corsoy 79 population, 150 F2:3 families were obtained, and from the Williams

79 x Cajeme population, 90 F2:3 families.

Phenotypic evaluation

Phenotyping of F2:3 families was done in Isabela, Puerto Rico in 2001, however, the environmental conditions that year did not favor whitefly infestation and use of this location was discontinued. In 2003 and 2004 phenotypic evaluations of F2:3 families were done at the

International Maize and Wheat Improvement Center (CYMMIT) - Ciudad Obregon, Sonora,

Mexico. The location was selected because whitefly is a common and recurrent pest of soybean in this region. Parents and F2:3 lines were assigned to plots, arranged in a randomized complete block design with three replications, with each population grown as a separate test. The lines were grown in 4-row plots, 3 m long, 0.6 m between rows, planted at

8-seeds/0.3048 m.

Assessment of whitefly infestation was done weekly (7-10 times) during the pod- filling period, which is usually when the heaviest infestation occurs in soybean (Arioglu,

1987). From five random plants in each plot, five fully expanded leaflets were cut from near the top of the plant and the number of nymphs in an area of 6.5 cm2 (1 square inch) was 49

recorded using a magnifying glass. This number represented the variable nymph density.

Nymph density at the peak of whitefly infestation was used for QTL analysis.

DNA marker analysis

From each F2 population and the parental lines, genomic DNA was extracted from leaf tissue. Both populations were genotyped with polymorphic simple sequence repeat

(SSR) markers covering the 20 chromosomes (Schmutz et al., 2010) of the soybean genome, according to the soybean consensus map (Cregan et al., 1999; Song et al., 2004). Genotyping was performed by the Monsanto molecular laboratory in Huxley, IA. The SSR marker fragments were amplified in a 10 µl reaction mixture consisting of 10x polymerase chain reaction (PCR) buffer, 1 µl containing 1.25 mM of each dNTP, 1 µl forward and reverse primers (20 uM), 20 ng DNA, 0.3 µl Taq DNA polymerase (5 U/µl), and sterile deionized

H2O to volume. Amplification conditions were: 95°C for 2 min, 35 cycles of 95°C for 45s;

53°C for 45 s; 72°C for 1 min; and a 10 min extension at 72°C. PCR products were size separated on 4% SFR agarose (Amresco) gel.

Statistical analysis

Statistical analyses for phenotypic data were performed using SAS 9.1 (SAS Institute,

2003). A combined analysis was performed for the years 2003 and 2004. Homogeneity of variances across years was evaluated using the Bartlett test (Snedecor and Cochran, 1983). A mixed-effects linear model was used to estimate nymph density values using PROC MIXED in SAS 9.1 (SAS Institute, 2003). The model for both populations was Yijk = µ + Gi + Yj + 50

R(Y)jk + GYij + εijk where µ is the overall mean for nymph density; Gi is the effect of genotypes (i = 1, 2, 3, ..., n, where n is the number of F2:3 lines plus parental lines); Yj is the effect of year (j=1,2); R(Y)jk is the effect of the k replication within each year (k=1,2); GYij is genotype x year interaction, and εijk is the error term. Lines were considered fixed effects and the other effects were considered random. Mean whitefly infestation averaged across all replications and years were calculated with PROC LSMEANS in SAS 9.1 (SAS Institute,

2003). Differences among the parents and the mean of the F2:3 lines were assessed using

Tukey’s multiple-comparison test. The whitefly resistance level of each F2:3 line was compared with the resistance exhibited by the resistant parent. Transgressive segregants for whitefly resistance were defined as F2:3 lines which exceeded the most resistant parent by at least two standard deviations.

Genotypic information was used to construct a linkage map with the polymorphic

SSR markers for each population, using Mapmaker 3.0 (Lander et al., 1987). Linkage was declared with minimum likelihood of odds (LOD) of 3.0 and maximum recombination frequency of 40 cM using the Kosambi's mapping function (Kosambi, 1944).

Genomic regions significantly associated with nymph density in both populations, in

2003, 2004, and in the combined analysis across years, were detected using standard interval mapping (SIM), and composite interval mapping (CIM). QTL analysis was conducted with

R/qtl version 1.14-2 (Broman and Sen, 2009; Broman et al., 2003), an add-on package for the general statistical software R (The R Project for Statistical Computing, 2005). The Haley-

Knott regression methods were used to identify QTL. This method uses the regression of phenotypes on multi-point genotype probabilities (Haley and Knott, 1992). The statistical 51

significance of the results was evaluated by permutation tests. One thousand permutation replicates were used to calculate genome-wide significant LOD thresholds (Lander and

Botstein, 1989). QTL with LOD values higher than the genome-wide threshold were considered significant. A QTL was considered suggestive when the LOD value was larger than two and smaller than the genome-wide threshold value. In the CIM, molecular markers near a large-effect QTL, detected by SIM, were used as covariates, in order to reduce residual variation and to identify other QTL. This analysis uses forward selection to a fixed number of markers, followed by interval mapping, omitting any marker covariates within some fixed distance of the position under test. We used three marker covariates, and a window size on either side of the flanking loci was defined as 20 cM.

RESULTS

Phenotypic analysis

Significant differences (P<0.05) among lines for nymph density were observed in both populations. Given that the effects year and genotype x year interaction were not significant, we combined phenotypic data from 2003 and 2004. In addition, homogeneity of variance for nymph density across years was observed. Both populations showed segregation for the resistance trait (Figures 1-2). Transgressive segregation was observed for both populations. For Williams 79 x Corsoy 79, there were 17 transgressive segregants two standard deviations more resistant than the resistant cultivar, Corsoy 79, in 2004. For population Williams 79 x Cajeme, two transgressive segregants were two standard deviations more resistant than the resistant cultivar, Cajeme, in 2003. 52

In both populations, average nymph density of F2:3 lines showed intermediate values of whitefly susceptibility (Table 1). For Williams 79 x Corsoy 79, the average nymph density value of the F2:3 lines was not significantly different from Williams 79, in both years and in the combined analysis. Conversely, the average nymph density value of F2:3 lines was significantly (P<0.05) different from Corsoy 79, in 2003 and in the combined analysis. For population Williams 79 x Cajeme, in both years, average nymph density across the F2:3 lines was not significantly different from that of the resistant parent. However, the average of the

F2:3 lines was significantly (P<0.05) different from that of the susceptible parent, Williams

79. In both populations, significant (P<0.01) differences between parental lines were observed every year and in the combined analysis, except in 2004 where no significant differences were detected between Williams 79 and Corsoy 79 (Table 1).

Genotypic analysis

For population Williams 79 x Corsoy 79, 181 F2 plants were genotyped with 101 polymorphic SSR markers. In the linkage analysis, 61 markers grouped into 15 linkage groups, while 50 remained unlinked. The length of the linkage map was 725 cM and the average distance between markers was 16.8 cM. For population Williams 79 x Cajeme, 90 polymorphic SSR markers were detected, from these, 60 mapped to 14 linkage groups and 30 markers remained unlinked. The length of the linkage map was 706 cM with an average distance between markers of 16.6 cM. The linkage groups detected in both populations corresponded to soybean chromosomes and the order of the majority of SSR markers coincided with the soybean consensus map (Song et al., 2004). 53

QTL analysis

Eight putative QTL were detected across the two populations (Table 2). QTL detected by SIM also were identified by CIM within the same significant marker interval. Phenotypic variance for whitefly resistance, explained by individual QTL detected in population

Williams 79 x Corsoy 79, ranged from 6 to 13%. In population Williams 79 x Cajeme, QTL for whitefly resistance explained 5 to 18% of the phenotypic variation.

For population Williams 79 x Corsoy 79, the genome-wide 95% LOD thresholds were 3.16, 3.18, and 3.13 for 2003, 2004, and combined analysis, respectively. QTL were detected on chromosomes 10, 12, 16, 17, and 18 (formerly molecular linkage groups O, H, J,

D2, and G respectively) and when combined explained 20.4, 20.2, and 19.9% of the phenotypic variation for 2003, 2004, and combined, respectively (Table 2). A major QTL for whitefly resistance was detected on chromosome 18 every year and in the combined analysis

(Table 2; Figure 3). This QTL was highly significant in 2004 (LOD = 4.5) and in the combined analysis (LOD = 4.88). In 2003 the same QTL was detected as suggestive (LOD between 2 and the genome-wide threshold).

For population Williams 79 x Cajeme, the genome-wide 95% LOD thresholds for

QTL detection were 3.07, 3.28, and 3.19 for 2003, 2004, and combined analysis, respectively. In 2003, QTL for whitefly resistance were detected on chromosomes 7, 12, 18, and 19 (formerly molecular linkage groups M, H, G, and L respectively), and in the combined analysis explained 52.6 % of the phenotypic variation (Table 2). Two significant

QTL were detected on chromosomes 12 (LOD = 3.81) and 19 (LOD = 3.11), and two additional QTL on chromosomes 7 and 18 were declared suggestive (LOD between 2 and the 54

genome-wide threshold) (Table 2; Figure 3). No significant QTL were detected in 2004, however, a suggestive QTL (LOD = 2.82) was detected on chromosome 1 (formerly molecular linkage group D1a), explaining 5% of the phenotypic variation. In the combined analysis, a suggestive QTL (LOD = 2.0) was detected on chromosome 12, explaining 4.9% of the phenotypic variation.

The two SSR markers flanking the most significant QTL region in the two populations were examined for allele-specific effects on phenotype. For population Williams

79 x Corsoy 79, the two SSR markers examined were Satt394 and Satt594. For population

Williams 79 x Cajeme, Satt181 and Satt317 were examined for allele-specific effects on whitefly resistance. For the two populations, consistent with the phenotypes observed in the parental lines, the presence of Williams 79 alleles at the SSR markers correlated with an increase in nymph density (Figures 4-5). Conversely, the presence of alleles from the resistant parents represented a decrease in nymph density. For Williams 79 x Corsoy 79, intermediate phenotypes for whitefly resistance, were observed in individuals that were heterozygous for the two loci, consistent with additive effects of alleles at each locus. On the other hand, heterozygous individuals for the two SSR markers examined in population

Williams 79 x Cajeme showed lower nymph density than the resistant parent.

DISCUSSION

Inheritance to whitefly resistance in soybean was studied in two F2 mapping populations evaluating nymph density on the lower leaf surface. We confirmed that whitefly 55

resistance was a polygenic trait, controlled by QTL on different chromosomes. Xu et al.

(2009) reported that genetic resistance to whitefly was controlled by one major gene and polygenes in soybean. However, these authors did not report specific chromosome locations associated with the resistance trait. In our study, eight putative QTL for whitefly resistance were detected on different soybean chromosomes, across the two segregating populations.

QTL detection varied between populations and years, although whitefly resistance QTL on chromosomes 12 and 18 were detected in more than one year and in both populations. For population Williams 79 x Corsoy 79, whitefly resistance was controlled by a major QTL on chromosome 18 and four minor QTL on chromosomes 10, 12, 16, and 17. These findings are in agreement with the major gene + polygenic mixed inheritance model suggested by Xu et al. (2009). Conversely, population Williams 79 x Cajeme had two major QTL on chromosomes 12 and 19, with three QTL of minor effect on chromosomes 1, 7, and 18.

Most of the significant and suggestive whitefly resistance QTL detected in our study where located on chromosomes that are known to contain several disease and insect resistance genes. Insect resistance in soybean is controlled by two mechanisms, antibiosis, in which feeding causes a disruption of growth and development of the insect; and antixenosis, where the insect is repelled by the host plant (Rector et al., 2000). In our study, we found major whitefly resistance QTL on chromosomes 12 and 18. These two chromosomes are known to have insect-resistance QTL associated with antibiosis and antixenosis (Komatsu et al., 2008; Narvel et al., 2001; Rector et al., 1998, 2000). We also observed one minor whitefly resistance QTL on chromosome 7. Insect resistance QTL associated with antibiosis and antixenosis previously have been mapped to this chromosome (Komatsu et al., 2004; 56

Narvel et al., 2001; Rector et al., 1998, 2000). Epistatic interaction between insect resistance

QTL on chromosomes 7, 12, and 18 also has been documented. Zhu et al. (2006) reported epistatic interactions between a QTL on chromosome 7 controlling antixenosis and other resistance QTL on chromosomes 12 and 18. Similarly, Komatsu et al. (2008) observed interactions between common cutworm-resistance QTL on chromosomes 7 and 18. We detected another significant QTL associated with whitefly resistance on chromosome 19.

Resistance QTL for sclerotinia (Arahana, 2001) and soybean cyst nematode (Guo et al.,

2006) also have been mapped to this chromosome.

In addition to the QTL previously associated with insect resistance, we also detected

QTL on chromosomes 1, which is known to have QTL for pubescence density (Komatsu et al., 2007). Pubescence type and density are associated with whitefly resistance in soybean.

Cultivars with dense pubescence have higher whitefly infestation than glabrous cultivars. On the other hand, cultivars that have trichomes laying flat on the leaves were very resistant to whitefly infestation (Arioglu et al., 1989; Kilen and Lambert, 1993; Lambert et al., 1995;

McAuslane, 1996; McAuslane et al., 1995). Future research will help to elucidate the relationship between the whitefly-resistance QTL on chromosome 1 and the previously reported pubescence density QTL on the same chromosome. Furthermore, it is necessary to clarify the relationship between whitefly resistance of cultivars Corsoy 79 and Cajeme and pubescence density, for a better identification of the resistance mechanisms of whitefly resistance in soybean.

We concluded that the presence of major-effect QTL and several small-effect QTL indicates complex inheritance of whitefly resistance in soybean. Development of insect- 57

resistance cultivars is difficult because plant resistance to insects is most often a quantitatively inherited trait (Kogan and Turnipseed, 1987; McAuslane, 1996). However, identification of molecular markers linked to whitefly-resistance QTL could be used in marker-assisted selection (MAS) programs to develop whitefly-resistance soybean cultivars.

To implement MAS in a breeding program, however, a few QTL with consistent large effects across breeding populations have to be identified (Holland, 2004).

Location of whitefly-resistance QTL in soybean were detected in this study.

Knowledge of the approximate locations of QTL is a starting point for fine mapping or for studying candidate genes that are close to the identified QTL (Bernardo, 2008). Fine- mapping of the large-effect whitefly-resistance QTL detected in this study would be helpful to identify tightly linked markers for MAS programs for cultivar development of whitefly- resistant soybean cultivars.

Acknowledgements

We thank Dr. John Tamulonis for facilitating the genotyping performed by the Monsanto molecular laboratory in Huxley, IA. This is a joint contribution of the Iowa Agriculture and

Home Economics Experiment Station, Ames, IA. Project No. 4403 and the USDA-ARS,

Corn Insects and Crop Genetics Research Unit, and supported by Hatch Act and State of

Iowa. The research was in part financed by funds provided by the Iowa Soybean Association

(ISA), and the North Central Soybean Research Program (NCSRP).

58

REFERENCES

Arahana, V., G.L. Graef, J.E. Specht, J.R. Steadman, and K.M. Eskridge. 2001. Identification

of QTLs for resistance to Sclerotinia sclerotiorum in soybean. Crop Sci. 41:180–188.

Arioglu, H.H. 1987. Screening of some soybean varieties for resistance to whitefly (Bemisia

tabaci Genn.) Soybean Genet. Newsl. 14:136-139.

Arioglu, H.H., A.F Ozgur, and N. Isler. 1989. Influence of soybean pubescence type and

density on whitefly (Bemisia tabaci Genn.) resistance. In: A.J. Pascale (Ed.) Proc.

World Soybean Res. Conf. IV, Argentina, 5-9 March 1989. pp. 1235.

Bernard, R.L. and C.R. Cremeens. 1988a. Registration of 'Corsoy 79' soybean. Crop Sci.

28:1027.

Bernard, R.L. and C.R. Cremeens. 1988b. Registration of 'Williams 79' soybean. Crop Sci.

28:1027.

Bernardo, R. 2008. Molecular markers and selection for complex traits in plants: Learning

from the last 20 years. Crop Sci. 48:1649–1664.

Broman, K.W., H. Wu, S. Sen, and G.A. Churchill. 2003. R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889-890.

Broman, K.W. and S. Sen. 2009. A Guide to QTL Mapping with R/qtl (Springer, Berlin).

Byrne, D.N. and T.S. Bellows. 1991. Whitefly biology. Ann. Rev. Entomol. 36: 431-457. 59

Cregan P.B., T. Jarvik, A.L. Bush, R.C. Shoemaker, K.G. Lark, A.L. Kahler, N. Kaya, T.T.

VanToai, D.G. Lohnes, J. Chung, and J.E. Specht. 1999. An integrated genetic

linkage map of the soybean genome. Crop Sci. 39:1464-1490.

Costa, A.S. 1976. Whitefly-transmitted plant diseases. Ann. Rev. Phytopathol. 16: 429-449.

Faust, R.M. 1992. Conference report and 5-year national research and action plan for

development of management and control methodology for the sweet potato whitefly.

USDA-ARS Publ. 107:165 pp.

Gray, M. 1999. Two-spotted spider mites, whiteflies, and potato leafhoppers: common

inhabitants of many soybean fields. University of Illinois – Integrated Pest

Management Bulletin. 21:192.

Guo, B., D.A. Sleper, P.R. Arelli, J.G. Shannon, and H.T. Nguyen. 2005. Identification of

QTL associated with resistance to soybean cyst nematode races 2, 3, and 5 in soybean

PI 90763. Theor. Appl. Genet. 111:965-971.

Haley, C.S. and S.A. Knott. 1992. A simple regression method for mapping quantitative trait

loci in line crosses using flanking markers. Heredity 69: 315–324.

Holland, J.B. 2004. Implementation of molecular markers for quantitative traits in breeding

programs-challenges and opportunities "New directions for a diverse planet." Proc.

4th Intern. Crop Sci. Congress, Brisbane, Australia, Published on CDROM. Web site

www.cropscience.org.au.

Jones, D. 2003. Plant viruses transmitted by whiteflies. Euro. J. Plant Pathol. 109:195-219. 60

Kilen, T.C. and L. Lambert. 1994. Registration of D89-9121 soybean germplasm line

resistant to foliar feeding insects, phytophthora rot, and stem canker. Crop Sci.

34:1427.

Kogan, M. and S.G. Turnipseed. 1987. Ecology and management of soybean arthropods,

Ann. Rev. Entomol. 32:507-538.

Kogan, M. and D.C. Herzog. 1980. Sampling Methods in Soybean Entomology. Springer,

New York, 587pp.

Kosambi, D. D. 1944. The estimation of map distance from recombination values. Ann.

Eugen. 12:185-199.

Komatsu K., S. Okuda, M. Takahashi, and R. Matsunaga. 2004. Antibiotic effect of insect-

resistant soybean on common cutworm (Spodoptera litura) and its inheritance. Breed.

Sci. 54:27–32.

Komatsu, K., S. Okuda, M. Takahashi, R. Matsunaga, and Y. Nakazawa. 2007. QTL

mapping of pubescence density and flowering time of insect-resistant soybean. Genet.

Mol. Biol. 30:635–639.

Komatsu, K., M. Takahashi, and Y. Nakazawa. 2008. Antibiosis resistance of QTL

introgressive soybean lines to common cutworm (Spodoptera litura Fabricius). Crop

Sci. 48:527–532.

Lambert, A.L., R.M. McPherson, and K.E. Espelie. 1995. Soybean host plant resistance

mechanisms that alter abundance of whiteflies (Homoptera: Aleyrodidae), Environ. 61

Entomol. 24:1381-1386.

Lander, E.S., P. Green, J. Abrahamson, A. Barlow, M.J. Daly, E.S. Lincoln, and L. Newburg.

1987. MAPMAKER: an interactive computer package for constructing primary

genetic linkage maps of experimental and natural populations. Genomics 1: 174–181.

Lander, E. S., and D. Botsein. 1989. Mapping Mendelian factors underlying quantitative

traits using RFLP linkage maps. Genetics 121:185–199.

McAuslane, H.J., F.A. Johnson, D.L. Colvin, and B. Sojack. 1995. Influence of foliar

pubescence on abundance and parasitism of Bemisia argentifolii (Homoptera: Aley-

rodidae) on soybean and peanut. Environ. Entomol. 24:1135-1143.

McAuslane, H. J. 1996. Influence of leaf pubescence on ovipositional preference of Bemisia

argentifolii (Homoptera: Aleyrodidae) on soybean. Environ. Entomol. 25:834-841.

McPherson, R.M. and G.K. Douce, G.K. 1992. Summary of losses from insect damage and

costs of control in Georgia, 1991. Ga. Agric. Expt. Stn. Spec. Publ. 81:65 pp.

Morales, F.J. and P.K. Anderson. 2001. The emergence and dissemination of whitefly-

transmitted geminiviruses in Latin America. Arch. Virol. 146: 415-441.

Narvel, J.M., D.R. Walker, B.G. Rector, J.N. All, W.A. Parrott, and H.R. Boerma. 2001. A

retrospective DNA marker assessment of the development of insect resistant soybean.

Crop Sci. 41:1931–1939.

O’Neal, M. and M.E. Rice. 2006. Soybean aphids look-alikes; don’t be confused by 62

alternatives. Iowa State University. Integrated Crop Management. IC-496(15)

165-166.

Perring, T.M. 2001. The Bemisia tabaci species complex. Crop Prot. 20: 725-737.

R Development Core Team. 2005. R: A language and environment for statistical computing.

R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL

http://www.R-project.org.

Rector, B.G., J.N. All, W.A. Parrott, and H.R. Boerma. 1998. Identification of molecular

markers linked to quantitative trait loci for soybean resistance to corn earworm.

Theor. Appl. Genet. 96:786–790.

Rector, B.G., J.N. All, and H.R. Boerma. 2000. Quantitative trait loci for antibiosis resistance

to corn earworm in soybean. Crop Sci. 40:233–238.

Rice, M.E. 2007. Whiteflies in Iowa soybeans. Iowa State University. Integrated Crop

Management. IC-498(22) pp. 249.

SAS System. SAS Institute Inc. Version 9.1. 2003. Cary, NC.

Schumtz, J., S.B. Cannon, J. Schlueter, J. Ma, T. Mitros, W. Nelson, D.L. Hyten, Q. Song,

J.J. Thelen, J. Cheng, D. Xu, U. Hellsten, G.D. May, Y. Yu, T. Sakurai, T. Umezawa,

M.K. Bhattacharyya, D. Sandhu, B. Valliyodan, E. Lindquist, M. Peto, D. Grant, S.

Shu, D. Goodstein, K. Barry, M. Futrell-Griggs, B. Abernathy, J. Du, Z. Tian, L Zhu,

N. Gill, T. Joshi, M. Libault, A. Sethuraman, X.-C. Zhang, K. Shinozaki, H.T.

Nguyen, R.A. Wing, P. Cregan, J. Specht, J. Grimwood, D. Rokhsar, G. Stacey, R.C. 63

Shoemaker, and S.A. Jackson. 2010. Genome sequence of the paleopolyploid

soybean. Nature 463: 178–183.

Snedecor, G.W. and W.G. Cochran. 1989. Statistical methods. 8th edition. Iowa State

University Press, Ames, pp. 158-160.

Song, Q.J., L.F. Marek, R.C. Shoemaker, K.G. Lark, V.C. ConcibidO, X. Delannay, J.E.

Specht, and P.B. Cregan. 2004. A new integrated genetic linkage map of the soybean.

Theor. Appl. Genet. 109:122–128.

University of Illinois - Integrated Pest Management Bulletin. 2004. Greenhouse whitefly

Trialeurodes vaporariorum.

Vaishampayan, S.M., M. Kogan, G.P. Waldbauer, and J. T. Woolley. 1975. Spectral

specific responses in the visual behavior of the greenhouse whitefly, Trialeurodes

vaporariorum (Homoptera: Aleyrodidae). Ent. Exp. Appl. 18:344–356.

Xu, R., W. Li, L. Zhang, C. Wang, B. Qi, and H. Xing. 2009. Analysis of the inheritance of

resistance to whitefly in soybean. Soybean Sci. 28:954-958.

Zhu, S., D.R. Walker, H.R. Boerma, J.N. All, and W.A. Parrott. 2006. Fine mapping of a

major insect resistance QTL in soybean and its interaction with minor resistance

QTLs. Crop Sci. 46:1094–1099. 64

Table 1. Population mean, standard error, and tests of significance of the difference between genotypic classes, for whitefly infestation in populations Williams 79 x Corsoy 79, and Williams 79 x Cajeme, in 2003, 2004, and the combined analysis of both years.

t test (P-value) of the difference between Mean nymph density† populations

Williams 79 Mean F Mean F F Williams 79 2:3 2:3 2 :3 Resistant† vs. vs. vs. population (susceptible) Resistant Williams 79 Resistant Williams 79 X Corsoy 79

2003 24.06 ± 0.8 30.93 ± 3.9 10.34 ± 4.5 -3.41(0.001) -1.72(0.086) 2.94(0.003)

2004 24.54 ± 1.1 27.83 ± 4.9 22.36 ± 6.3 -0.68(0.499) -0.65(0.519) 0.34(0.736)

Combined 24.29 ± 0.6 29.38 ± 3.8 15.85 ± 3.1 -2.72(0.007) -1.58(0.114) 2.16(0.031)

Williams 79 X Cajeme

2003 24.93 ± 1.1 39.87 ± 5.8 17.55 ± 5.8 -2.69(0.008) -2.50(0.013) 1.24(0.217)

2004 21.39 ± 1.2 35.22 ± 6.2 11.44 ± 6.2 -2.67(0.008) -2.16(0.032) 1.55(0.122)

Combined 23.16 ± 0.8 37.54 ± 4.3 14.50 ± 4.3 -3.78(0.000) -3.28(0.001) 1.97(0.049)

† Nymph density is the number of nymphs in an area of 6.5 cm2 65

Table 2. Summary of significant† and suggestive ‡ whitefly resistance quantitative trait loci (QTL) for populations Williams 79 x Corsoy 79, and Williams 79 x Cajeme, in 2003, 2004, and in the combined analysis across years.

R2 Population Chromosome Position (cM) LOD§ P-value (%) Williams 79 x Corsoy 79

2003 18 22.0 2.03 0.441 5.9 16 4.0 2.33 0.276 6.9 10 17.8 2.33 0.276 7.6 2004 18 23.0 4.50 0.001 12.7 12 17.0 2.42 0.222 7.5 Combined 18 22.0 4.88 0.002 13.3 17 2.0 2.03 0.457 6.6

Williams 79 x Cajeme

2003 12 7 3.81 0.009 18.0 19 90 3.11 0.049 15.6 18 37 2.17 0.338 10.2 7 74 2.14 0.356 8.8 2004 1 Satt129 2.82 0.121 5.03 Combined 12 10 2.00 0.123 4.92

† A QTL was declared significant when the Logarithm of Odds (LOD) value was higher than the genome-wide LOD threshold for significance level of 5%. ‡ A QTL was declared suggestive if LOD was higher than 2 but lower than genome-wide LOD threshold for a significance level of 5%. § For Williams 79 x Corsoy 79, genome-wide LOD thresholds were 3.16, 3.18, and 3.13 for 2003, 2004, and combined analysis, respectively. For Williams 79 x Cajeme, genome-wide LOD thresholds were 3.07, 3.28, and 3.19 for 2003, 2004, and combined analysis, respectively. 66

Number of Lines

Figure 1. Distribution of whitefly infestation among 150 F2:3 lines from the population Williams 79 x Corsoy 79. Histograms show distribution of mean nymph density (number of nymphs in an area of 6.5 cm2) in 2003, 2004, and in the combined analysis of both years. Arrows indicate parental mean values (continuous arrows for Williams 79 and dashed arrows for Corsoy 79).

67

Number of Lines

Figure 2. Distribution of whitefly infestation among 90 F2:3 lines from the population Williams 79 x Cajeme. Histograms show distribution of mean nymph density (number of nymphs in an area of 6.5 cm2) in 2003, 2004, and in the combined analysis of both years. Arrows indicate parental mean values (continuous arrows for Williams 79 and dashed arrows for Cajeme).

68

LOD score

Chromosome

LOD score

Chromosome

Figure 3. Logarithm of Odds (LOD) curves from interval mapping analysis in the combined analysis of two years for Williams 79 x Corsoy 79, and in 2003 for Williams 79 x Cajeme. The dashed horizontal lines indicate the genome-wide 95% LOD threshold. 69

Figure 4. Nymph density (number of nymphs in an area of 6.5 cm2) for different genotypes at the two SSR markers flanking a QTL on chromosome 18 (formerly molecular linkage group G) in the combined analysis. Error bars represent 95% confidence intervals. A: alleles from Williams 79, and B: alleles from Corsoy 79. 70

Figure 5. Nymph density (number of nymphs in an area of 6.5 cm2) for different genotypes at the two SSR markers flanking a QTL on chromosome 12 (formerly molecular linkage group H) in 2003. Error bars represent 95% confidence intervals. A: alleles from Williams 79 and B: alleles from Cajeme. 71

CHAPTER 4: NECTAR COMPOSITION AND FLOWER CHARACTERISTICS OF

WILD PERENNIAL GLYCINE SPECIES1

A paper to be submitted to American Journal of Botany

Paola T. Perez2,7, Anthony H. D. Brown6, Dianne Cook3, Harry T. Horner4, M. Paul Scott2,5,

and Reid G. Palmer2,5

2Department of Agronomy, 3Department of Statistics, 4Department of Genetics, Development and Cell Biology and Microscopy and NanoImaging Facility, Iowa State University, Ames,

IA 50011 USA; 5USDA-ARS-CICGR, Ames, Iowa 50011 USA; 6CSIRO Plant Industry,

GPO Box 1600, Canberra, ACT 2601, Australia

1 Manuscript received ______; revision accepted ______. This is a joint contribution of the Iowa Agriculture and Home Economics Experiment Station, Ames, IA. Project No. 4403 and the USDA-ARS, Corn Insects and Crop Genetics Research Unit, and supported by Hatch Act and State of Iowa. The mention of a trademark or proprietary product does not constitute a guarantee or warranty of the project by Iowa State University or the USDA, and the use of the name by Iowa State University or the USDA, and the use of the name by Iowa State University or the USDA implies no approval of the product to the exclusion of others that may also be suitable. All experiments conducted in this manuscript comply with the current laws of the country in which they were performed.

7 Author for correspondence (e-mail: [email protected]) 72

ABSTRACT

The genus Glycine contains the cultivated annual soybean G. max, the wild annual soybean

G. soja, and about 25 wild perennial Glycine species. Outcrossing by insects occurs at a very low frequency (<1%) in the cultivated species, between <1% to 20% in the wild annual species, and can exceed 50% in certain wild perennial species. Unfortunately, there is a paucity of data on the plant traits that contribute to outcrossing in the wild perennial species.

The objective was to determine nectar composition (sucrose, glucose, fructose, and total carbohydrates) and flower size from greenhouse-grown wild perennial plants. Nectar samples were collected from 70 accessions representing 19 perennial species. Sugars were analyzed by spectrophotometric enzyme assays. Sugar ratios were the most meaningful for comparative purposes. We found variation for nectar sugar composition, and flower size across all the species. Sugar proportions varied across species and on average, the nectar was dominated by sucrose (62.5%), followed by fructose (25.1%) and glucose (12.4%).

Multivariate analysis was used to detect groups of species with similar nectar sugar characteristics. Two groups were detected, one group was formed by the two species G. falcata and G. canescens, which were similar in nectar sugar concentration and composition.

The other 17 species conformed the second group. The information gained from the floral nectar composition of the wild perennial species will guide research for nectar composition- outcrossing associations in the cultivated species. The overall goal is to produce large quantities of F1 hybrid seed for plant breeding studies, and for commercial hybrid seed production.

Key words: flower size; Glycine; nectar sugar composition; soybean wild relatives. 73

INTRODUCTION

The soybean [Glycine max (L.) Merr.] is an economically important crop rich in seed protein (about 40%) and oil (about 20%), and is ranked number one in world oil production in the international trade markets among the major oil seed crops (Wilson, 2004). Given the importance of this crop, breeding efforts have been made to enhance soybean as a food, feed, medicinal, and industrial crop (Palmer, 1999). The gene pools of the soybean, or groups of genes useful for genetic improvement, comprise species of the genus Glycine. This genus is divided into two subgenera Glycine (perennials) and Soja (Moench) F.J. Herm. (annuals).

The progenitor of the soybean is the wild annual species Glycine soja (Sieb. and Zucc.)

(Hymowitz, 2004). Despite the rich diversity in the soybean gene pools, the genetic base of soybean cultivars is extremely narrow. Gizlice et al. (1993) found that fewer than 15 progenitors constituted the major portion of the genetic base for cultivars in U.S soybean production. Soybean breeding programs could exploit the genetic diversity present in exotic germplasm as G. soja and the wild perennial Glycine species. Limited numbers of interspecific crosses between G. max and G. soja have been made in attempts to broaden the genetic base of the cultivated soybean, however, the G. soja accessions harbor many undesirable genetic traits as lodging, vining growth habit, susceptibility to biotic and abiotic stresses, lack of complete leaf abscission, pod shattering, and black seed coat (Palmer and

Hymowitz, 2004). Nevertheless, a successful example of the introgression of G. soja germplasm into the cultivated soybean is the development of small-seeded cultivars, where

G. soja has been used as the donor parent (Fehr et al., 1990a, 1990b; Carter et al., 1995;

Singh and Hymowitz, 1999). 74

The unexploited valuable pool of genetic diversity of the wild perennial species is largely unexplored by soybean breeders for improvement of the cultivated soybean (Singh and Hymowitz, 1999). These species offer a potential source of new germplasm, but embryo rescue techniques are necessary to obtain viable hybrids (Singh et al., 1990, 1993;

Hymowitz, 2004; Palmer and Hymowitz, 2004). The wild perennial Glycine species have genes that could contribute to the improvement of the cultivated soybean.

Hybrids have proven to be a practical method of crop improvement. According to

Palmer et al. (2001), five components are crucial for the successful development of commercial hybrid soybean. The first is that parent combinations must produce heterosis levels superior to the best pure-line cultivars. The second is that a stable male-sterile, female-fertile sterility system is needed. The third is that a selection system to obtain 100% female (pod parent) plants that set seed normally and can be harvested mechanically is necessary. The fourth is that the pollen transfer mechanism has to be efficient in transfer of pollen from the male to pod parent. The fifth requires that hybrid seed production process be profitable for the seed industry and for the farmers. Pollen transfer is the most limiting factor for hybrid soybean production because natural cross-pollination in the cultivated soybean is generally less than 1% (Carlson and Lersten, 2004, Palmer et al., 2009), and manual cross- pollination to produce large quantities of hybrid seed is difficult, time consuming, and expensive. The small size of the soybean flowers, the low success rate and the few seeds obtained per hybrid pod contribute to the difficulty of manually producing large quantities of hybrid seed (Fehr, 1991). Insect cross-pollination of male-sterile soybean plants facilitates the production of hybrid seed (Nelson and Bernard, 1984; Lewers et al., 1996; Ortiz-Perez et 75

al., 2007, Perez et al., 2009a, 2009b). The cultivated soybean attracts relatively few bees

(Blickenstaff and Huggans, 1962); however, bees do visit the flowers for both pollen and nectar.

Glycine species are predominantly self-pollinated, although they possess chasmogamous flowers capable of outcrossing. An example of this phenomenon, two species of the wild perennial Glycine, G. argyrea and G. clandestina, have a dual flowering strategy of both self-fertilized cleistogamous flowers and chasmogamous flowers on the same plant

(Brown et al., 1986; Schoen and Brown, 1991). In a study of selfing in these two species,

Schoen and Brown (1991) showed differences in the levels of selfing due to the contrast in the variability of environmental conditions for insect-mediated pollination between the habitats evaluated. In a study of seed set of chasmogamous and cleistogamous flowers of G. clandestina, Hempel (2004) found that pollinator limitation in chasmogamous flowers is an important factor limiting seed production.

The identification and characterization of sources of insect attraction among the wild perennial species of the genus Glycine could be useful in soybean breeding programs for hybrid production, because in contrast with the low natural out-crossing in the cultivated soybean, out-crossing in the wild perennial soybean species can exceed 50%. We believe that flowers of the perennial Glycine have characteristics that attract and reward insect- pollinators. Insect pollinator attraction is determined by flower color, flower accommodation, flower anthesis and pollen dehiscence, volatile production, and nectary structures and secretion. Nectar is a rich source of sugars and amino acids and provides reward to pollinators, which in turn help to increase the fecundity of those plants that provide nectar 76

(Perret et al., 2001; Carter et al., 2006). The nectar of soybean is a complex mixture of many compounds, consisting of sugars, amino acids, proteins, lipids, and other compounds that provide nutritional and protective functions (Horner et al., 2003). In soybean, honeybees visit the flowers mainly for nectar collection, which has sugar content between 37 and 45%

(Erickson, 1975; Chiari et al., 2005). The quantification and qualification of some of these characteristics, such as the content of glucose, fructose, sucrose (and their ratios), and total carbohydrates would give us an insight into the factors that may affect cross-pollination in the wild perennial soybean.

Thus far, no studies have been done in the search for useful genes for insect attraction in the wild perennial soybean. We performed a comprehensive study of a sample of the wild perennial soybean species to elucidate their variations in nectar composition and in flower characteristics that could be determinants for pollinator-insect attraction and reward in the cultivated soybean. The objective of this study was to determine nectar composition (sucrose, glucose, fructose, and total carbohydrates), and floral morphology from 19 greenhouse- grown wild perennial Glycine species.

MATERIALS AND METHODS

Study plants

Species belonging to subgenus Glycine are perennial and largely indigenous to

Australia and the South Pacific Islands and have been a focus of attention for over two decades because of their relationship to the annual cultivated soybean (Hymowitz, 2004).

The perennial Glycine species have intricate and overlapping combinations of genomes, 77

multiple origins, and ecological variation. In this study, 70 accessions, representing 19 species of the wild perennial soybean were surveyed for nectar composition and flower size

(Table 1).

Study site and methods

Representative accessions of 19 perennial Glycine species were grown in a greenhouse at the Commonwealth Scientific and Industrial Research Organization (CSIRO) in Canberra, Australia from April 2000 to April 2001. The 19 species were accessions from

Australia, Papua New Guinea, Indonesia, the Philippines, and Taiwan (Table 1).

Flowers were collected daily from November 2000 to April 2001, between 9:00 am

(after watering) until 12:00 pm (noon). The anthers were removed before collecting nectar to minimize pollen contamination of the nectar samples. A 1 µL micropipette was inserted at the base of the standard petal and nectar was taken up by capillary action, approximately 0.50

µL of nectar were collected per flower and nectar samples were stored at -20ºC until assayed.

Nectar analysis and flower evaluation

Quantification of glucose, sucrose, fructose, and total carbohydrates were done with

Kits for Carbohydrate Analysis (Sigma-Aldrich). Standard curves for glucose, fructose, sucrose, and total carbohydrates were prepared and used to calculate the concentration of each sample; each sample was measured three times. Assay results were obtained using a 96- well plate reader and absorbance was read using a Schimadzu spectrophotometer. 78

Soybean has typical papilonaceous flowers with a standard petal, two lateral wing petals, and an anterior keel petal (Carlson and Lersten, 2004) (Figure 1). In order to estimate floral structure of the perennial soybean, we recorded length and width of standard, wing, and keel petals. Petal perimeter was estimated as the sum of petal length and width (Dafni and Neal, 1997). In order to simplify the analysis of complete set of petal variables, a principal component analysis (PCA) was performed to reduce the number of variables and summarize the sources of variation. We found most of the variation was explained by the first component, and each petal size had a similar contribution to the variance. For this reason, we concluded it was appropriate to add standard perimeter, wings perimeter, and keel perimeter, to create a new variable called flower size.

Data analysis

All the statistical analyses were done with the general statistical software R (The R

Project for Statistical Computing, 2005). Analyses were performed for three sets of variables.

The first set included sugar (sucrose, glucose, and fructose) and carbohydrates content. The second set was formed by sugar proportions (sucrose, glucose, and fructose percentages) and the variables sugar ratio (sucrose/hexoses). The third set of variables included flower size measurements (petal perimeters and total flower size). Initially, summary statistics were obtained for each of the 19 species, and Pearson correlation coefficients were determined for each pair of variables in a given set. Subsequently, we focused the analyses on species with a sample size equal or larger than ten. Multivariate analysis of variance (MANOVA) for each pair of species was done using nectar information as response variables. Given that we had fewer measurements of flower size per species, and that MANOVA requires more degrees of 79

freedom, we performed a single factor analysis of variance (ANOVA) for flower size. Since the ANOVA was significant for this variable, t-tests were done for each pair of species.

RESULTS

Summary statistics

Sugar contents were estimated for each accession in each species by averaging the results from three measurements per flower. Across the 19 species, the mean value of sucrose was 196 g L-1 (Table 2). Mean hexose sugars values were 44.9 g L-1 and 82.3 g L-1 for glucose and fructose, respectively. Carbohydrates had a mean value of 117.2 g L-1. Pearson correlation coefficients were obtained for each pair of variables (Table 2). A high positive association was observed for fructose and glucose (r=0.9); sucrose content had a similar correlation with glucose and fructose (r=0.6); and carbohydrate content was not associated with sugar contents.

Sugar ratio (sucrose/hexoses) was estimated for each species and the average value across species was 2.2 (Table 3). This result suggested that the nectar composition of the perennial Glycine species was dominated by sucrose. Sugar proportions varied across species and on average, the nectar was dominated by sucrose (62.5%), followed by fructose

(25.1%) and then glucose (12.4%). Negative correlations between sucrose and hexoses percentages were observed (Table 3).

Petal perimeter of the perennial Glycine species was estimated as the sum of petal length and width (Dafni and Neal, 1997). Linear measurements of the floral display could 80

give an accurate estimate of the true area; however, these measurements could be bias depending on the shape and size differences of the flowers under study (Dafni and Neal,

1997). In our study, we used linear measurement estimates only for comparative purposes among the perennial species. On average, the standard petal was larger than the keel and wing petals, although, the perimeter of the two wings combined was larger than the standard petal (Table 4). Total flower size was estimated by adding the measurements of the petals.

Flower size had a high positive correlation with standard and wing petal perimeters (Table

4).

Nectar and flower characteristics of perennial Glycine species

Sample size per species varied from 2 to 134 flowers due to differences in availability of accessions and flowers per species. Sugar and carbohydrates contents showed a large variation among species (Figure 2). Sucrose, glucose, and fructose content varied from 45.0 to 281.3 g L-1, from 4.5 to 76.6 g L-1, and from 22.3 to 149.2 g L-1, respectively (Table 5).

Total carbohydrates per flower varied from 94.5 to 170.9 g L-1. The hybrid G. max x G. tomentella showed the highest content of all sugars, although, only two plants were evaluated. On the other hand, G. pescadrensis had the highest carbohydrate content, with four plants evaluated (Table 5).

Sucrose/hexose ratio showed that the nectar of the perennial Glycine species evaluated was sucrose-dominant, ranging from 1.2 to 8.9 (Table 6, Figure 3). The lowest sucrose/hexoses ratios were observed in G. pullenii, G. falcata, G. pindanica, and G. canescens. Sucrose percentages ranged from 44.7 to 85.5 %, while the ranges of the hexoses were from 4.1 to 32 %, and from 7.7 to 33.4 %, for glucose and fructose, respectively. 81

Considering hexoses only, the proportion of fructose was always higher than that of glucose, except in G. pindanica and G. albicans.

Flower size varied across species (Figure 4). The number of flowers evaluated for this variable was lower than the number of flowers evaluated for nectar composition, and for many species only one accession per species was evaluated. The hybrid G. max x G. tomentella had the largest flowers, however, only two flowers were evaluated (Table 7). G. tabacina and G. falcata also exhibited large flowers and the number of flowers evaluated was

8 and 16, respectively. The smallest flower perimeter was observed in G. aphynota, where only two flowers were evaluated. With a larger sample size (n=38), G. tomentella (diploid) also exhibited small flower size.

Analysis of variance

A multivariate analysis of variance (MANOVA) was performed using information from species where more than 10 flowers were evaluated. Pair-wise comparisons between species were performed using different sets of data. The first multivariate analysis included the response variables: sucrose, glucose, fructose, and carbohydrates content. Wilks lambda test and p-values were reported for each pair of species. Species that were not significantly different were grouped together. Using this approach, two groups were detected; one group was formed by G. falcata and G. canescens, and the other species formed a second group

(Table 8). Another multivariate analysis was performed using the variables sucrose/hexoses ratio, and sucrose, glucose, and fructose percentages. The same group patterns were detected using this approach (Table 9). After the multivariate analysis was done, a single factor analysis of variance (ANOVA) was performed for each sugar nectar variable. The variables 82

that contributed the most to the differences among wild perennial Glycine species were sugar ratio (F9:247= 29.23, P<0.0001), sucrose percentage (F9:247= 28.7, P<0.0001), glucose percentage (F9:247= 14.34, P<0.0001), fructose content (F9:256= 9.7, P<0.0001), sucrose content (F9:247= 8.6, P<0.0001), and total carbohydrates content (F9:256= 2.4, P=0.013).

An ANOVA for flower perimeter was performed with species with more than six flowers evaluated. Species were significantly (F9:69 = 8.05, P < 0.0001) different for flower size. Pair-wise comparisons between species, using t-test, were done. Glycine tomentella

(diploid) had significantly smaller flower size than that of the other species, except for G. lactovirens, where no significant differences were observed (Table 10). Glycine tabacina flowers were larger than those of the other perennial Glycine species. However, these differences were significant only for G. canescens, and G. tomentella (diploid). The second largest flowers were observed for G. falcata. Its flowers were significantly larger than those of G. canescens, and G. tomentella (diploid and tetraploid).

DISCUSSION

Nectar composition and flower size of 19 wild perennial Glycine species were assessed. Variation was found for nectar sugar composition and flower size across all species. Perennial Glycine species have predominantly sucrose-rich nectars. Species with sucrose-rich nectar usually are pollinated by hummingbirds, butterflies, moths, and long- tongued bees; conversely, flowers with hexose-rich nectar are pollinated by short-tongued bees, flies, perching birds, and bats (Baker and Baker, 1983; Freeman et al., 1985; Lammers 83

and Freeman, 1986; Baker et al., 1998; Perret et al., 2001; Dupont et al., 2004; KrÖmer et al.,

2008). To date, no studies on pollinator syndrome of the perennial Glycine have been reported. The cultivated soybean is pollinated by short-tongued bees from the families

Apidae, Megachilidae, Halictidae, Anthophoridae, and Adrenidae (Erickson, 1975; Chang and Kiang, 1987; Ortiz-Perez et al., 2007; Perez et al., 2009). Pollinators from the

Lepidoptera family also have been observed on soybean (Chiari et al., 2005).

In the present study, none of the species evaluated had sucrose/hexoses ratios lower than one. The perennial Glycine species G. falcata had one of the lowest sugar ratios of all the species studied. This species also had larger flowers than those of most of the other species. This could be indicative of a different pollinator syndrome. Based on our current knowledge about phylogenetic relationships in perennial Glycine species, G. falcata belongs to a unique taxonomic/phylogenetic group (Brown et al., 2002; Doyle et al., 2002, 2003,

2004). This species, along with G. canescens, was different in nectar composition from the other Glycine species. Another similarity between these two species is their geographic distribution pattern, both collected from dry interior areas, in contrast to the other species, which were collected mainly in the east coast of Australia (Figure 5). Brown et al. (1990) reported that G. canescens is found in zones of high genetic diversity, which are associated with valleys in the lower slopes of mountain ranges. On the other hand, these two species differed in flower size, G. falcata had significantly larger flowers than those of G. canescens.

The G. tomentella complex was of particular interest given that it is the most widely distributed and diverse of the perennial Glycine species. The G. tomentella complex has overlapping combinations of genomes (Brown et al., 2002; Doyle et al., 2002). We found 84

that both G. tomentella ploidy levels grouped with most of the species under study for nectar composition. However, diploid G. tomentella had significantly smaller flowers than those of the other Glycine species and its tetraploid counterpart.

Several approaches are being explored by us to assess nectar output and composition, nectary structure and development, and increased seed set in the wild perennial Glycine and wild annual Glycine soja species to determine feasibility to introduce them into the cultivated soybean. We report a wide range of variation for nectar composition and flower size across the perennial Glycine species. This variation could be used to improve pollinator-insect attraction in the cultivated soybean. Although, embryo rescue techniques are necessary to obtain viable interspecific hybrids (Singh et al., 1990, 1993; Hymowitz, 2004; Palmer and

Hymowitz, 2004). In the present study, one interspecific hybrid, G. tomentella x G. max, was evaluated. Interestingly, this hybrid had very high sugar contents and large flowers. We could speculate that interspecific hybrids exhibit beneficial insect-pollinator attraction traits.

However, only two plants from this hybrid were evaluated.

Commercial hybrid soybean production could be feasible if insect-pollinator attraction is increased. Perennial Glycine species have a great range of genetic variability for nectar and flower characteristics and eventually, desirable insect-pollinator attraction traits could be introduced into the cultivated soybean. Also, the information gained from the floral nectar composition of the wild perennial species will guide research for nectar composition- outcrossing associations in the cultivated soybean. The overall goal is to produce large quantities of F1 hybrid seed for plant breeding studies, and for commercial hybrid seed 85

production. Unfortunately, there is a paucity of data on the native insect-pollinators of the wild perennial Glycine species.

Acknowledgements

We are most grateful to Dr. Anthony H.D. Brown and to CSIRO Plant Industry, Canberra,

Australia for maintenance of the Glycine species and for their hospitality during Dr. Reid G.

Palmer's study visit.

REFERENCES

Baker, H. G., and I. Baker. 1983. A brief historical review of the chemistry of floral nectar.

In B. Bentley and T. Elias [eds.], The biology of nectaries, 126–152. Columbia

University Press, New York, New York, USA.

Baker, H. G., I. Baker, and S. A. Hodges. 1998. Sugar composition of nectars and fruits

consumed by birds and bats in the tropics and subtropics. Biotropica 30: 559–586.

Blickenstaff, C. C. and J. L. Huggans. 1962. Soybean insects and related arthropods in

Missouri. Missouri Agricultural Experiment Station Research Bulletin 803.

Brown, A. H. D., J. J. Burdon, and J. P. Grace. 1990. Genetic structure of Glycine canescens,

a perennial relative of soybean. Theoretical and Applied Genetics 79: 729-736.

Brown, A. H. D., J. L. Doyle, J. P. Grace, and J. J. Doyle. 2002. Molecular phylogenetic

relationships within and among diploid races of Glycine tomentella (Leguminosae).

Australian Systematic Botany 15: 37– 47. 86

Brown, A. H. D., J. E. Grant, and R. Pullen. 1986. Outcrossing and paternity in Glycine

argyrea by paired fruit analysis. Biological Journal of the Linnean Society 29: 283-

294.

Carlson, J. B., and N. R. Lersten. 2004. Reproductive morphology. In H. R. Boerma and J. E.

Specht [eds.], Soybeans: Improvement, production, and uses, 59–93. Publ. No.16. 3rd

ed. American Society of Agronomy, Madison, Wisconsin, USA.

Carter, T. E., Jr., E. B. Huie, J. W. Burton, F. S. Farmer, and Z. Gizlice. 1995. Registration of

‘Pearl’ soybean. Crop Science 35: 1713.

Carter, C., S. Shafir, L. Vaknin, R. G. Palmer, and R. W. Thornburg. 2006. Proline

accumulating in plant nectars is a determinant of preferential honeybee feeding.

Naturwissenschaften 93: 72-79.

Chang, Y. C., and Y. T. Kiang. 1987. Geometric position of genotypes, honeybee foraging

patterns and out-crossing in soybean. Botanical Bulletin of Academia Sinica 28: 1–11.

Chiari, W. C., V. Alencar, M. C. Colla, V. M. Attencia, F. M. Costa, C. S. Kotaka, E.

Shiguero, and H. R. Magalhaes. 2005. Floral biology and behavior of Africanized

honeybees Apis mellifera in soybean (Glycine max L. Merrill). Brazilian Archives of

Biology and Technology 48: 367-378.

Dafni, A. and P. Neal. 1997. Size and shape in floral advertisement: measurement, concepts

and implications. Acta Horticulturae 437: 121-140.

Doyle, M. J., and A. H. D. Brown. 1985. Numerical analysis of isozyme variation in Glycine

tomentella. Biochemical Systematics and Ecology 13: 413-419. 87

Doyle, J. J., J. L. Doyle, A. H. D. Brown, and R. G. Palmer. 2002. Genomes, multiple

origins, and lineage recombination in the Glycine tomentella (Leguminosae)

polyploid complex: Histone H3-D gene sequences. Evolution 56: 1388-1402.

Doyle, J. J., J. L. Doyle, and C. Harbison. 2003. Chloroplast-expressed glutamine synthetase

in Glycine and related Leguminosae: phylogeny, gene duplication, and ancient

polyploidy. Systematic Botany 28: 567-577.

Doyle, J. J., J. L. Doyle, J. T. Rauscher, and A. H. D. Brown. 2004. Evolution of the

perennial soybean polyploid complex (Glycine subgenus Glycine): a study of

contrasts. Biological Journal of the Linnean Society 82: 583-597.

Dupont, Y. L., D. M. Hansen, J. T. Rasmussen, and J. M. Olesen. 2004. Evolutionary

changes in nectar sugar composition associated with switches between bird and insect

pollination: the Canarian bird-flower element revisited. Functional Ecology 18: 670–

676.

Erickson, E. H. 1975. Variability of floral characteristics influences honey bee visitation to

soybean blossoms. Crop Science 15: 767–771.

Fehr, W. R. 1991. Principles of cultivar development. Theory and technique. Macmillan

Publishing Company, Ames, Iowa, U.S.A.

Fehr, W. R., S. R. Cianzio, and G. A. Welke. 1990a. Registration of ‘SS201’soybean. Crop

Science 30: 1361.

Fehr, W. R., S. R. Cianzio, G. A. Welke, and A. R. LeRoy. 1990b. Registration of SS202’

soybean. Crop Science 30: 1361.

Freeman, C. E., W. H. Reid, and R. D. Worthington. 1985. Patterns of floral nectar-sugar

composition of Ipomopsis longiflora (Polemoniaceae) near the contact zone of its 88

subspecies longiflora and australis. American Journal of Botany 72: 1662–1667.

Gizlice, Z., T. E. Carter, Jr., and J. W. Burton. 1993. Genetic diversity in North American

soybean: I. Multivariate analysis of founding stock and relation to coefficient of

parentage. Crop Science 33: 614–620.

Hempel, K.A. 2004. Advantages of chasmogamy and cleistogamy in a perennial Glycine

clandestina Wendl. (). Ph.D. Diss., Australian National University, Canberra,

Australia.

Horner, H. T., R. A. Healy, T. Cervantes-Martinez, and R. G. Palmer. 2003. Floral nectary

fine structure and development in Glycine max L. (Fabaceae). International Journal

of Plant Science 164: 675–690.

Hymowitz, T. 2004. Speciation and cytogenetics. In H. R. Boerma and J. E. Specht [eds.],

Soybeans: Improvement, production, and uses, 59–93. Publ. No.16. 3rd ed. American

Society of Agronomy, Madison, Wisconsin, USA.

KrÖmer, T., M. Kessler, G. Lohaus, and A. N. Schmidt-Lebuhn. 2008. Nectar sugar

composition and concentration in relation to pollinator syndrome in Bromeliaceae.

Plant Biology 10: 502-511.

Lammers, T. G., and C. E. Freeman. 1986. Ornithophily among the Hawaiian Lobelioideae

(Campanulaceae): evidence from floral nectar sugar composition. American Journal of

Botany 73: 1613-1619.

Lewers, K. S. 1996. Production, evaluation, and utilization of hybrid soybean [Glycine max

(L.) Merr.], Ph.D. Diss., Iowa State University, Ames (Diss. Abstr. 96-26046).

Nelson, R. L., and Bernard, R. L. 1984. Production and performance of hybrid soybeans.

Crop Science 24: 549-553. 89

Ortiz-Perez, E., S. R. Cianzio, H. Wiley, H. T. Horner, W. H. Davis, and R. G. Palmer. 2007.

Insect-mediated cross-pollination in soybean [Glycine max (L.) Merr.]: I. Agronomic

performance. Field Crops Research 101: 259-268.

Palmer, R. G. 1999. Germplasm and crop improvement for breaking yield barriers: Examples

from soybean. In R. K. Behl, M. S. Punia, and B. P. S. Lather [eds.], Crop

improvement for food security, 72-115. SSARM, Hisar, India.

Palmer, R. G., J. Gai, H. Sun, and J. W. Burton. 2001. Production and evaluation of hybrid

soybean. Plant Breeding Reviews 21: 263-307.

Palmer, R. G. and T. Hymowitz. 2004. Soybean Genetics and Breeding. In C. Wrigley, H.

Corke, and C. Walker [eds.], Encylopedia of grain sciences, 136-146. Academic

Press, London, England.

Palmer, R. G., P. T. Perez, E. Ortiz-Perez, F. Maalouf, and M. J. Suso. 2009. The role of

crop-pollinator relationships in breeding for pollinator-friendly legumes: from a

breeding perspective. Euphytica 170: 35-52.

Perez. P. T., S. R. Cianzio, and R. G. Palmer. 2009. Evaluation of soybean [Glycine max (L.)

Merr.] F1 hybrids. Journal of Crop Improvement 23: 1-18

Perret M., A. Chautems, R. Spichiger, M. Peixoto, and V. Savolainen. 2001. Nectar sugar

composition in relation to pollination syndromes in Sinningieae (Gesnerianceae).

Annals of Botany 87: 267-273.

R Development Core Team. 2005. R: A language and environment for statistical computing.

R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL

http://www.R-project.org. 90

Schoen, D, and A. H. D. Brown. 1991. Whole- and part-flower self-pollination in Glycine

clandestina and G. argyrea and the evolution of autogamy. Evolution 45: 1651-1664.

Singh, R. J., K. P. Kollipara, and T. Hymowitz. 1990. Backcross-derived progeny from

soybean and Glycine tomentella Hayata intersubgeneric hybrids. Crop Science 30:

871–874.

Singh, R. J., K. P. Kollipara, and T. Hymowitz. 1993. Backcross (BC2–BC4)-derived fertile

plants from Glycine max and G. tomentella intersubgeneric hybrids. Crop Science 33:

1002–1007.

Singh, R. J. and T. Hymowitz. 1999. Soybean genetic resources and crop improvement.

Genome 42: 605-616.

Wilson, R.F. 2004. Seed composition. In H. R. Boerma and J. E. Specht [eds.], Soybeans:

Improvement, production, and uses, 621–677. Publ. No.16. 3rd ed. American Society

of Agronomy, Madison, Wisconsin, USA. 91

Table 1. Species evaluated of the genus Glycine Willd., subgenus Glycine, and the interspecific hybrid G. max x G. tomentella, ploidy level, genome type, accession number, and geographic distribution*.

Geographic Subgenus Glycine 2n Genome Accession number** distribution

G. albicans Tind. and Craven 40 II G 2049, G 2937 Australia G. aphyonota B.E. Pfeil 40 ? G 2589 Australia G. arenaria Tind. 40 HH G 3136 Australia G. argyrea Tind. 40 A2A2 G 1420, G 1626 Australia G 1232, G 1301, G 1302, G G. canescens F.J. Herman 40 AA 1485, G 3004, G 3014, G 3088, Australia G 3094 G 1126, G 1171, G 2347, G G. clandestina Wendl. 40 A1A1 Australia 2860, G 3114 G. curvata Tind. 40 C1C1 G 1396 Australia G 1184, G 1236, G 1842, G G. cyrtoloba Tind. 40 CC Australia 3127 G 1824, G 1825, G 2087, G G. falcata Benth. 40 FF 2088, G 2440, G 2658, G 3019, Australia G 3079, G 3105 G. lactovirens Tind. and Craven 40 I1I1 G 2597, G 2598, G 2720 Australia G 1137, G 1456, G 1580, G G. latifolia (Benth.) Newell and 2037, G 2115, G 2139, G 2253, 40 B1B1 Australia Hymowitz G 2287, G 2345, G 2756, G 2785, G 2815, G 3082, G 3121 G. max x G. tomentella 40 G 113519 G 1143, G 1477, G 1703, G G. microphylla (Benth.) Tind. 40 BB Australia 1830, G 2704 G. pescadrensis G 1142, G 1298, G 3123 G. pindanica Tind. and Craven 40 H2H2 G 2951 Australia G. pullenii B.E Pfeil, Tind. and ? ? G 2728, G 2786 Australia Craven G. stenophita B.E. Pfeil and 40 B3B3 G 2600 Australia Tind.

Australia, West B2B2 - G 1075, G 1314, G 1706, G G. tabacina (Labill.)Benth. 40-80 Central and South Complex 2475, G 2761, G 2875, G 3116, Pacific Islands G 1300, G 1303, G 1316, G 1366, G 1413, G 1749, G 1932, G 1941, G 1944, G 1946, G Australia, Papua G. tomentella Hayata 40 DD 2053, G 2055, G 2076, G 2723, New Guinea G 3098, G 3101, G 3108, G 3119, G 3128, G 3129 Australia, Papua G 1133, G 1134, G 1398, G New Guinea, 1487, G 1763, G 1929, G 1945, 78 - 80 Complex Indonesia, G 2098, G 2437, G 2469, G Philippines, 2474, G 2491, G 2557 Taiwan * Adapted from Hymowitz (2004) ** Accession number, CSIRO Plant Industry, Canberra, Australia. 92

Table 2. Nectar sugar content and correlations across perennial Glycine species and the interspecific hybrid G. max x G. tomentella.

Correlation Variable Mean s.d Sucrose Glucose Fructose Carbohydrates Sucrose (g L-1) 196 76.7 1 0.6 0.59 0.31 Glucose (g L-1) 44.9 30.3 1 0.90 0.15 Fructose (g L-1) 82.3 50 1 0.10 Carbohydrates (g L-1) 117.2 48.1 1

93

Table 3. Sugar ratio, sugar percentages, and correlations across perennial Glycine species and the interspecific hybrid G. max x G. tomentella,

Correlation Variable Mean s.d. Sugar ratio Sucrose Glucose Fructose Sugar ratio 2.2 1.9 1 0.81 -0.61 -0.75 Sucrose (%) 62.5 12.1 1 -0.79 -0.90 Glucose (%) 12.4 5.9 1 0.44 Fructose (%) 25.1 8.2 1

94

Table 4. Flower size across species and correlations across perennial Glycine species and the interspecific hybrid G. max x G. tomentella.

Correlation Mean Variable s.d. Flower Standard Wings Keel (mm) size perimeter perimeter perimeter Flower size 3.9 0.6 1 0.94 0.96 0.71 Standard perimeter 1.5 0.2 1 0.85 0.57 Wings perimeter 1.7 0.3 1 0.59 Keel perimeter 0.7 0.1 1

95

Table 5. Nectar sugar content [mean (s.d.)] of perennial Glycine species and the interspecific hybrid G. max x G. tomentella.

Sucrose Glucose Fructose Carbohydrates Species n (g L-1) (g L-1) (g L-1) (g L-1) G. max x G. tomentella 2 281.3 (65) 76.6 (33.7) 149.2 (67.7) 146.1 (38.6) G. arenaria 2 252.1 (11.9) 53.8 (3.0) 104.4 (6.9) 122.4 (18.2) G. aphyonota 4 242.2 (31.7) 66.3 (19.7) 126.0 (27.1) 95.1 (18.4) G. curvata 4 231.5 (48.7) 27.0 (19.0) 50.1 (18.8) 137.9 (13.9) G. clandestina 15 223.6 (55.4) 41.3 (26.4) 69.8 (33.4) 120.4 (44.3) G. albicans 7 219.8 (75.6) 27.5 (39.4) 22.9 (16.9) 129.4 (29.3) G. cyrtoloba 12 215.2 (91.1) 19.1 (15.1) 37.0 (11.4) 140.4 (74.1) G. canescens 67 204.6 (75.9) 63.7 (32.1) 116.6 (51.4) 119.2 (29.9) G. latifolia 33 196.6 (77.5) 28.0 (20.4) 62.0 (44.7) 120.4 (33.6) G. lactovirens 9 195.3 (62.3) 9.3 (4.3) 22.3 (11.0) 121.2 (35.2) G. tomentella (diploid) 55 193.1 (78.8) 50.5 (25.7) 92.1 (47.1) 118.6 (31.8) G. argyrea 8 192.7 (40.1) 39.3 (15.8) 73.1 (18.8) 102.3 (36.1) G. pescadrensis 4 181.7 (148.6) 40.7 (40.5) 54.5 (36.3) 170.9 (75.6) G. microphylla 12 177.5 (57.9) 31.9 (35.6) 62.6 (41.7) 131.4 (31.3) G. tomentella (tetraploid) 31 170.5 (68.9) 39.3 (25.4) 71.3 (39.3) 121.9 (31.6) G. tabacina 15 169.8 (78.1) 34.3 (28.3) 62.0 (46.2) 122.8 (35.8) G. pullenii 5 163.8 (34.7) 47.9 (9.7) 88.1 (16.2) 106.1 (39.4) G. falcata 17 161.5 (86.1) 48.3 (29.8) 93.8 (47.5) 116.2 (36.5) G. stenophita 1 71.6 4.5 27.7 94.5 G. pindanica 2 45.0 (35.3) 29.0 (30.8) 22.6 (6.8) 117.4 (23)

96

Table 6. Sugar ratios and sugar percentages [mean (s.d.)] of perennial Glycine species and the interspecific hybrid G. max x G. tomentella.

Fructose Species n S/H* Sucrose (%) Glucose (%) (%) G. pullenii 5 1.2 (0.3) 54.5 (6.1) 16.0 (2.4) 29.5 (3.8) G. falcata 17 1.2 (0.4) 52.8 (6.3) 13.8 (5.6) 33.4 (8.2) G. pindanica 2 1.2 (1.2) 44.7 (31.3) 32.1 (35.6) 23.1 (4.4) G. canescens 67 1.24 (0.5) 53.67 (7.6) 15.79 (3.9) 30.54 (6.0) G. aphyonota 4 1.3 (0.2) 56.1 (2.9) 15.0 (1.9) 28.9 (2.0) G. max x G. tomentella 2 1.3 (0.3) 56.4 (5.7) 14.8 (1.8) 28.8 (3.9) G. arenaria 2 1.6 (0.2) 61.5 (2.6) 13.1 (0.8) 25.4 (1.8) G. tomentella (diploid) 51 1.6 (0.9) 58.8 (8.2) 13.4 (3.9) 27.9 (5.5) G. argyrea 8 1.8 (0.4) 63.6 (5.3) 12.5 (2.9) 23.9 (3.0) G. tomentella 29 1.8 (0.6) 62.0 (8.8) 12.3 (5.3) 25.6 (6.9) (tetraploid) G. stenophita 1 2.2 68.9 4.4 26.7 G. tabacina 13 2.2 (0.7) 67.3 (8.4) 10.2 (4.0) 22.5 (4.7) G. pescadrensis 4 2.3 (1.0) 67.4 (8.3) 11.6 (6.4) 21.1 (4.8) G. clandestina 15 2.4 (1.0) 68.0 (9.2) 11.4 (4.5) 20.6 (6.1) G. latifolia 33 2.8 (1.3) 71.1 (71.1) 8.6 (3.9) 20.3 (5.9) G. microphylla 12 3.2 (2.4) 69.2 (14.7) 9.2 (7.6) 21.6 (9.4) G. curvata 4 3.4 (1.3) 75.5 (9.1) 8.4 (4.9) 16.1 (4.2) G. cyrtoloba 12 4.2 (2.2) 78.1 (8.0) 6.7 (2.5) 15.2 (6.8) G. lactovirens 8 6.7 (2.9) 85.4 (5.3) 4.1 (1.6) 10.6 (4.1) G. albicans 7 8.9 (6.3) 83.6 (13.9) 8.6 (11.4) 7.7 (4.5)

* S/H = [%Sucrose]/[%Glucose + %Fructose] 97

Table 7. Flower size and petal perimeter [mean (s.d.)] of perennial Glycine species and the interspecific hybrid G. max x G. tomentella.

Standard Wings Keel Flower size Species n perimeter perimeter perimeter (mm) (mm) (mm) (mm) G. max x G. tomentella 2 5.0 1.8 2.4 0.8 G. tabacina 8 4.7 (0.5) 1.8 (0.2) 2.2 (0.3) 0.7 (0.0) G. pescadrensis 2 4.7 1.8 2.2 0.7 G. stenophita 2 4.6 1.9 2.1 0.7 G. falcata 16 4.4 (0.4) 1.6 (0.2) 2.0 (0.2) 0.7 (0.1) G. cyrtoloba 8 4.3 (0.5) 1.6 (0.2) 2.0 (0.3) 0.7 (0.0) G. clandestina 8 4.1 (0.6) 1.5 (0.3) 1.9 (0.2) 0.7 (0.1) G. microphylla 10 4.1 (0.2) 1.6 (0.0) 1.8 (0.2) 0.7 (0.1) G. tomentella 24 4.0 (0.5) 1.5 (0.2) 1.8 (0.2) 0.7 (0.1) (tetraploid) G. argyrea 4 4.0 (0.4) 1.7 (0.2) 1.6 (0.2) 0.7 (0.0) G. latifolia 24 4.0 (0.3) 1.5 (0.2) 1.8 (0.2) 0.7 (0.1) G. albicans 4 3.8 (0.5) 1.3 (0.1) 1.6 (0.3) 0.8 (0.1) G. canescens 16 3.8 (0.2) 1.5 (0.1) 1.7 (0.2) 0.6 (0.0) G. lactovirens 6 3.7 (0.6) 1.3 (0.3) 1.6 (0.2) 0.8 (0.1) G. curvata 2 3.6 1.4 1.5 0.6 G. pullenii 4 3.5 (0.1) 1.4 (0.0) 1.6 (0.1) 0.6 (0.1) G. tomentella (diploid) 38 3.2 (3.2) 1.2 (1.2) 1.4 (1.4) 0.6 (0.6) G. pindanica 2 3.1 1.2 1.4 0.5 G. arenaria 2 3.1 1.3 1.3 0.5 G. aphyonota 2 2.7 1.0 1.1 0.5

98

Table 8. Multivariate analysis of variance (MANOVA) for nectar sugar content, using species with a sample size larger than ten. Wilk's lambda test and p-value (in parentheses) are reported. Wilk's lambda test and p-value of a pair of species that are non-significantly different in nectar sugar content are reported in bold.

tomentella tomentella Glycine falcata latifolia clandestina tabacina microphylla cyrtoloba lactovirens (diploid) (tetraploid) canescens 0.85 0.82 0.95 0.55 0.65 0.78 0.74 0.53 0.52 (0.00) (0.00) (0.40) (0.00) (0.00) (0.00) (0.00) (0.00) (0.00) tomentella 0.92 0.83 0.70 0.75 0.91 0.77 0.54 0.56 (diploid) (0.17) (0.09) (0.00) (0.00) (0.23) (0.00) (0.00) (0.00) tomentella 0.81 0.81 0.78 0.94 0.81 0.55 0.52 (tetraploid) (0.07) (0.01) (0.04) (0.70) (0.09) (0.00) (0.00) falcata 0.44 0.52 0.55 0.67 0.27 0.24 (0.00) (0.00) (0.00) (0.04) (0.00) (0.00) latifolia 0.89 0.94 0.91 0.75 0.66 (0.27) (0.63) (0.41) (0.02) (0.00) clandestina 0.86 0.84 0.62 0.58 (0.46) (0.40) (0.03) (0.04) tabacina 0.85 0.63 0.44 (0.51) (0.05) (0.01) microphylla 0.73 0.62 (0.18) (0.10) cyrtoloba 0.71 (0.25)

99

Table 9. Multivariate analysis of variance (MANOVA) for nectar sugar ratio and percentages, using species with a sample size larger than ten. Wilk's lambda test and p-value (in parentheses) are reported. Wilk's lambda test and p-value of a pair of species with non- significant differences in nectar sugar proportions are reported in bold.

tomentella tomentella Glycine falcata latifolia clandestina tabacina microphylla cyrtoloba lactovirens (diploid) (tetraploid) canescens 0.88 0.79 0.94 0.46 0.66 0.69 0.68 0.39 0.25 (0.00) (0.00) (0.15) (0.00) (0.00) (0.00) (0.00) (0.00) (0.00) tomentella (diploid) 0.94 0.85 0.62 0.73 0.80 0.82 0.53 0.35 (0.20) (0.01) (0.00) (0.00) (0.00) (0.01) (0.00) (0.00) tomentella (tetraploid) 0.74 0.77 0.82 0.88 0.81 0.53 0.30 (0.01) (0.00) (0.05) (0.16) (0.05 (0.00) (0.00) falcata 0.37 0.48 0.49 0.61 0.20 0.12 (0.00) (0.00) (0.00) (0.01) (0.00) (0.00) latifolia 0.88 0.95 0.86 0.84 0.52 (0.14) (0.52) (0.10) (0.07) (0.00) clandestina 0.84 0.83 0.68 0.39 (0.25) (0.22) (0.03) (0.00) tabacina 0.86 0.67 0.33 (0.34) (0.03) (0.00) microphylla 0.83 0.63 (0.29) (0.05) cyrtoloba 0.70 (0.12)

100

Table 10. Pair wise comparisons of flower size, using species with sample size three or larger. Student t-test and p-values (in parentheses) are reported. t-test and p-value of a pair of species with non-significant differences in flower size are reported in bold.

tomentella tomentella Glycine falcata latifolia clandestina tabacina microphylla cyrtoloba lactovirens (diploid) (tetraploid) canescens 4.00 -1.01 -3.51 -2.11 -1.23 -3.10 -2.23 -1.67 2.2 (0.00) (0.31) (0.00) (0.10) (0.31) (0.00) (0.14) (0.20) (0.90) tomentella (diploid) -4.12 -6.00 -5.61 -3.02 -4.89 -5.51 -3.56 -1.31 (0.00) (0.00) (0.00) (0.00) (0.00) (0.00) (0.00) (0.32) tomentella (tetraploid) -2.22 -0.51 -0.56 -2.32 -0.74 -1.01 0.62 (0.00) (0.61) (0.62) (0.11) (0.46) (0.29) (0.56) falcata 2.12 0.89 -0.67 1.78 0.40 1.67 (0.13) (0.41) (0.45) (0.11) (0.67) (0.21) latifolia -0.31 -2.22 -0.30 -0.78 0.78 (0.78) (0.10) (0.78) (0.45) (0.52) clandestina -1.33 0.21 -0.42 0.89 (0.23) (0.89) (0.67) (0.44) tabacina 2.01 1.02 2.01 (0.12) (0.44) (0.12) microphylla -0.67 0.89 (0.52) (0.45) cyrtoloba 1.20 (0.33)

101

Figure 1. Typical soybean flower. A) Standard petal. B) Two wing petals. C) Fused keel petals.

A

B C B

102

Figure 2. Distribution of nectar sugar content of perennial Glycine species* with a sample size larger than ten. Mean value represented by X symbol.

* G. tomentella_D represents diploid species of Glycine tomentella and G. tomentella_T represents tetraploid species of Glycine tomentella 103

Figure 3. Distribution of nectar sugar ratio and sugar percentages of perennial Glycine species* with a sample size larger than ten. Mean value represented by X symbol.

* G. tomentella_D represents diploid species of Glycine tomentella and G. tomentella_T represents tetraploid species of Glycine tomentella 104

Figure 4. Distribution of petal perimeter and flower size of perennial Glycine species* with a sample size larger than ten. Mean value represented by X symbol.

* G. tomentella_D represents diploid species of Glycine tomentella and G. tomentella_T represents tetraploid species of Glycine tomentella 105

Figure 5. Geographic distribution of G. clandestina, G. cyrtoloba, G. latifolia, G. microphylla, G. tomentella, G. tabacina, G. canescence, G. falcata, and G. lactovirens.

106

CHAPTER 5: GENERAL CONCLUSIONS

The purpose of this research project was to investigate different strategies for soybean improvement, from genetic resistance to biotic stresses, to surveying perennial wild relatives for insect-pollinator attraction characteristics that could be useful in hybrid seed production in the cultivated species.

The first study evaluated four plant introductions from south-central China to identify new sources of resistance to BSR. Four F2 segregating populations were developed and BSR resistance/susceptibility was evaluated under growth chamber conditions. Using single marker regression with SSR markers linked to previously identified BSR resistance genes, high association was found between the SSR markers used and BSR resistance in populations

PI 594638B X Century 84, PI 594650A X Century 84, and PI 594658B X Century 84. It was concluded that the BSR resistance genes present in these PIs could be allelic to previously identified BSR resistance genes on chromosome 16. Population PI 594637 X Century 84 showed non-association with SSR Satt547. This could be indicative of the presence of a new

BSR resistance gene in PI 594637. In the same population, highly resistant F2:3 individuals were present. This highly resistant progeny can be used for development of BSR resistant cultivars. In order to test the hypothesis of a different BSR resistance gene present in PI

594637, allelism tests between PI 594637 and genotypes possessing Rsb1, Rbs2, and Rbs3 genes need to be performed.

The objective of the second study was to investigate the genetic factors associated with whitefly resistance in soybean. Two F2 segregating populations for whitefly resistance 107

were developed for QTL analysis. It was found that whitefly resistance was a polygenic trait, controlled by QTL on different chromosomes. Eight putative whitefly-resistance QTL were detected on different soybean chromosomes across the two populations. QTL detection varied between populations and years, although whitefly-resistance QTL on chromosomes 12 and 18 were detected in more than one year and in both populations. Most of the whitefly- resistance QTL detected in the study were located on chromosomes that are known to contain several disease and insect resistance genes. The presence of major-effect QTL and several small-effect QTL indicated complex inheritance of whitefly resistance in soybean. Future research should focus on fine-mapping of the large-effect whitefly-resistance QTL detected in this study. This will help to identify tightly linked markers to whitefly-resistance QTL, which can be use for marker-assisted selection.

The third study surveyed the floral characteristics of wild perennial Glycine species.

The nectar sugar composition and flower size varied among the 19 species evaluated. The study found that perennial Glycine species have predominantly sucrose-rich nectars. These nectars are usually attractive to pollinators such as hummingbirds, butterflies, moths, and long-tongued bees. Studies of pollinator syndrome of wild perennial Glycine species are needed. In this study, multivariate analysis of nectar sugar characteristics identified one group consisting of G. falcata and G. canescens, which had different nectar sugar characteristics in comparison with the other perennial species evaluated. These two species have similar geographic distribution, both collected from dry interior areas. The information gained from the floral nectar composition and flower morphology of the wild perennial species will guide research for nectar composition-outcrossing associations in the cultivated 108

soybean. Future research should focus on the native insect-pollinators of the wild perennial

Glycine species.